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INTRODUCTION 


ICASE*  is  operated  at  the  Langley  Research  Center  (LaRC)  of  NASA  by  the  Universities  Space  Research 
Association  (USRA)  under  a  contract  with  the  Center.  USRA  is  a  nonprofit  consortium  of  major  U.S.  colleges 
and  universities. 

The  Institute  conducts  unclassified  basic  research  in  applied  mathematics,  numerical  analysis  and  algo¬ 
rithm  development,  computer  science,  fluid  mechanics,  and  structures  and  materials  in  order  to  extend  and 
improve  problem-solving  capabilities  in  science  and  engineering,  particularly  in  the  areas  of  aeronautics  and 
space  research. 

ICASE  has  a  small  permanent  staff.  Research  is  conducted  primarily  by  its  permanent  staff  and  visiting 
scientists  from  universities  and  industry  who  have  resident  appointments  for  limited  periods  of  time  as  well 
as  by  visiting  and  resident  consultants.  Members  of  NASA’s  research  staff  may  also  be  residents  at  ICASE 
for  limited  periods. 

The  major  categories  of  the  current  ICASE  research  program  are: 

•  Applied  and  numerical  mathematics,  including  multidisciplinary  design  optimization; 

•  Applied  computer  science:  system  software,  systems  engineering,  and  parallel  algorithms; 

•  Theoretical,  computational,  and  experimental  research  in  fluid  mechanics  in  selected  areas  of  interest 
to  LaRC,  such  as  transition,  turbulence,  flow  control,  and  acoustics;  and 

•  Theoretical,  computational,  and  experimental  research  in  structures  and  material  sciences  with 
emphasis  on  smart  materials  and  nanotechnologies. 

ICASE  reports  are  primarily  considered  to  be  preprints  of  manuscripts  that  have  been  submitted  to 
appropriate  research  journals  or  that  are  to  appear  in  conference  proceedings.  A  list  of  these  reports  for 
the  period  April  1,  2000  through  September  30,  2000  is  given  in  the  Reports  and  Abstracts  section,  which 
follows  a  brief  description  of  the  research  in  progress. 


*ICASE  is  operated  at  NASA  Langley  Research  Center,  Hampton,  VA,  under  the  National  Aeronautics  and  Space  Adminis¬ 
tration,  NASA  Contract  No.  NAS1-97046.  Financial  support  was  provided  by  NASA  Contract  Nos.  NASl-97046,  NAS1-19480, 
NAS1-18605,  NASl-18107,  NAS1-17070,  NAS1-17130,  NAS1-15810,  NAS1-16394,  NAS1-14101,  and  NAS1-14472. 
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RESEARCH  IN  PROGRESS 


APPLIED  AND  NUMERICAL  MATHEMATICS 


Active  flow  control  research  tool  development 

Brian  G.  Allan 

The  use  of  secondary  flow  control  devices  on  inlets,  such  as  vortex  generators  and  small  synthetic  jets, 
has  recently  attracted  interest.  Significant  improvements  in  engine  flow  distortion  and  pressure  loss  have 
been  demonstrated  using  these  small  flow  control  devices.  This  technology  has  shown  a  great  potential  for 
improving  existing  inlet  configurations  and  enabling  more  advanced  inlet  designs,  allowing  for  reduced  inlet 
lengths  and  sharper  turning  angles.  In  order  to  take  advantage  of  these  devices,  computational  tools  need 
to  be  developed  in  order  to  evaluate  potential  engine  inlet  designs.  The  goal  of  this  research  is  to  develop 
computational  boundary  conditions  and  techniques  for  a  CFD  code,  in  order  to  model  vortex  generators  and 
synthetic  jets  for  an  inlet  application. 

Using  a  compressible  Navier-Stokes  code,  developed  at  NASA,  numerical  simulations  of  micro-vortex 
generators  and  synthetic  jets  on  a  flat  plate  and  inlet  duct  flows,  will  be  performed.  These  simulations  are  to 
be  performed  in  conjunction  with  wind  tunnel  experiments  conducted  at  NASA  Langley  Research  Center. 
From  these  simulations  and  experiments,  a  better  understanding  of  the  flow  physics  of  the  flow  control  devices 
will  be  gained.  Using  these  insights  into  the  flow  physics  of  the  control  devices,  a  computational  boundary 
condition  can  then  be  developed.  This  boundary  condition  wall  capture  the  general  behavior  of  the  flow 
control  devices  without  having  to  simulate  them  directly.  This  will  result  in  a  reduced  computational  cost 
by  eliminating  the  need  to  include  fine  grids  for  the  vortex  generators  and  synthetic  jets.  By  reducing  the 
computational  cost,  the  inlet  design  cycle  time  can  also  be  significantly  reduced.  High  resolution  calculations 
of  the  flow  about  micro- vortex  generator  fins  and  jets  on  a  flat  plate,  have  been  performed  using  a  Reynolds 
averaged  Navier-Stokes  code. 

Future  work  will  include  a  comparison  of  the  numerical  results  to  experimental  data  for  a  better  un¬ 
derstanding  of  the  flow  physics  of  these  devices.  Once  an  understanding  of  the  flow  physics  is  known,  the 
development  of  computational  boundary  conditions,  which  model  the  flow  control  devices,  can  be  performed. 

This  work  was  done  in  collaboration  with  Pieter  Buning  (NASA  Langley). 

Anisotropic  mesh  adaption 

Carlo  L.  Bottasso 

Many  problems  in  fluid  and  solid  mechanics  present  solution  features  that  are  inherently  directional. 
For  example  shock  waves  are  characterized  by  gradients  that  are  very  high  in  one  local  direction,  while 
substantially  smaller  in  the  other  two.  Similarly,  a  vortex  shed  by  the  tip  of  a  wing  presents  gradients  in  the 
radial  direction  that  are  much  higher  than  in  the  direction  of  its  core.  In  all  these  cases,  anisotropic  meshes 
are  highly  desirable  for  their  evident  economic  advantages. 

We  are  developing  a  methodology  for  introducing  regions  of  high  anisotropy  in  existing  isotropic  unstruc¬ 
tured  grids  in  complex,  curved,  three-dimensional  domains.  The  new  procedures  can  be  used  for  refining 
solution  features  internal  to  the  computational  domain  (e.g.,  shock  waves)  or  in  the  proximity  of  its  bound¬ 
aries  (e.g.,  boundary  layers).  In  both  cases,  suitable  voids  are  created  in  the  existing  grid  in  the  regions  of 
localization  using  a  mesh  motion  algorithm  that  solves  a  fictitious  elasticity  problem.  The  voids  are  then 
filled  with  stacks  of  prisms  that  are  subsequently  tetrahedronized  to  yield  a  simplicial  mesh.  The  mesh 
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motion  algorithm  allows  dealing  in  a  simple  and  effective  manner  with  the  problem  of  self-intersection  of 
elements  in  concave  regions  of  the  model  boundaries  and  in  the  case  of  closely  spaced  model  faces,  avoiding 
the  need  for  cross-over  checks  and  complex  grid  correction  procedures. 

The  capabilities  and  performance  of  the  new  methodology  are  being  tested  with  the  help  of  practical 
examples  involving  the  growth  of  boundary  layers  and  the  refinement  of  shock  waves. 

Discontinuous  dual-primal  mixed  finite  elements 

Carlo  L.  Bottasso 

Discontinuous  Galerkin  methods  are  the  subject  of  active  research  for  the  solution  of  partial  differential 
equations  in  a  number  of  problems  in  mathematical  physics.  Great  insight  on  the  characteristics  of  such 
methods  can  be  gained  from  the  study  of  model  problems,  for  example  the  Laplacian. 

We  have  developed  a  novel  discontinuous  mixed  finite-element  formulation  for  the  solution  of  second- 
order  elliptic  problems.  Fully  discontinuous  piecewise  polynomial  finite-element  spaces  are  used  for  the 
trial  and  test  functions.  The  discontinuous  nature  of  the  test  functions  at  the  element  interfaces  allows 
the  introduction  of  new  boundary  unknowns  that,  on  the  one  hand  enforce  the  weak  continuity  of  the  trial 
functions,  and  on  the  other  avoid  the  need  to  define  a  priori  algorithmic  fluxes  as  in  standard  discontinuous 
Galerkin  methods.  Static  condensation  is  performed  at  the  element  level,  leading  to  a  solution  procedure 
based  on  the  sole  interface  unknowns.  In  the  one-dimensional  case,  we  can  show  the  equivalence  of  the 
method  with  implicit  Runge-Kutta  schemes  of  the  collocation  type  exhibiting  optimal  behavior.  This  clearly 
provides  a  strong  incentive  towards  the  generalization  to  multiple  space  dimensions.  Numerical  experiments 
in  one  and  two  dimensions  demonstrate  the  order  accuracy  of  the  new  method,  confirming  the  results  of  the 
analysis. 

Future  efforts  will  concentrate  on  the  extension  of  these  ideas  to  other  model  problems,  and  on  compar¬ 
isons  of  this  class  of  methods  with  alternative,  well  established  solution  procedures. 

This  work  is  done  in  collaboration  with  S.  Micheletti  and  R.  Sacco  (Politecnico  di  Milano). 

Geometric  integration  theory  for  multibody  dynamics 

Carlo  L.  Bottasso 

The  partial  differential  equations  governing  the  dynamics  of  nonlinear  multibody  systems  composed  of 
beams,  shells  and  rigid  bodies  are  known  to  present  a  rich  mathematical  structure.  In  particular,  the  resulting 
models  are  Hamiltonian  systems  characterized  by  a  symplectic  nature  and  associated  with  conservation  laws 
that  stem  from  symmetries  of  the  Hamiltonian.  The  linear  and  angular  momentum  as  well  as  the  total 
mechanical  energy  are  conserved  for  free  motions  of  such  systems.  Classical  integration  methods  rarely 
preserve  the  underlying  structure  of  the  problem  being  solved,  and  hence,  such  structure  is  lost  in  the 
numerical  solution. 

We  have  developed  a  novel  integration  scheme  for  the  nonlinear  dynamics  of  geometrically  exact  shells 
based  on  the  inextensible  director  assumption.  This  method  extends  previous  work  done  on  the  dynamics 
of  rigid  bodies,  geometrically  exact  beams  and  flexible  multibody  systems.  The  new  algorithm  is  designed 
so  as  to  imply  the  strict  decay  of  the  system  total  mechanical  energy  at  each  time  step,  and  consequently 
unconditional  stability  is  achieved  in  the  fully  nonlinear  regime.  Furthermore,  the  scheme  features  tunable 
high  frequency  numerical  damping  and  it  is,  therefore,  stiffly  accurate.  The  method  was  implemented  in  a 
general-purpose  multibody  code.  Numerical  test  examples  have  been  solved  for  verifying  the  characteristics 
of  the  new  method. 
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Future  efforts  will  concentrate  on  the  use  of  the  code  for  solving  more  complex  problems,  in  particular 
for  the  modeling  of  helicopters  and  tiltrotors. 

This  work  is  done  in  collaboration  with  O.A.  Bauchau  (Georgia  Institute  of  Technology). 

Textbook  multigrid  efficiency  for  CFD  applications 

Boris  Diskin 

State-of-the-art  multigrid  methodologies  for  large-scale  compressible  flow  applications  use  a  block-matrix 
relaxation  and/or  a  pseudo-time-dependent  approach  to  solve  the  equations;  significant  improvements  have 
been  demonstrated  using  multigrid  approaches,  but  the  methods  are  not  optimally  convergent.  The  Reynolds- 
averaged  Navier-Stokes  (RANS)  equation  sets  are  systems  of  coupled  nonlinear  equations  which  are  not,  even 
for  subsonic  Mach  numbers,  fully  elliptic,  but  contain  hyperbolic  partitions.  There  is  a  potential  gain  of 
several  orders  of  magnitude  in  operation  count  reduction  if  the  textbook  multigrid  efficiency  (TME)  could  be 
attained  for  the  RANS  equations.  The  objective  of  the  research  is  to  develop  TME  methods  for  large-scale 
CFD  applications. 

The  method  proposed  for  achieving  TME  is  the  distributed  relaxation  approach  of  Brandt  which  de¬ 
composes  the  system  of  equations  into  separable,  usually  scalar,  factors  that  can  be  treated  with  optimal 
methods.  The  distributed  relaxation  is  applied  in  regular  (smoothly  varying)  flow  regions.  Near  boundaries 
and  discontinuities,  the  general  approach  is  to  relax  the  governing  equations  directly  in  terms  of  primitive 
variables.  Several  sweeps  of  robust  (but  possibly  slowly  converging)  relaxation,  such  as  Newton-Kacmarcz 
relaxation,  can  be  made  in  these  regions.  The  additional  sweeps  will  not  affect  the  overall  complexity  because 
the  number  of  boundary  and/or  discontinuity  points  is  usually  negligible  in  comparison  with  the  number  of 
interior  points.  Recently,  TME  has  been  demonstrated  for  viscous  compressible  flow  applications  and  for 
high-Reynolds-number  incompressible  viscous  flow.  The  main  accomplishment  of  this  half-year  period  is  the 
formulation  of  the  general  framework  expected  to  be  required  for  achieving  TME  for  large-scale  compressible 
flow  applications.  This  framework  was  tested  by  solving  the  quasi-one-dimensional  Euler  equations. 

Currently  efforts  are  directed  to  extension  of  the  proposed  TME  methodology  to  different  (subsonic, 
transonic,  and  supersonic;  compressible  and  incompressible;  free,  wake,  boundary-layer,  stagnation,  lifting, 
etc.)  flow  regimes  on  different  structured  (staggered  and  non-st aggered;  Cartesian  and  body-fitted)  grids  for 
different  (conservative  and  non-conservative;  low  and  high  order)  discretizations  of  Navier-Stokes  equations. 

The  research  in  these  directions  is  conducted  in  collaboration  with  with  J.L.  Thomas  (NASA  Langley) 
and  A.  Brandt  (The  Weizmann  Institute  of  Science). 

Materials  science:  Optimization  of  chemical  vapor  infiltration  with  simultaneous  powder 
formation 

Adi  Ditkowski 

A  variety  of  materials  are  produced  by  infiltration  processes.  In  these  techniques  a  fluid  phase  (i.e.,  a 
gas  or  a  liquid)  is  transported  into  a  porous  structure,  where  it  then  reacts  to  form  a  solid  product.  These 
methods  are  particularly  important  for  producing  composite  materials,  where  the  initial  porous  perform 
is  composed  of  the  reinforcement  phase  (i.e.,  fibers,  whiskers,  or  particles)  and  infiltration  produces  the 
matrix.  A  detailed  assessment  of  the  relevant  reaction  and  mass  transport  rates  during  infiltration  requires 
mathematical  modeling,  using  a  minimum  of  two  coupled  partial  differential  equations  which  describe  changes 
in  the  reactant  concentration  and  the  solid  structure  as  a  function  of  both  position  and  time.  A  key  difficulty 
in  chemical  vapor  infiltration  (CVI)  is  the  long  processing  times  that  are  typically  required.  With  this  in 
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mind,  it  is  important  to  minimize  infiltration  times.  The  goal  of  the  research  is  to  find  the  optimal  conditions 
(i.e.,  pressure  and  temperature)  to  run  this  process. 

This  study  specifically  considers  optimization  for  a  set  of  two  equations  which  describe  isothermal, 
isobaric  CVI.  This  work  was  divided  into  the  following  stages.  (1)  A  set  of  PDEs  which  model  the  process 
was  derived.  (2)  An  analytic  analysis  was  made  on  this  model  to  determine  its  proper  ties.  (3)  Asymptotic 
solutions  were  found  and  an  optimization  procedure  was  established  on  these  asymptotic  solutions.  These 
results  are  an  estimate  of  the  optimal  conditions.  (4)  Two  numerical  methods  were  used  to  verify  the 
asymptotic  results  and  to  study  their  limits.  (5)  The  numerical  optimization  procedure  was  developed  and 
used  for  this  problem.  Our  work  presents  a  framework  for  treating  the  problem  of  determining  the  optimal 
pressure  and  temperature  which  corresponds  to  the  minimum  infiltration  time. 

In  the  future  the  methodology  developed  here  will  be  used  to  find  the  optimal  operation  condition  for 
other  CVI  processes  as  well  as  other  problems  in  materials  science,  such  as  silicone  oxidation  and  trench- 
filling. 

This  work  was  done  in  collaboration  with  David  Gottlieb  and  Brian  W.  Sheldon. 

On  error  bounds  of  finite  difference  approximations  to  partial  differential  equations  temporal 
behavior  and  rate  of  convergence 

Adi  Ditkowski 

The  question  of  the  role  of  numerically  imposed  boundary  conditions  in  the  solution  of  parabolic  and 
hyperbolic  PDE’s  has  been  with  us  for  many  years.  Many  investigators  have  studied  the  effect  of  boundary 
conditions  on  the  stability  of  the  overall  scheme,  (i.e.,  the  inner  algorithm’  -I-  boundary  conditions).  In 
this  context  stability  implies  convergence  of  the  scheme,  at  a  fixed  time  £,  as  the  mesh  is  refined.  The 
question  of  the  temporal  behavior  of  the  error  was  usually,  not  considered.  When  constructing  higher  order 
schemes,  third-order  accuracy  and  above,  it  turns  out  that  it  is  difficult  to  state  boundary  conditions  such 
that  the  overall  scheme  remains  stable.  The  question  then  arises:  What  happens  to  the  overall  accuracy 
of  the  numerical  solution  if  the  order  of  accuracy  of  the  inner  stencil,  (m),  is  higher  than  the  order  of  the 
boundary  conditions,  (m  —  s).  This  problem  has  been  tackled  by  Gustafsson  in  his  papers  from  1975  and 
1981.  His  main  result,  for  both  parabolic  and  hyperbolic  PDE’s,  is  that  if  the  accuracy  of  the  extra  boundary 
conditions,  required  for  ‘numerical  closure’  of  the  problem,  are  one  less  than  that  of  the  inner  scheme,  then 
the  overall  accuracy  is  not  affected.  The  physical  boundary  conditions,  however,  must  be  approximated  to 
the  same  order  as  the  inner  scheme.  The  goal  of  this  work  was  to  derive  error  bounds  which  are  functions 
of  grid  space  and  time . 

In  this  work  we  had  considered  a  form  of  differentiation  matrices,  both  hyperbolic  and  parabolic,  which 
represent  a  fairly  wide  family  of  boundary  condition  formulations  plus  central  inner  schemes.  We  investigate 
the  dependence  of  the  error  on  time  as  well  as  on  mesh  size.  The  main  results  are  as  follows: 

•  In  the  parabolic  case,  the  overall  convergence  rate  is  of  order  m  if  s  =  0, 1,  and  m  -  s  4-  3/2  if  s  >  2. 
The  error  is  uniformly  bounded  independent  of  t  for  allt,  t  >  0. 

•  In  the  hyperbolic  case,  the  overall  convergence  rate  is  of  order  of  min  (m,  m  —  s  +  1)  for  all  s,  in 
agreement  with  the  results  given  by  Gustafsson.  For  s  =  0,1,  the  temporal  bound  of  the  error 
behaves  as  y/t  for  t  <  1,  and  bends  over  smoothly  to  a  linear  bound  as  t  increases.  For  s  >  2,  the 
temporal  behavior  is  ~  yft  for  all  f,  t  >  0. 

The  parabolic  and  hyperbolic  results  for  the  case  s  >  2  were  found  by  using  standard  energy  methods.  In 
order  to  derive  the  hyperbolic  results  for  the  case  s  =  1,  however,  a  new  method  had  to  be  developed.  The 
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approach  we  took  here  is  to  define  a  new  problem,  which  is  somewhat  similar  to  the  original  one  and  to  show 
that  the  solution  of  this  auxiliary  problem  bounds  the  error  of  the  original  problem. 

In  the  near  future  we  would  like  to  refine  the  bounds,  derive  error  bound  for  different  numerical  schemes 
and  to  use  the  new  approach  to  other  applications  such  as  stability  analysis  of  flows. 

This  work  was  done  in  collaboration  with  Saul  Abarbanel  and  Bertil  Gustafsson. 

High  order  finite  element  methods  for  electromagnetics 

D.  Gottlieb,  J.S.  Hesthaven,  and  T.  Warburton 

It  is  becoming  increasingly  apparent  that  there  is  a  need  in  the  electromagnetics  community  for  improved 
numerical  methods  for  simulation  of  electromagnetic  wave  propagation  in  complex  media.  With  this  in  mind, 
we  are  proposing  the  use  of  high  order  methods  that  can  be  used  on  unstructured  grids.  This  enables  the  use 
of  standard  finite  element  meshes  for  complex  geometries  with  the  additional  feature  that  field  profiles  can 
be  accurately  and  efficiently  propagated  through  the  domain  without  the  phase  errors  typically  seen  with 
lower  order  methods. 

We  have  created  a  two-dimensional  code  that  uses  triangular  elements  which  are  up  to  13th  order  in 
space  and  fourth  order  in  time  accurate.  We  have  tested  this  against  standard  test  cases  and  have  shown 
exponential  accuracy.  For  a  specific  case  we  see  that  increasing  the  order  of  the  element  decreases  the  error 
in  the  calculated  radar  cross  section  (RCS)  by  a  factor  of  ten.  Additionally,  we  have  shown  theoretically 
that  this  kind  of  convergence  will  be  expected  whenever  the  solution  we  are  trying  to  represent  is  smooth 
enough.  We  are  currently  collaborating  with  Sherwin  (Imperial  College)  and  Morgan  (University  of  Swansea) 
to  compare  the  efficiency  of  this  code  against  an  optimized  finite  volume  code. 

We  have  also  created  a  three-dimensional  code  that  uses  tetrahedral  elements  which  are  up  to  tenth 
order  in  space  and  fourth  order  accurate  in  time.  We  have  spent  the  summer  interfacing  this  code  with 
an  industrial  mesh  generator  (Gambit  from  Fluent).  We  have  also  verified  this  code  against  results  from 
Hypercomp  Inc.  and  have  shown  reasonable  agreement  with  bistatic  RCS  data  for  scattering  from  a  zero¬ 
thickness  business  card  geometry.  Additionally,  we  have  been  computing  results  for  geometries  (from  the 
EMCC  test  case  database)  suggested  .by  the  Electromagnetics  Research  Branch  at  NASA  Langley,  under 
the  supervision  of  Truong  Nguyen  and  Fred  Beck  (now  retired).  We  have  compared  results  from  the  3D 
code  with  results  from  an  axisymmetric  multi-domain  code  of  Hesthaven  for  scattering  from  a  geometry 
constructed  as  the  union  of  a  cone  and  a  hemisphere.  They  have  shown  good  agreement  for  incident  fields  at 
certain  angles,  however  this  test  case  has  indicated  that  further  work  is  required  to  understand  and  account 
for  the  singularities  in  fields  caused  by  the  sharp  point  of  the  cone. 

We  intend  to  continue  optimizing  the  three-dimensional  code  as  well  as  calculating  results  for  the  cone- 
sphere  scattering  problem  and  for  the  cavity  penetration  problems  suggested  by  Truong  Nguyen,  Fred  Beck 
and  M.C.  Bailey  from  the  Electromagnetics  Research  Branch.  These  geometries  have  proved  to  be  extremely 
interesting  in  the  way  that  they  stretch  certain  features  of  these  time-domain  methods.  These  difficulties  are 
certainly  not  unique  to  our  method  but  pose  problems  for  almost  all  numerical  methods  used  to  solve  the 
time-domain  Maxwell’s  equations.  We  also  intend  to  investigate  mixed  implicit /explicit  time  integration  to 
treat  the  thinly  coated  cone-sphere  problem  that  was  proposed  by  the  Electromagnetics  Research  Branch. 

Robust  optimization  including  model  uncertainties  and  reliability  constraints 

Luc  Huyse 

Current  formal  aerodynamic  optimization  procedures  assume  that  airfoil  performance  is  perfectly  pre¬ 
dicted  by  CFD  analysis  and  that  all  operating  conditions  (e.g.,  cruise  Mach  number)  are  fixed.  During  the 
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early  stages  of  the  design  process  many  parameters  are  merely  estimates  of  desired  operating  conditions. 
Single  point  designs  using  these  uncertain  parameters  can  lead  to  overly  optimistic  projections  of  as-built 
performance.  Such  assumptions  can  lead  to  designs  with  poor  performance  at  off-design  conditions.  This 
research  will  develop  new  optimization  methods  that  anticipate  parameter  uncertainty. 

A  theoretical  framework  for  robust  (or  stochastic)  optimization  has  been  developed  and  its  superiority 
over  existing  multi-point  optimization  techniques  has  been  demonstrated  through  analytic  examples.  Based 
on  a  second-order  expansion  of  a  mean- value  analysis,  substantial  computational  savings  can  also  be  achieved 
when  only  moderate  levels  of  uncertainty  exist  for  some  of  the  variables  or  parameters  in  the  mathematical 
model. 

The  newly  developed  techniques  were  applied  to  an  airfoil  optimization  problem,  using  the  FUN2D  code. 
The  minimization  of  the  drag  in  cruise  regime  while  maintaining  a  constant  lift  is  chosen  as  the  objective 
function.  Only  the  Mach  number  is  assumed  to  be  uncertain  in  this  preliminary  analysis.  The  performance 
of  the  resulting  robust  design  is  considerably  better  than  the  expected  value  of  the  performance  of  an  optimal 
design,  obtained  at  selected  deterministic  or  fixed  operating  conditions.  Currently,  we  are  combining  the 
FUN2D  code  with  the  iSight  optimization  framework  so  that  constrained  optimization  with  a  more  realistic 
set  of  design  variables  can  be  studied.  The  ultimate  goal  is  to  develop  general-purpose  robust  optimization 
methods  for  use  in  multidisciplinary  design. 

This  work  was  conducted  in  collaboration  with  Sharon  Padula  (NASA  Langley). 

Active  shielding  and  control  of  environmental  noise 

Josip  Loncaric ,  Victor  S.  Ryaben’kii ,  and  Semyon  V.  Tsynkov 

Rejection  of  exterior  noise  caused  by  periodic  sources  such  as  propellers  or  turbines  would  significantly 
enhance  passenger  comfort  and  reduce  noise  fatigue  on  long  flights.  Passive  sound  absorbing  materials  help  at 
high  frequencies,  but  to  be  effective  below  about  1  kHz  their  weight  penalty  becomes  significant.  Active  noise 
control  can  reduce  low  frequency  noise  with  less  weight  penalty.  We  present  the  mathematical  foundations 
of  a  new  active  technique  for  control  of  the  time-harmonic  acoustic  disturbances. 

Unlike  many  existing  methodologies,  the  new  approach  provides  for  the  exact  volumetric  cancellation 
of  the  unwanted  noise  in  a  given  predetermined  region  of  space  while  leaving  those  components  of  the  total 
sound  field  deemed  as  friendly  unaltered  in  the  same  region.  Besides,  the  analysis  allows  us  to  conclude 
that  to  eliminate  the  unwanted  component  of  the  acoustic  field  in  a  given  area,  one  needs  to  know  relatively 
little;  in  particular,  neither  the  locations  nor  the  structure  of  the  external  noise  sources  need  to  be  known. 
We  constructed  the  general  solution  for  the  aforementioned  noise  control  problem.  The  apparatus  used  for 
deriving  the  general  solution  is  closely  connected  to  the  concepts  of  generalized  potentials  and  boundary 
projections  of  Calderon’s  type.  To  prove  that  the  new  technique  is  appropriate,  we  thoroughly  worked  out 
a  simple  two-dimensional  model  example  that  allows  full  analytical  consideration,  including  optimalization 
of  the  control  effort. 

In  order  to  develop  numerically  computable  solutions,  we  plan  to  describe  the  discrete  framework  for  the 
noise  control  problem  parallel  to  the  continuous  one.  This  discrete  framework  is  obtained  using  difference 
potentials  method;  in  the  future  it  is  going  to  be  used  for  analyzing  complex  configurations  that  originate 
from  practical  designs.  Once  we  have  computed  the  solution  for  a  particular  configuration,  we  intend  to 
investigate  the  possibilities  of  optimizing  it  according  to  the  different  criteria  that  would  fit  different  practical 
requirements.  We  expect  to  discuss  the  applicability  of  the  technique  to  quasi-stationary  problems,  future 
extensions  to  the  cases  of  the  broad-band  spectra  of  disturbances,  as  well  as  other  possible  applications 
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which  may  include  different  physics,  such  as  electrodynamics,  and  different  formulations  of  the  boundary- 
value  problems,  such  as  scattering. 

Performance  of  different  communication  patterns  on  the  Coral  system 

Ignacio  M.  Llorente 

The  evaluation  of  the  message  passing  performance  of  a  parallel  system  should  not  only  consider  the 
performance  of  the  interconnection  system,  but  also  the  impact  of  the  memory  hierarchy  exploitation.  The 
communication  cost  in  a  parallel  application  depends  not  only  on  the  amount  of  data  exchanged  between  the 
processes,  but  also  on  how  they  are  structured  in  the  local  memories.  In  fact,  the  bandwidth  reduction  due 
to  local  memory  access  with  poor  data  locality  may  be  more  important  than  the  reduction  due  to  contention 
in  the  network.  This  situation  might  influence,  for  example,  the  optimal  data  partitioning  of  domain  decom¬ 
position  applications.  Although  3-D  decompositions  exhibit  lower  inherent  communication-to-computation 
ratios,  lower  dimensional  decompositions,  where  boundaries  with  poor  spatial  locality  are  not  needed,  may 
be  more  efficient.  This  behavior  has  been  previously  reported  for  some  current  parallel  architectures  as  the 
Cray  T3E  and  the  SGI  Origin  2000.  The  aim  of  this  research  is  to  evaluate  the  performance  of  the  Coral 
system  for  different  communication  patterns  using  LAM/MPI  6.3.  Since  the  study  is  focused  on  the  Coral 
subsystem  based  on  dual-CPU  nodes  (500MHz  Pentium  III),  we  have  also  studied  the  influence  of  the  pro¬ 
cessor  mapping.  In  this  case,  an  optimal  mapping  does  not  only  reduce  network  contention,  but  also  allows 
the  exploitation  of  shared  memory  communications. 

We  have  performed  the  following  communication  tests: 

•  Unidirectional  Communications:  First,  we  have  used  the  classical  Ping-Pong  test  to  establish  the 
achievable  bandwidths.  Then,  this  base  experiment  has  been  modified  by  increasing  the  network  load 
and  by  decreasing  the  spatial  locality  properties  of  the  messages.  For  contiguous  messages,  prelim¬ 
inary  results  show  an  asymptotic  bandwidth  of  80  Mbytes/s  between  processors  in  the  same  node 
(shared  memory  communication)  and  11  Mbytes/s  for  processors  in  different  nodes.  For  strided 
data,  the  bandwidth  drops  down  to  20  (25%  the  asymptotic  bandwidth)  and  8  Mbytes/s  (73%) 
respectively.  With  network  saturation,  the  bandwidth  between  different  nodes  for  contiguous  and 
strided  data  is  6  (55%)  and  5  Mbytes/s  (45%)  respectively.  Summing  up  these  results,  for  inter-node 
communications,  the  bandwidth  reduction  due  to  network  saturation  is  more  important  than  the  re¬ 
duction  due  to  non-contiguous  messages,  while  for  intra-node  communications,  message  performance 
strongly  depends  on  data  locality  and  is  not  affected  by  the  network  contention. 

•  Processor  Mapping:  It  is  clear  that  an  optimal  processor  mapping  should  exploit  the  faster  commu¬ 
nication  of  contiguous  messages  inside  the  nodes  which  also  helps  to  reduce  the  network  contention. 
For  example,  the  processor  mapping  of  a  1-D  topology  should  be  made  to  assign  neighboring  pro¬ 
cesses  to  the  same  node.  The  communication  bandwidth  has  been  measured  to  be  30%  higher  in 
such  a  situation. 

•  Domain  Decomposition  Technique:  The  relation  between  computation  and  communication  times  is 
much  higher  on  the  Coral  system  than  in  current  parallel  architectures,  so  the  3-D  decomposition 
achieves  the  best  performance,  requiring  30%  less  time  than  the  1-D  decomposition  and  8%  less  than 
the  2-D  one. 

•  Set  of  Tridiagonal  Systems  of  Equations:  A  set  of  tridiagonal  systems  of  equations  has  been 
solved  using  different  approaches:  pipelined  sequential  algorithms  (Pipelined  Gaussian  Elimination 
method),  matrix  transpose  algorithms  (Mapping  Transposition  algorithm)  and  different  parallel  im- 
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plicit  solvers  (Cyclic  Reduction  method,  Mattor’s  algorithm,  and  Wang’s  partition  method).  Each 
one  of  those  methods  exhibits  a  different  communication  pattern:  pipelined  communication,  col¬ 
lective  communications,  binary  tree  communications,  etc.  The  Mattor’s  algorithm  combined  with 
Cyclic  Reduction  shows  the  best  results  in  terms  of  execution  time  and  so  parallel  efficiency. 

The  work  performed  will  provide  guidance  for  efficient  programming  on  the  Coral  system. 

The  results  will  be  compared  with  those  obtained  on  some  current  parallel  architectures  (Cray  T3E  and 
SGI  Origin  2000)  considering  not  only  performance  issues,  but  also  the  performance-cost  ratio. 

This  research  was  conducted  in  collaboration  with  M.  Prieto-Matias  and  E.  Huedo-Cuesta  (Universidad 
Complutense,  Spain). 

A  robust  and  parallel  multigrid  method  for  the  simulation  of  a  yawed  flat  plate 

Ignacio  M.  Llorente 

The  combination  of  implicit  plane  relaxation  with  partial  and  full  coarsening  has  been  found  to  be  fully 
robust  for  solving  the  incompressible  Navier-Stokes  equations  on  highly  stretched  grids  in  the  simulation  of 
the  driven  cavity  and  a  boundary  layer  over  a  flat  plate.  In  particular,  the  combination  of  plane  smoothing 
and  semi-coarsening  achieves  a  convergence  rate  independent  of  grid  size,  stretching  and  Reynolds  number 
for  the  boundary  layer  simulation.  Moreover,  a  tri-plane  ordering  in  the  smoothing  sweep  exhibits  similar 
convergence  to  the  lexicographic  ordering  and  allows  the  efficient  parallel  implementation  of  the  algorithm. 
The  objective  of  this  research  is  to  study  whether  those  convergence  properties  are  maintained  in  the  sim¬ 
ulation  of  a  boundary  layer  at  high  Reynolds  numbers  over  a  yawed  flat  plate,  where  the  leading  edge  is 
not  perpendicular  to  the  stream.  This  simulation  exhibits  some  basic  problems  that  may  prevent  optimal 
multigrid  efficiencies  from  being  achieved,  namely,  a  highly  stretched  grid  near  the  yawed  plate  in  order  to 
capture  small  scale  physical  phenomena  and  an  entering  flow  not  aligned  with  the  discretization  grid. 

A  FAS  multigrid  code  has  been  used  to  solve  the  incompressible  Navier-Stokes  equations  on  a  single¬ 
block  grid  with  stretching.  The  discrete  equations  are  obtained  using  a  finite  volume  approach  in  a  staggered 
grid.  For  the  convective  terms  a  first  order  upwind  discretization  is  used.  Second-order  accuracy  is  achieved 
with  a  defect- correction  procedure  based  on  a  QUICK  scheme  inside  the  multigrid  cycle.  The  smoothing 
operator  is  a  cell-implicit  Symmetric  Coupled  Gauss  Seidel  method,  where  the  momentum  and  continuity 
equations  are  relaxed  in  a  coupled  manner.  The  combination  of  plane  smoothing  (where  each  plane  consists 
of  a  slab  of  cells)  and  semi-coarsening  has  been  found  to  be  robust  in  the  simulation  of  the  yawed  flat  plate 
for  angles  up  to  45  degrees.  The  smoothing  process  is  not  computationally  expensive  because  each  slab  of 
cells  is  approximately  solved  with  a  2-D  robust  multigrid  algorithm  consisting  of  one  FAS  F(l,l)  cycle  that 
combines  a  cell-implicit  smoother  and  semi-coarsening.  Textbook  multigrid  convergence  is  attained  for  the 
model  problem,  i.e.,  preliminary  results  show  that  the  convergence  rate  is  independent  of  the  grid  size,  grid 
stretching,  Reynolds  number,  and  angle  of  yaw.  In  fact,  the  full  multigrid  algorithm  converges  the  solution 
to  below  the  truncation  error  with  one  F(2,l)  cycle  per  level.  The  parallel  version  of  the  code  is  being 
implemented  on  the  Coral  system  by  using  the  MPI  communication  library  in  order  to  obtain  results  in  finer 
grids  up  to  256  cells  per  side. 

We  intend  to  continue  working  on  parallel  and  robust  multigrid  methods  for  block-structured  applica¬ 
tions.  In  particular,  we  will  compare  the  numerical  and  architectural  properties  of  coupled  and  distributive 
relaxation  for  the  incompressible  Navier-Stokes  equations. 

This  research  was  conducted  in  collaboration  with  R.S.  Montero  and  M.  Prieto-Matias  (Universidad 
Complutense,  Spain),  and  M.D.  Salas  (ICASE). 
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Large  eddy  simulation  using  a  parallel  multigrid  solver 

Dimitri  J.  Mavriplis 

The  failure  to  develop  a  universally  valid  turbulence  model  coupled  with  recent  advances  in  computa¬ 
tional  technology  have  generated  a  greater  interest  in  the  large-eddy  simulation  approach  for  computing  flows 
with  large  amounts  of  separation.  This  approach  involves  resolving  the  large-scale  unsteady  turbulent  eddies 
down  to  a  universally  valid  range  in  the  hope  of  yielding  a  more  generally  valid  simulation  tool.  The  purpose 
of  this  work  is  to  develop  a  large-eddy  simulation  capability  based  on  an  existing  unstructured  grid  Navier- 
Stokes  solver.  The  use  of  unstructured  grids,  which  facilitates  the  discretization  of  complex  geometries  and 
adaptive  meshing  techniques,  is  expected  to  enhance  the  flexibility  of  the  resulting  simulation  capability. 

An  unsteady  Reynolds-averaged  Navier-Stokes  (RANS)  flow  solver  based  on  unstructured  meshes  has 
been  developed  and  validated  on  the  case  of  a  circular  cylinder.  A  Detached  Eddy  Simulation  (DES)  model 
based  on  modifications  to  the  one-equation  Spalart-Allmaras  RANS  turbulence  model  has  been  implemented 
and  validated  for  the  case  of  flow  over  a  sphere.  Flow  over  a  wing  at  post-stall  incidences  has  also  been 
computed  using  the  unsteady  RANS  and  the  DES  methods.  These  computations  are  being  performed  on 
the  ICASE  PC  cluster,  Coral. 

The  simulation  of  decaying  isotropic  turbulence  in  a  periodic  domain  is  underway  in  order  to  enable  the 
calibration  of  the  model  constant  and  to  study  the  effect  of  artificial  dissipation  on  the  simulation  of  eddies 
in  the  inertial  range.  Higher  resolution  simulations  for  the  stalled  wing  and  flow  over  a  bluff  landing  gear 
geometry  are  to  be  performed  in  the  near  future. 

This  work  is  being  carried  out  in  collaboration  with  Juan  Pelaez  (Old  Dominion  University). 

Unstructured  multigrid  algorithms  for  the  solution  of  radiation  diffusion  problems 

Dimitri  J.  Mavriplis 

Under  the  ASCI  program,  the  simulation  of  radiation  transport  phenomena  has  been  identified  as  one 
of  the  most  time  consuming  elements  within  large  simulation  codes.  The  objective  of  this  research  is  to 
investigate  the  effectiveness  of  unstructured  multigrid  algorithms  in  efficiently  solving  unsteady  radiation 
transport  problems,  in  the  radiation  diffusion  limit. 

A  two-dimensional  unstructured  grid  radiation  diffusion  solver  has  been  developed  which  employs  ag¬ 
glomeration  multigrid  to  converge  the  nonlinear  problem  at  each  physical  time  step.  In  a  first  approach, 
a  nonlinear  FAS  multigrid  scheme  is  used  to  directly  converge  the  nonlinear  problem  at  each  time  step, 
whereas  in  the  second  approach,  a  Newton  iteration  scheme  is  applied  to  the  nonlinear  problem  and  a  linear 
multigrid  algorithm  is  employed  to  solve  the  linear  system  at  each  Newton  iteration.  A  hybrid  scheme  based 
on  a  nonlinear  FAS  multigrid  scheme  but  using  a  linear  solver  on  each  grid  level  has  also  been  implemented. 
Results  indicate  that  while  all  three  schemes  converge  at  the  same  rate  based  on  the  number  of  (linear  or 
nonlinear)  grid  sweeps,  the  linear  multigrid  solver  is  more  efficient  due  to  the  smaller  number  of  expensive 
nonlinear  function  evaluations. 

These  results  are  to  be  extended  into  the  three-dimensional  setting  and  effects  of  scalability  on  large 
numbers  of  processors  will  be  studied  in  addition  to  algorithmic  efficiency.  A  complete  validation  of  the 
solver  in  terms  of  accuracy  will  also  be  completed  in  the  near  future. 
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On  the  extension  of  compressible  flow  solvers  to  incompressible  flows 

C.-C.  Rossow 

The  extension  of  methods  designed  to  solve  compressible  flows  near  the  incompressible  flow  regime  still 
remains  a  challenge.  The  usual  way  to  make  compressible  codes  applicable  for  the  solution  of  low  Mach 
number  flows  is  to  precondition  the  system  of  compressible  equations.  This  is  done  in  such  a  way  that  the 
disparity  in  the  eigenvalues  is  reduced  and  the  condition  number  of  the  resulting  system  becomes  of  the  order 
of  unity.  However,  in  regions  of  very  low  Mach  number,  e.g.,  near  stagnation  points  and  in  recirculation 
regions,  the  preconditioning  matrix  may  become  singular  and  the  robustness  of  the  computation  may  be 
impaired.  Therefore,  to  make  compressible  codes  applicable  to  incompressible  flow  problems  with  the  same 
reliability  as  for  compressible  flows,  further  effort  seems  to  be  required  to  understand  the  basic  mechanisms 
of  computing  incompressible  flows. 

Due  to  the  robustness  problems  encountered  when  using  preconditioning  techniques,  the  question  may  be 
asked  whether  some  elements  are  missing  in  this  approach.  It  is  well  known  that  in  the  incompressible  limit, 
the  velocity  field  of  the  flow  must  be  divergence  free.  This  is  an  essential  requirement,  and  it  may  be  argued 
whether  this  condition  has  to  be  satisfied  on  a  discrete  level,  too.  Analysis  of  unsuccessful  computations 
with  preconditioning  revealed  that  this  constraint  was  violated  in  these  cases.  Based  on  this  observation, 
a  possibility  was  sought  to  respect  this  constraint  for  the  discrete  equations.  Based  on  the  Mach  number 
expansion  of  the  classical  Flux  Difference  Splitting  scheme,  the  discrete  mass  flux  was  investigated.  In  the 
incompressible  limit,  the  discrete  mass  flux  is  mainly  governed  by  pressure  differences  arising  in  the  artificial 
dissipation  due  to  the  solution  of  a  Riemann  problem.  Using  these  pressure  differences,  an  equation  for  a 
correction  of  the  pressure,  such  that  a  divergence  free  velocity  field  will  be  enforced,  could  be  derived.  In 
the  incompressible  limit,  this  equation  is  of  purely  elliptic  nature.  Solving  this  equation  in  the  framework  of 
an  unstructured,  compressible  code  allowed  the  computation  of  completely  incompressible  (M=0)  flows. 

It  could  be  shown  that  this  formulation  of  a  pressure  correction  equation  is  equivalent  to  classical  methods 
for  solving  incompressible  flows.  However,  in  the  present  formulation,  no  heuristic  arguments  are  required, 
but  the  derivation  is  based  on  the  solution  of  the  corresponding  Riemann  problem  in  combination  with  the 
divergence  free  velocity  field  constraint.  As  a  result,  the  discrete  incompressible  method  arises  consistently 
from  the  compressible  formulation. 

During  the  research,  the  observation  was  made  that  the  pressure  equation  is  not  well  suited  for  the 
computation  of  compressible  flows  with  shocks.  This  may  be  consistently  understood  with  the  vanishing 
pressure  differences  in  the  mass  flux  when  the  Mach  number  approaches  unity.  An  attempt  was  made 
to  turn  the  basic  elliptic  equation  into  a  hyberbolic  equation  for  supersonic  flows,  as  is  practice  for  other 
pressure-based  methods  used  for  the  computation  of  compressible  flows.  However,  the  experience  was  made 
that  very  small  timesteps  were  required  to  obtain  a  stable  integration  towards  steady  state,  and  this  approach 
was  recognized  as  not  efficient  for  compressible  flows.  Therefore,  a  blending  technique  was  derived  which 
allows  a  transition  from  the  pressure-based  incompressible  formulation  to  the  density-based  compressible 
formulation. 

With  this  formulation,  completely  incompressible  and  compressible,  transonic  flows  around  airfoils  were 
computed.  The  method  was  used  to  simulate  inviscid  as  well  as  viscous  flowfields,  and  the  implementation 
into  the  multigrid  framework  of  the  basic  unstructured  code  was  straightforward.  A  comparison  with  the 
standard  preconditioning  technique  used  in  the  same  code  showed  an  increased  robustness  of  the  pressure 
correction  method  for  high  lift  flows  close  to  separation.  This  may  be  attributed  to  the  fact  that  the  pressure 
correction  equation  enforces  the  divergence  free  velocity  field  on  a  discrete  level  during  each  timestep.  It  was 


found  that  for  highly  stretched  meshes,  the  elliptic  equation  became  stiff  due  to  the  high  aspect  ratios,  and 
the  number  of  Jacobi-iterations  had  to  be  increased. 

The  results  obtained  during  this  work  indicate  that  it  is  crucial  for  convergence  to  respect  the  divergence 
free  constraint.  Furthermore,  care  has  to  be  taken  when  regarding  incompressible  or  nearly  incompressible 
unsteady  flow:  this  constraint  applies  here  also,  and  the  violation  of  it  may  lead  to  wrong  interpretations. 
This  holds  especially  for  the  computation  of  flows  at  maximum  lift  conditions,  where  unsteadiness  is  likely 
to  occur,  and  recirculation  zones  with  very  low  speeds  are  present. 

Future  work  will  be  directed  towards  an  enhanced  method  for  the  solution  of  the  pressure  correction 
equation  to  ensure  convergence  for  more  severe  conditions.  Furthermore,  it  will  be  investigated  whether  the 
additional  evaluation  of  the  continuity  equation  could  be  avoided  to  reduce  computational  cost. 

High-order  discontinuous  Galerkin  method  and  WENO  schemes 

Chi- Wang  Shu 

Our  motivation  is  to  have  high  order  non-oscillatory  methods  for  structured  and  unstructured  meshes 
which  are  easy  to  implement  for  parallel  machines.  The  objective  is  to  develop  and  apply  high  order  dis¬ 
continuous  Galerkin  finite  element  methods  and  weighted  ENO  schemes  for  convection  dominated  problems. 
The  applications  will  be  problems  in  aeroacoustics  and  other  time  dependent  problems  with  complicated 
solution  structure. 

Jointly  with  Harold  Atkins  at  NASA  Langley,  we  are  continuing  in  the  investigation  of  developing  the 
discontinuous  Galerkin  method  to  solve  the  convection  dominated  convection-diffusion  equations.  We  have 
been  continuing  in  our  investigation  of  local  preconditioners  which  hopefully  would  dramatically  increase  the 
time  step  restriction  on  the  time  dependent  simulations  and  make  convergence  to  steady  states  much  faster, 
both  for  convection  problems  and  for  convection-diffusion  problems.  Jointly  with  Cockburn,  Luskin  and  Suli, 
we  have  been  developing  a  local  post  processing  technique  for  the  discontinuous  Galerkin  method  applied 
to  linear  PDEs  which  will  effectively  double  the  order  of  accuracy  of  the  scheme  with  minimal  additional 
computational  cost. 

Research  will  be  continued  for  high  order  discontinuous  Galerkin  methods  and  weighted  ENO  methods 
and  their  applications. 

High-order  two-way  artificial  boundary  conditions  for  nonlinear  wave  propagation  with 
backscattering 

Semyon  Tsynkov 

When  solving  linear  scattering  problems,  one  typically  first  solves  for  the  impinging  wave  in  the  absence 
of  obstacles.  Then,  by  linear  superposition,  the  original  problem  is  reduced  to  one  that  involves  only  the 
scattered  waves  driven  by  the  values  of  the  impinging  field  at  the  surface  of  the  obstacles.  In  addition,  when 
the  original  domain  is  unbounded,  special  artificial  boundary  conditions  (ABCs)  that  would  guarantee  the 
reflectionless  propagation  of  waves  have  to  be  set  at  the  outer  boundary  of  the  finite  computational  domain. 
The  situation  becomes  conceptually  different  when  the  propagation  equation  is  nonlinear.  In  this  case  the 
impinging  and  scattered  waves  can  no  longer  be  separated,  and  the  problem  has  to  be  solved  in  its  entirety. 
In  particular,  the  boundary  on  which  the  incoming  field  values  are  prescribed,  should  transmit  the  given 
incoming  waves  in  one  direction  and  simultaneously  be  transparent  to  all  the  outgoing  waves  that  travel  in 
the  opposite  direction.  We  call  this  type  of  boundary  conditions  two-way  ABCs. 

In  the  paper,  we  construct  the  two-way  ABCs  for  the  nonlinear  Helmholtz  equation  that  models  the  laser 
beam  propagation  in  a  medium  with  nonlinear  index  of  refraction.  In  this  case,  the  forward  propagation  is 
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accompanied  by  backscattering,  i.e.,  generation  of  waves  in  the  direction  opposite  to  that  of  the  incoming 
signal.  Our  two-way  ABCs  generate  no  reflection  of  the  backscattered  waves  and  at  the  same  time  impose 
the  correct  values  of  the  incoming  wave.  The  ABCs  are  obtained  for  a  fourth-order  accurate  discretization  to 
the  Helmholtz  operator;  the  fourth-order  grid  convergence  is  corroborated  experimentally  by  solving  linear 
model  problems.  We  also  present  solutions  in  the  nonlinear  case  using  the  two-way  ABC  which,  unlike  the 
traditional  Dirichlet  boundary  condition,  allows  for  direct  calculation  of  the  magnitude  of  backscattering. 

Future  work  will  include  building  a  more  efficient  iteration  procedure  than  the  one  currently  employed, 
numerical  studies  of  both  subcritical  and  supercritical  regimes,  comparing  the  results  obtained  in  the  frame¬ 
work  of  the  nonlinear  Helmholtz  equation  with  those  obtained  for  the  nonlinear  Schroedinger  equation, 
and  finally  casting  the  ABCs  that  we  have  constructed  into  the  general  framework  of  Calderon’s  pseudo- 
differential  boundary  equations. 

This  research  was  conducted  in  collaboration  with  Gadi  Fibich  (Tel  Aviv  University). 
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COMPUTER  SCIENCE 


Preconditioning  techniques  for  solving  the  particulate  flow  problem  in  parallel 

Abdelkader  Baggag 

To  extract  the  implicit  information  that  is  in  the  equations  of  the  motion  of  the  flow  of  particulates,  it  is 
necessary  to  numerically  solve  the  coupled  system  of  PDEs  consisting  of  the  equations  of  the  fluid  motion,  and 
the  equations  of  the  rigid-body  motion  governing  the  particles.  These  equations  are  coupled  through  the  no¬ 
slip  boundary  condition  on  the  particle  surfaces,  and  through  the  hydrodynamic  forces  and  torques  exerted 
by  the  fluid  on  the  particles.  To  solve  numerically,  the  particulate  flow  problem  is  discretized  by  a  standard 
Galerkin  finite  element  scheme,  in  which  both  the  fluid  and  particle  equations  of  motion  are  incorporated 
into  a  single  variational  equation  where  the  fluid  and  particle  velocities  appear  as  primitive  unknowns.  The 
computation  is  then  performed  on  an  unstructured  grid,  and  an  arbitrary  Lagrangian-Eulerian  moving  mesh 
technique  is  adopted  to  deal  with  the  motion  of  the  particles. 

The  numerical  simulation  of  the  particulate  flow  problem  is  extremely  computationally  intensive,  since 
it  is  nonstationary,  i.e.,  it  evolves  over  many  time  steps,  and  within  each  time  step  a  nonlinear  system  of 
equations  must  be  solved,  and  within  each  nonlinear  iteration  several  linear  systems  are  solved.  The  time 
stepping  is  handled  by  a  backward  Euler  scheme,  and  a  Newton  method  is  used  to  deal  with  the  nonlinearity. 
The  various  algebraic  systems  (Jacobians)  which  arise  thereafter  are  large,  sparse,  nonsymmetric,  and  often 
indefinite,  and  their  solution  can  consume  up  to  90%  of  the  total  CPU  time  of  the  entire  simulation.  The 
main  focus  of  this  research  is  to  address  the  efficient  solution  of  the  linear  systems  and  to  implement  them 
on  parallel  computers,  and  to  optimize  them  in  particular  on  the  Origin  2000. 

Although  the  methods  based  on  the  Krylov  subspaces  have  excellent  parallelization  properties,  they 
lack  robustness  without  good  preconditioning.  In  this  research  some  novel  techniques  are  investigated  and 
devised.  These  techniques  include  the  balance  scheme,  the  multilevel  method,  and  the  augmentation  method. 
These  algorithms  are  implemented  in  a  Particle  Mover  object-oriented  code,  called  PM+- {-,  and  extensively 
validated  for  the  particulate  flow  problem  of  the  Grand- Challenge  project. 

Parallel  implicit  solvers  for  PDEs  and  PDE-constrained  optimization 

David  E.  Keyes 

The  development  and  application  of  parallel  implicit  solvers  for  multiscale  phenomena  governed  by  PDEs 
are  our  chief  objectives.  Newton- Krylov-Schwarz  (NKS)  methods  have  proven  to  be  broadly  applicable, 
architecturally  versatile,  and  tunable  for  high  performance  on  today’s  high-end  commercial  parallel  platforms 
(e.g.,  Cray  T3E,  SGI  Origin,  IBM  SP,  Compaq  AlphaServer).  The  challenge  of  obtaining  high  performance 
is  likely  to  be  increased  on  coming  generations  of  high-end  machines,  as  represented,  for  instance,  by  the 
ASCI  “white”  machine  at  Lawrence  Livermore  National  Laboratory,  and  also  on  Beowulf  clusters,  such  as 
ICASE’s  “Coral.”  Our  primary  efforts  are  concentrated  on  algorithmic  adaptations  of  NKS  methodology 
appropriate  for  emerging  architectures  and  on  evaluation  of  new  software  tools  and  methodology. 

Since  we  employ  Newton’s  method  as  the  solver,  we  can  also  apply  our  techniques  to  nonlinear  systems 
that  augment  the  original  PDE  discrete  equations  by  additional  equations  that  represent  the  first-order 
necessary  conditions  (the  Karush-Kuhn- Tucker  equations)  for  optimality.  Working  through  the  Lagrange 
formulation  for  PDE-constrained  optimization,  we  designate  this  approach  Lagrange-NKS,  or  simply  LNKS. 

Both  structured-grid  and  unstructured-grid  CFD  legacy  codes  have  been  ported  to  advanced  platforms 
and  reasonable  objectives  for  algorithmic  convergence  rate,  parallel  efficiency,  and  raw  floating  point  perfor- 
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mance  have  been  met.  The  general  approach  embodied  in  the  NKS  family  of  algorithms  is  documented  in 
previous  ICASE  technical  reports,  among  other  places.  Newton’s  method,  robustified  by  pseudo-transient 
continuation,  generates  global  linear  systems  that  are  solved  in  a  matrix-free  manner  by  Krylov  iteration, 
and  preconditioned  locally.  In  the  absence  of  timestepping,  a  non-local  two-level  preconditioner  based  on  a 
coarsening  of  the  problem  is  also  required  to  maintain  algorithmic  scalability  at  a  modest  cost  to  asymptotic 
parallel  efficiency,  but  most  nonlinear  problems  require  some  form  of  timestepping,  even  for  the  computations 
of  steady  states. 

During  the  past  year,  our  1999  Bell  Prize-winning  work  —  an  implementation  of  the  NASA  code  FUN3D 
running  on  up  to  3,072  dual-processor  nodes  of  the  ASCI  Red  machine  at  227  Gigaflops/s  —  has  been 
extended  by  algorithmic  performance  tuning  and  memory  system  modeling,  while  we  await  access  to  the 
next  larger  machine. 

Domain-decomposition  algorithms  for  PDEs  and  Lagrange-based  algorithms  for  PDE-constrained  opti¬ 
mization  are  special  cases  of  partitioned  solvers,  where  the  partitions  are  geometrically  or  functionally  based, 
respectively.  We  are  examining  a  variety  of  partitioned  solvers,  where  the  variety  arises  from  the  type  of 
partitioning  and  also  the  order  of  nestedness  of  the  partitioning  with  respect  to  the  linearization  required 
in  any  nonlinear  problem.  A  key  decision  in  the  trade-off  between  globalized  convergence  properties  and 
parallel  computational  cost  per  iteration  is  whether  a  global  Newton  linearization  is  employed,  or  multiple 
independent  Newton  linearizations  within  partitions.  We  have  recently  begun  to  explore  nonlinear  precon¬ 
ditioning  which  is  expected  to  be  well  suited  for  machines  with  10,000  processors  or  more.  In  this  range 
conventional  NKS  methods  begin  to  lose  efficiency  to  the  cost  of  global  synchronization. 

We  will  continue  to  develop  NKS  methods  in  implicit  parallel  CFD  and  radiation  transport,  examining 
a  variety  of  algorithmic,  programming  paradigm,  and  architectural  issues. 

The  Coral  Project 

Josip  Loncaric,  Thomas  W.  Crockett ,  Piyush  Mehrotra,  Peter  Kearney ,  and  Manuel  D.  Salas 

The  cost  of  developing  complex  computer  components  such  as  CPUs  has  become  so  high  that  scientific 
applications  alone  cannot  carry  the  full  burden.  In  the  future,  scientific  computing  will  have  to  use  mass 
market  leverage  to  overcome  the  cost  barrier.  A  cost-effective  alternative  to  high-end  supercomputing  was 
pioneered  by  Beowulf,  a  cluster  of  commodity  PCs.  By  now,  high  performance  Beowulf  clusters  can  be  built 
using  fast  commodity  PCs  and  switched  Fast  Ethernet.  We  are  exploring  the  benefits  and  the  limitations  of 
this  approach,  based  on  applications  of  interest  to  ICASE. 

The  initial  phase  of  the  Coral  project,  consisting  of  32  Pentium  II  400  MHz  nodes  and  a  dual-CPU 
server,  demonstrated  aggregate  peak  performance  in  excess  of  10  Gflops/s,  with  sustained  performance  on 
CFD  applications  of  about  1.5  Gflops/s.  In  order  to  provide  a  richer  environment  for  further  experimentation, 
a  dual-CPU  configuration  was  chosen  for  the  second  phase  of  the  Coral  project.  We  have  added  16  dual 
Pentium  III  500  MHz  machines  and  two  dual-CPU  file  servers.  The  third  stage  of  this  project  added  16  dual 
Pentium  III  800  MHz  machines  and  a  32-node  low  latency  1.25  Gbps  Giganet  cLAN  network  fabric.  The 
resulting  system  contains  96  compute  CPUs  with  an  aggregate  of  36  GB  of  RAM  and  981  GB  of  local  disk 
space.  Coral’s  three  dual-CPU  servers  provide  an  additional  1.5  GB  of  RAM  and  246  GB  of  disk  space. 

Coral  has  an  excellent  price/performance  ratio,  almost  an  order  of  magnitude  better  than  an  equivalent 
proprietary  supercomputer  design.  This  conclusion  is  based  on  our  experience  with  a  variety  of  applications, 
ranging  from  coarse-grained  domain  decomposition  codes  to  communication-intensive  parallel  Tenderers. 

For  a  computer  system  with  this  much  memory  and  disk  capacity,  the  ability  to  move  data  around 
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within  the  cluster,  as  well  as  to  and  from  other  systems,  becomes  an  important  consideration-  In  conjunc¬ 
tion  with  the  second  phase  of  the  Coral  project,  we  have  built  a  small  Gigabit  Ethernet  testbed  which 
connects  the  Coral  servers,  one  of  the  dual-CPU  compute  nodes,  and  one  of  ICASE’s  graphics  workstations. 
Based  on  our  measurements,  the  Gigabit  Ethernet  based  on  Alteon’s  chipset  is  best  suited  for  bulk  data 
movement.  Latency  sensitive  parallel  codes  would  be  better  served  by  another  network  technology.  This 
was  confirmed  by  our  test  of  Giganet  cLAN  in  May,  which  also  demonstrated  a  strong  influence  of  the  MPI 
library  implementation  on  dual  CPU  nodes. 

The  Coral  cluster  usage  has  dramatically  increased  this  year.  In  August  of  2000,  Coral’s  existing  64 
compute  CPUs  delivered  15,890  CPU*hours  of  computing.  We  will  continue  to  use  this  cluster  to  develop 
and  run  research  codes  of  interest  to  ICASE  and  NASA  Langley,  and  to  evaluate  price/performance  tradeoffs 
among  various  hardware,  software,  and  networking  configurations. 

Integrating  process  algebra  and  linear-time  temporal  logic 

Gerald  Lilttgen 

This  research  is  devoted  to  the  development  of  novel  techniques  for  heterogeneous  formal  specifications 
of  reactive  systems  such  as  avionics  systems  and  embedded  software.  Traditional  formal-methods  research 
in  this  area  has  focused  on  homogeneous  specification  methodologies  in  which  all  specifications  are  given  in 
the  same  notation.  In  practice,  however,  heterogeneous  methodologies  supporting  multi-paradigm  specifica¬ 
tions  are  desired  for  the  following  reasons.  First,  different  system  components  have  different  characteristics. 
Second,  components  of  real-world  systems  are  often  developed  within  different  teams.  Third,  system  require¬ 
ments,  which  are  often  stated  using  assertions  (e.g.,  via  temporal  logics),  should  be  refinable  into  system 
designs  that  are  operationally  concrete  (e.g.,  represented  in  process  algebra ). 

The  taken  approach  extends  the  process  algebra  Calculus  of  Communicating  Systems  (CCS)  by  two 
operators  that  permit  the  inclusion  of  Linear-time  Temporal  Logic  (LTL)  formulas:  (i)  embed{<j>),  where  <f)  is 
a  LTL  formula;  this  operator  facilitates  the  embedding  of  LTL  formulas  within  process-algebraic  terms;  and 
(ii)  constraint(P,(f)),  where  P  is  a  process  term  and  4>  is  a  LTL  formula;  this  operator  allows  one  to  restrict 
the  behavior  of  P  as  specified  by  4>.  The  resulting  mixed  language  was  intended  to  be  given  a  semantics  in 
terms  of  Biichi  testing ,  a  semantic  framework  which  was  previously  developed  under  ICASE/NASA  funding. 
In  the  context  of  compositionality,  however,  it  turned  out  that  Streett  automata  provide  a  notationally  more 
convenient,  albeit  equally  expressive,  semantic  domain  than  Biichi  automata.  Consequently,  our  Biichi- 
testing  theory  has  been  re-developed  for  Streett  automata,  and  the  resulting  framework  has  been  interfaced 
to  our  mixed  language  via  structural-operational  rules. 

Regarding  future  work,  we  plan  to  investigate  concurrency-theoretic  issues  of  our  heterogeneous  language. 
Also,  an  efficient  algorithm  for  determining  whether  one  specification  refines  another  within  our  semantic 
theory  should  be  developed. 

This  work  is  done  in  collaboration  with  Ranee  Cleaveland  (SUNY  at  Stony  Brook,  New  York)  and  is 
supported  by  NSF  Grant  CCR-9988489. 

Relating  asynchronous  systems  with  respect  to  speed 

Gerald  Liittgen 

Classical  process  algebras  model  concurrent  systems  whose  components  are  assumed  to  have  arbitrary 
relative  speeds.  To  consider  the  temporal  behavior  of  systems,  several  timed  process  algebras  have  been 
proposed.  Remarkably,  most  of  these  assume  a  synchronous  system  model  where  all  system  components  have 
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fixed  speeds.  The  objective  of  our  research  is  to  develop  a  timed  process  algebra  for  evaluating  the  temporal 
worst-case  efficiency  of  asynchronous  concurrent  systems,  in  which  every  system  component  potentially 
performs  at  a  different  speed,  but  guarantees  an  upper  time  bound  on  all  of  its  actions. 

Our  approach  extends  Milner’s  Calculus  of  Communicating  Systems  (CCS)  by  a  clock-prefix  operator. 
This  operator  allows  for  specifying  a  discrete  amount  of  time  a  process  is  being  permitted  to  delay  before 
it  must  engage  in  an  enabled  communication  with  its  environment.  Since  the  semantic  theory  of  CCS 
is  based  on  the  behavioral  equivalence  bisimulation,  we  aimed  at  refining  this  relation  to  a  preorder  that 
compares  functionally  equivalent  processes  with  respect  to  their  speed.  To  this  end,  we  studied  several 
possible  definitions  for  such  a  preorder,  all  of  which  are  intuitively  justifiable,  and  formally  proved  them  to 
coincide.  For  the  resulting  preorder,  we  characterized  the  fully  abstract  pre-congruence  contained  in  it  and 
axiomatized  this  pre-congruence  for  finite  sequential  processes.  Finally,  we  also  investigated  a  corresponding 
‘weak’  pre-congruence  which  abstracts  from  internal  computation.  Our  work  concludes  long  overdue  research 
in  process  algebra.  The  obtained  results  testify  to  the  elegance  of  bisimulation-based  behavioral  relations 
when  compared  to  related  preorders  built  on  top  of  DeNicola  and  Hennessy’s  testing  theory. 

Regarding  future  work,  we  plan  to  complete  the  concurrency-theoretic  study  of  our  faster-than  preorder 
and  to  publish  our  findings  in  a  technical  report. 

This  work  is  done  in  collaboration  with  Walter  Vogler  (University  of  Augsburg,  Germany). 

Saturation:  An  efficient  iteration  strategy  for  implicit  state-space  generators 

Gerald  Liittgen 

Many  state-of-the-art  verification  techniques  rely  on  the  automated  construction  of  the  reachable  state 
space  of  the  system  under  consideration.  In  previous  work  we  developed  an  algorithm  for  the  MDD-based 
construction  of  state  spaces  of  asynchronous  system  models,  such  as  Petri  nets.  By  avoiding  to  encode  a 
given  global  next-state  function  implicitly  as  MDD,  but  splitting  it  in  several  explicitly  represented  local 
next-state  functions  instead,  we  gained  the  freedom  to  choose  the  sequencing  of  firing  events,  which  controls 
the  fixed-point  iteration  calculating  the  desired  global  state  spaces.  The  objective  of  this  research  is  to 
develop  a  suitable  strategy  for  firing  events  which  increases  the  time  and  space  efficiency  of  state-space 
generation. 

Our  approach  is  based  on  the  observation  that  an  MDD  node  is  not  modified  any  more,  when  all  events 
affecting  the  considered  node  have  been  fired  exhaustively,  i.e.,  when  the  node  is  saturated.  Having  the  specific 
representation  of  state-spaces  using  MDDs  in  mind,  this  results  in  a  depth-first  strategy  for  saturating  MDD 
nodes.  The  intuitive  benefits  of  this  strategy  are  threefold.  First,  it  brings  MDD-nodes  early  into  their 
final  shapes.  This  eliminates  much  of  the  overhead  necessary  in  our  previous  algorithm.  For  example,  it 
makes  many  checks  determining  whether  some  event  is  enabled  at  a  given  snap-shot  during  the  algorithm’s 
execution  unnecessary.  Second,  it  permits  an  easy  cache  management  in  which  cache  entries  do  not  need 
to  be  invalidated,  as  they  either  concern  MDD  nodes  which  are  not  updated  any  further  or  are  guaranteed 
never  to  be  accessed  again.  Third,  it  allows  for  considering  event  locality  naturally,  thereby  eliminating  the 
need  for  accessing  nodes  rather  randomly  and,  thus,  making  the  corresponding  infrastructure  obsolete.  Our 
novel  iteration  strategy  leads  to  a  new  state-space  generation  algorithm  which  has  been  implemented  in  the 
Petri  net  tool  SMART.  The  algorithm’s  utility  has  been  proved  by  experimental  studies  which  show  that 
it  often  performs  several  orders  of  magnitude  faster  than  existing  algorithms.  Equally  important,  it  keeps 
the  peak  sizes  of  MDDs  usually  quite  close  to  their  final  sizes.  Our  results  will  soon  enable  the  verification 
of  larger  systems  in  less  time  and,  thus,  will  help  to  employ  automated  formal  verification  techniques  into 
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engineering  practice. 

Regarding  future  work,  we  plan  to  employ  our  idea  of  saturation  for  implementing  an  MDD-based  CTL 
model  checker  within  smart.  Moreover,  we  intend  to  investigate  whether  our  new  algorithm  is  suitable  for 
parallelization. 

This  work  is  done  in  collaboration  with  Gianfranco  Ciardo  and  Radu  Siminiceanu  (The  College  of 
William  &  Mary). 

Arcade:  A  distributed  computing  environment  for  ICASE 

Piyush  Mehrotra 

Distributed  heterogeneous  computing  is  being  increasingly  applied  to  a  variety  of  large  size  computational 
problems.  Such  computations,  for  example,  the  multidisciplinary  design  optimization  of  an  aircraft,  generally 
consists  of  multiple  heterogeneous  modules  interacting  with  each  other  to  solve  the  problem  at  hand.  These 
applications  are  generally  developed  by  a  team  in  which  each  discipline  is  the  responsibility  of  experts  in  the 
field.  The  objective  of  this  project  is  to  develop  a  GUI-based  environment  which  supports  the  multi-user 
design  of  such  applications  and  their  execution  and  monitoring  in  a  heterogeneous  environment  consisting 
of  a  network  of  workstations,  specialized  machines  and  parallel  architectures. 

We  have  been  implementing  a  Java-based  three-tier  prototype  system  which  supports  a  thin  client 
interface  for  the  design  and  execution  of  multi-module  codes.  The  middle  tier  consists  of  logic  to  process  the 
user  input  and  also  to  manage  the  resource  controllers  which  comprise  the  third  tier.  We  have  continued 
implementation  of  our  Jini-based  resource  monitoring  and  management  system.  The  system  is  being  designed 
to  be  modular  with  well-defined  interfaces  so  that  it  can  be  used  with  other  frameworks  also.  For  the 
specification  of  the  resources,  we  are  closely  following  the  standards  being  set  by  the  Grid  Forum  Information 
services  group  and  are  implementing  the  specification  using  XML.  Similarly  we  have  also  designed  an  XML- 
based  project  specification  language  such  that  both  the  work-flow  and  the  inter-module  data  flow  can  be  easily 
specified.  We  have  implemented  the  translators  for  converting  back  and  forth  between  the  XML  specification 
and  our  internal  Java-based  project  specification.  We  have  also  started  a  new  subproject  focused  on  using 
XML  for  managing  scientific  data  used  in  multi-module  applications.  This  approach  will  allow  meta-data 
describing  the  data  to  be  specified  using  XML  thus  allowing  the  descriptions  to  be  manipulated  directly  by 
machines.  Thus,  translations  between  different  forms  of  the  same  data,  e.g.,  cell-centered  vs.  vertex-centered 
unstructured  grid  data,  can  be  automatically  translated. 

In  the  future  we  will  continue  to  work  on  the  enhancements  outlined  above.  In  addition  we  will  enhance 
the  current  Arcade  monitoring  system  by  adding  more  visualization  tools.  We  will  also  add  support  for 
dynamic  steering  of  executing  simulation  codes. 

This  work  is  being  performed  in  collaboration  with  A.  Al-Theneyan  and  M.  Zubair  (Old  Dominion 
University). 

Languages  for  high  performance  and  distributed  computing 

Piyush  Mehrotra 

There  are  many  approaches  to  exploiting  the  power  of  parallel  and  distributed  computers.  Under  this 
project,  our  focus  is  to  evaluate  these  different  approaches,  proposing  extensions  and  new  compilation  tech¬ 
niques  where  appropriate. 

OpenMP  is  a  set  of  directives  extending  C,  C++  and  Fortran  which  provide  a  shared-memory  paral¬ 
lel  programming  model.  Current  parallel  architectures  are  built  by  interconnecting  nodes  which  internally 
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provide  true-shared  memory  across  a  small  number  of  processors.  Both  hardware-  and  software-based  ap¬ 
proaches  are  used  to  provide  a  shared  address  space  across  the  physically  distributed  memories  of  the  nodes 
in  larger  systems.  OpenMP  provides  an  easy  and  incremental  approach  for  small  scale  shared  memory  sys¬ 
tems.  However,  controlling  and  exploiting  data  locality  becomes  an  issue  as  the  latency  for  data  transfers 
across  the  nodes  increases  for  larger  systems.  We  have  designed  a  set  of  directives  extending  OpenMP  which 
allow  users  to  directly  express  the  data-locality  of  their  programs.  Such  directives  can  then  be  translated 
by  a  compiler  to  take  advantage  of  the  underlying  target  architecture.  In  order  to  test  our  ideas,  we  hand- 
translated  directives  for  two  different  programs  targeting  two  different  architectures:  a)  a  simple  Jacobi 
program  on  an  IBM  SP2  using  the  Treadmarks  software  based  DSM  system  from  Rice  University,  and  b)  a 
program  implementing  lattice  Boltzmann  equations  for  solving  flows  (provided  by  L.-S.  Luo  and  group)  on 
the  SGI  Origin  2000.  For  both  the  programs  we  wrote  at  least  two  versions,  a  pure  OpenMP  version  and  one 
in  which  most  of  the  data  is  declared  private  with  only  buffers  being  shared  between  the  threads.  In  both 
cases  we  found  that  the  second  version  gave  much  better  performance  reinforcing  the  fact  that  data-locality 
is  important  even  on  shared  memory  systems. 

We  plan  to  continue  our  investigations  on  providing  data-locality  in  OpenMP  programs,  including  de¬ 
signing  and  implementing  compiler  transformations  required  to  translate  the  distribution  directives  so  as  to 
exploit  the  underlying  architectures. 

This  work  is  being  done  in  collaboration  with  B.  Chapman,  A.  Patil  and  A.  Prabhakar  (University  of 
Houston). 

On  the  initial  condition  for  LBGK  scheme  for  unsteady  flow  simulations 

Renwei  Mei 

In  the  computational  simulation  for  unsteady  flow  problems  using  the  method  of  lattice  Boltzmann 
equation,  it  is  critically  important  to  have  a  good  initial  condition  for  the  particle  distribution  functions. 
Otherwise,  the  error  in  the  initial  condition  could  contaminate  the  whole  solution.  A  typical  approach  is  to 
set  the  initial  particle  distribution  function  fi’s  equal  to  the  equilibrium  particle  distribution  function  based 
on  the  known  velocity  and  pressure  field.  This  method  is,  however,  shown  to  have  larger  initial  error. 

We  propose  an  initialization  scheme  to  obtain  improved  initial  condition  for  fi’s.  Let  the  equilibrium 
function  based  on  the  initial  velocity  and  density  field  be  /(e<?’0).  This  initialization  scheme  let  the  initial 
field  for  fi’s  to  stream  and  collide  based  on  the  /(eg’0).  During  this  initialization  procedure,  macroscopic 
flow  variables  remain  unchanged.  One  of  the  advantages  of  this  initialization  scheme  is  that  one  can  use  the 
existing  code  to  obtain  easily  an  improved  initial  condition  for  fi’s.  The  improvement  is  measured  for  velocity 
and  stress  field  in  two  cases:  decaying  Taylor  vortex  flow,  and  decaying  shear  flow.  Significant  improvement 
is  obtained  for  low  viscosity  (or  high  Reynolds  number). 

The  results  of  this  work  will  be  written  for  an  archival  journal. 

This  work  is  an  collaborative  effort  with  Dr.  Li-Shi  Luo,  Profs.  Domonique  d’Humieres  and  Lian-Ping 
Wang. 

Formal  analysis  of  the  AILS  alerting  algorithm 

Cesar  Munoz 

The  Airborne  Information  for  Lateral  Spacing  (AILS)  is  a  project  being  conducted  at  NASA  Langley 
Research  Center.  Its  global  objective  is  to  increase  airport  efficiency  by  enabling  approaches  to  closely  spaced 
parallel  runways  in  Instrument  Meteorological  Conditions.  One  of  the  most  critical  components  of  the  AILS 
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concept  is  the  alerting  algorithm.  It  provides  situational  awareness  to  pilots  involved  in  a  parallel  approach. 
The  aim  of  this  work  is  to  conduct  a  formal  analysis  of  the  AILS  alerting  algorithm  in  order  to  discover 
possible  misbehaviors  that  have  not  been  detected  during  simulation  and  testing. 

An  abstraction  of  the  AILS  alerting  algorithm  has  been  formalized  in  the  specification  language  of  the 
PVS  verification  tool.  Several  models  of  collision  trajectories,  which  vary  from  a  pure  discrete  model  to  a 
completely  continuous  one,  have  been  considered.  Discrete  models  allow  the  use  of  efficient  state  exploration 
techniques  to  validate  formal  properties.  However,  they  have  proven  to  be  too  inaccurate.  An  accurate 
continuous  model  of  trajectories  has  been  refined  by  gradually  removing  discrete  parameters  of  the  discrete 
models.  As  result  of  the  formal  analysis,  maximum  and  minimum  times  when  an  alarm  will  first  sound  prior 
to  a  collision  have  been  found. 

Techniques  and  models  developed  in  the  formal  analysis  of  the  AILS  concept  will  be  used  to  study  a  more 
general  concept  called  Distributed  Air-Ground  Traffic  Management  (DAG/TM).  DAG/TM  is  a  gate-to-gate 
concept  being  designed  as  part  of  the  NASA  Aviation  System  Capacity  Program  to  improve  the  overall 
performance  of  the  National  Airspace  System. 

This  work  is  done  in  collaboration  with  Victor  Carreno  (NASA  Langley). 

Verification  of  Java  byte-code 

Cesar  Munoz 

A  key  aspect  of  the  Java  language  is  the  definition  of  a  standard  execution  model,  called  the  Java  Virtual 
Machine  (JVM),  which  executes  machine-independent  byte-code.  In  this  work  we  apply  Formal  Methods 
techniques  to  formally  analyze  multi-thread  executions  of  Java  programs. 

A  formal  model  of  concurrent  Java  byte-code  based  on  transition  systems  has  been  developed.  Given 
a  particular  Java  program,  a  specification  of  its  execution  is  written  as  an  infinite  state  transition  system. 
This  model,  that  we  write  in  the  SAL  specification  language,  can  then  be  simplified  into  a  finite  system  by 
using  predicate  abstraction  techniques.  Properties  such  as  mutual-exclusion  have  been  verified  for  classical 
examples  of  Java  programs. 

Aspects  not  considered  in  the  first  approach,  such  as  exceptions  and  a  richer  set  of  instructions  of  the 
Java  byte-code,  will  be  considered  in  future  research. 

This  work  is  done  in  collaboration  with  Peter  Habermehl  (University  of  Paris  7). 

External  memory  sparse  direct  solvers 

Alex  Pothen 

We  are  developing  external  memory  algorithms  and  software  for  solving  large,  sparse  systems  of  equations 
by  means  of  direct  solvers.  Such  methods  will  enable  sparse  direct  solvers  to  make  effective  use  of  the 
Terabytes  of  storage  available  on  serial  computers;  parallel  algorithms  running  on  PC  clusters  like  Coral  and 
Teraflop  parallel  computers  with  multiple  levels  of  memory  hierarchy  will  also  benefit. 

We  formalize  two  problems  for  external  memory  sparse  matrix  factorizations:  minimizing  the  core 
memory  required  in  a  read-once/write-once  model,  and  minimizing  the  data  movement  needed  in  a  read- 
many /write-many  model.  We  compute  bounds  on  these  quantities  for  sparse  problems  whose  computational 
graphs  have  good  separators.  We  show  that  a  data  structure  called  the  elimination  tree  determines  the 
core  memory  requirements  and  the  data  movement  costs  of  a  sparse  factorization  algorithm.  The  relative 
performance  of  the  three  commonly  used  variants  of  the  factorization  algorithm,  viz.  left-looking,  right¬ 
looking,  and  multifrontal,  can  vary  a  great  deal  for  unbalanced  elimination  trees  that  occur  in  mathematical 
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programming  and  related  applications.  We  are  designing  fast  simulators  that  compute  the  data  movement 
costs  of  these  algorithms  in  time  proportional  to  the  size  of  the  factors  rather  than  the  flops  required  to 
compute  them. 

We  also  propose  to  extend  OBLIO,  our  object-oriented  sparse  direct  solver  library,  with  serial  and 
parallel  implementations  of  external  memory  solvers. 

This  is  joint  work  with  Florin  Dobrian  (Old  Dominion  University). 

Scalable  parallel  algorithms  for  incomplete  factorization  preconditioning 

Alex  Pothen 

The  parallel  computation  of  robust  preconditioners  is  necessary  for  solving  large  systems  of  equations 
iteratively  in  scientific  computing.  We  are  developing  scalable  parallel  algorithms  and  software  that  can 
compute  incomplete  factorization  preconditioners  for  high  levels  of  fill. 

We  partition  the  adjacency  graph  of  the  coefficient  matrix  into  subgraphs  of  roughly  equal  sizes  such 
that  each  subgraph  has  few  boundary  nodes  relative  to  the  number  of  interior  nodes.  We  map  the  subgraphs 
to  processors,  form  a  subdomain  interconnection  graph,  and  order  the  subdomains  by  a  graph  coloring  to 
reduce  global  dependences.  On  each  subdomain,  we  locally  reorder  the  interior  vertices  before  the  boundary 
vertices.  This  reordering  limits  the  fill  that  joins  a  subgraph  on  one  processor  to  a  subgraph  on  another,  and 
enhances  the  concurrency  in  the  computation.  The  preconditioner  computation  takes  place  in  two  phases:  in 
the  first  phase,  each  processor  computes  the  rows  of  the  preconditioner  corresponding  to  the  interior  vertices 
of  their  subdomains.  In  the  second  phase,  the  rows  corresponding  to  the  boundary  nodes  are  computed. 

We  have  reported  results  on  problem  sizes  of  up  to  20  million  rows  and  columns  on  up  to  216  processors  of 
the  SGI  Origin,  the  Cray  T3E,  the  Sun  HPC  10000,  and  Coral,  the  ICASE  Beowulf  cluster.  Our  results  show 
that  the  algorithm  is  scalable,  in  agreement  with  our  analysis  of  the  algorithm.  We  are  currently  implementing 
our  parallel  algorithm  on  parallel  computers  with  Teraflop  performances.  We  are  also  developing  multi¬ 
threaded  implementations  on  computers  with  shared  address-space  programming  models. 

This  is  joint  work  with  David  Hysom  (Lawrence  Livermore  National  Laboratory). 

Enhancing  the  performance  of  unstructured  mesh  codes 

Alex  Pothen 

Irregular  computations  such  as  unstructured  mesh  algorithms  or  sparse  matrix  algorithms  achieve  only 
a  small  fraction  of  the  available  peak  performance  on  modern  microprocessors.  Understanding  the  reasons 
for  this  low  performance  is  difficult  on  modern  microprocessors  since  they  issue  multiple  instructions  in 
each  clock  cycle,  and  permit  them  to  execute  out  of  issue-order  by  exploiting  multiple  functional  units  and 
multiple  levels  of  caches.  We  have  initiated  a  study  measuring  the  performance  of  these  codes  using  hardware 
performance  counters  and  simulation  tools  to  first  understand  and  then  improve  the  performance. 

We  have  studied  the  performance  of  the  kernel  of  an  unstructured  mesh  code  for  solving  Euler’s  equations. 
By  measuring  various  performance  metrics,  we  find  that  on  an  SGI  architecture  the  performance  is  limited 
primarily  by  the  number  of  loads  and  stores  that  can  be  issued  in  a  clock  cycle,  and  that  the  level  1  and 
level  2  caches  play  a  less  important  role  in  determining  performance  on  the  SGI.  We  have  also  studied  the 
role  of  scheduling  data  accesses  and  computations  by  reordering  the  irregular  data  structures.  We  have 
experimented  with  ordering  algorithms  such  as  nested  dissection,  Cuthill-McKee,  and  space-filling  curve 
based  orderings.  These  improve  performance  for  different  reasons:  some  of  them  improve  spatial  locality, 
while  others  improve  temporal  locality. 
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We  are  continuing  our  work  to  find  general  principles  for  improving  the  performance  of  irregular  com¬ 
putations.  This  is  joint  work  with  Jinghua  Fu  (Old  Dominion  University),  and  Dimitri  Mavriplis  (ICASE). 

Numerical  modeling  of  sound  radiation  and  propagation  in  a  non-uniform  flow 

Alex  Povitsky 

The  propagation  of  disturbances  in  a  non-uniform  mean  flow  in  the  presence  of  a  reflecting  wall  is 
investigated  by  high-order  numerical  scheme  using  compact  approximation  of  spatial  derivatives.  Simulations 
were  performed  on  serial  and  parallel  computers  for  three-dimensional  flows.  Stagnation  flows  and  flow 
around  a  cylinder  are  taken  as  prototypes  of  real  world  flows  with  strong  pressure  gradients.  Monopole  and 
dipole  acoustic,  vortical  and  entropy  pulses  are  embedded  in  an  incompressible  inviscid  background  flow. 
The  intensity  and  directivity  of  near-field  wave  patterns  appear  to  be  quite  different  from  the  benchmark 
solutions  obtained  in  a  uniform  flow  or  in  a  static  environment. 

Near-wall  stretching  of  a  vortical  monopole  pulse  by  the  stagnation  flow  leads  to  formation  of  a  wave 
with  two  maximums  apart  from  the  centerline  whereas  a  vortical  dipole  evolves  into  a  three-spot  wave 
with  maximum  acoustic  pressure  at  the  centerline.  For  an  entropy  pulse,  the  baroclinically  induced  vortex 
pair  produces  a  sound  wave  similar  to  that  originated  from  a  vorticity  dipole.  Initial  superposition  of 
vorticity  and  entropy  pulses  forms  non-symmetrical  wave  patterns.  For  an  acoustic  wave  in  the  flow  around 
a  cylinder,  the  observed  mean  acoustic  pressure  is  approximately  doubled  (upstream  pulse  position)  and 
halved  (downstream  pulse  position)  in  comparison  with  the  benchmark.  Angular  redistribution  of  acoustic 
energy  and  its  components  and  energy  transfer  between  the  mean  flow  and  the  wave  are  discussed  in  ICASE 
Report  No.  2000-35. 

Our  future  plans  include  the  study  of  sound  generation  and  propagation  in  non-uniform  compressible 
flows  and  flows  with  vortex  shedding  (cavity  noise).  To  get  realistic  distribution  of  acoustic  sources  in  a 
turbulent  flow,  we  plan  to  use  the  LES  approach. 

Optimization  of  jet  mixing  reactor  for  production  of  carbon  nanotubes 

A  lex  Povitsky 

The  high-pressure  jet  mixing  reactor  has  been  proposed  for  the  industrial  production  of  fullerene  carbon 
nanotubes  by  the  Nanotechnology  Lab  at  Rice  University.  The  C6o  molecules  combined  in  nanotubes  have 
unusual  and  very  useful  electronic  and  mechanical  properties.  However,  the  mass  production  of  nanotubes 
still  presents  a  challenge.  The  single-walled  carbon  nanotubes  may  be  produced  from  CO  with  the  help  of 
a  catalyst  whose  rapid  and  uniform  heating  is  critical  for  feasibility  of  this  technology.  Catalyst  particles 
get  into  the  reaction  chamber  via  the  central  CO  jet.  The  particles  must  be  heated  from  200° C  (central  jet 
temperature)  to  temperatures  close  to  1000°C  as  rapidly  as  possible,  since  slowly  heated  catalytic  particles 
quickly  loose  their  catalytic  ability.  To  increase  mixing,  hot  turbulent  jets  issuing  from  several  nozzles 
surrounding  the  central  nozzle  are  used.  The  goal  of  our  numerical  modeling  and  optimization  of  the  reactor 
is  to  obtain  a  uniformly  high  temperature  of  the  catalyst. 

The  maximum  strength  of  each  peripheral  jet  and  the  constraint  on  a  total  mass  flux  of  peripheral 
jets  suggests  a  configuration  of  two  peripheral  jets  with  high  intersection  angle  between  the  central  and 
peripheral  jets.  The  optimal  distance  between  peripheral  nozzles  and  the  central  jet  results  from  the  tradeoff 
between  the  strength  of  a  peripheral  jet  while  it  hits  the  central  jet  and  the  co-flowing  effect  delaying  mixing 
of  the  central  jet.  Then,  the  FLUENT  CFD  software  and  GAMBIT  grid  generator  were  used  in  a  detailed 
simulation  of  two-phase  flow  in  the  reactor.  In  addition  to  solving  transport  equations  for  the  continuous 
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phase,  FLUENT  has  an  option  to  simulate  a  discrete  second  phase  in  a  Lagrangian  frame  of  reference.  The 
analysis  of  velocity  field  and  particle  trajectories  shows  that  strong  recirculation  zone  at  the  centerline  may 
lead  to  undesirable  clustering  of  catalyst  particles.  Thus,  the  Eulerian  optimization  of  reactor  should  be 
refined  by  the  Lagrangian  simulation  of  the  discrete  phase. 

In  our  future  research,  the  refined  models  of  turbulence  (as  the  RNG  k  —  e)  and  a  precise  thermal 
radiation  model  should  be  evaluated  to  predict  a  correct  rate  of  chemical  reactions  leading  to  formation  of 
nano-tubes  in  presence  of  the  catalyst.  Dr.  Carl  Scott  (NASA  Johnson  Space  Center)  applied  his  chemical 
model  of  the  catalyst  behavior  and  CHEMKIN  software  to  the  particle  trajectories  obtained  by  our  CFD 
computations  and  his  results  are  indicative  of  the  time  it  takes  to  produce  iron  clusters  as  large  as  20  or  40 
atoms.  To  verify  our  results  by  physical  experiments,  we  collaborate  with  the  Nanotechnology  Laboratory 
at  Rice  University  (Prof.  Richard  Smalley). 

Solving  radiation  transport  equations  with  the  multigrid  method 

Linda  Stals 

Radiation  transport  equations  arise  in  the  study  of  many  different  fields  such  as  combustion,  astrophysics 
and  hypersonic  flow.  The  solution  of  these  equations  presents  interesting  challenges  due  to  large  jumps  in  the 
coefficients  and  strong  nonlinearities.  In  this  project  we  compare  the  efficiency  of  several  different  solution 
techniques.  Specifically,  we  focus  on  the  performance  of  an  inexact  Newton  multigrid  scheme  and  compare 
it  to  the  Full  Approximation  Scheme  (FAS). 

We  have  currently  investigated  the  case  where  all  temperatures  are  in  equilibrium  and  the  model  reduces 
to  a  single  equation.  This  equation  contains  many  of  the  features  inherent  in  the  full  model  such  as  the 
strong  nonlinearities  and  large  jumps  in  the  coefficients.  We  found  that  the  Newton-Multigrid  method 
worked  better  because  it  did  not  need  to  calculate  the  diffusion  coefficient  as  often.  Such  an  observation  is  of 
interest  because  including  more  physics  in  the  model  will  increase  the  complexity  of  the  diffusion  coefficient 
even  further  and  thus  the  Newton-Multigrid  method  is  the  method  of  choice.  A  flux  limiter  can  also  be 
included  in  the  equation  which  changes  its  behavior  from  being  locally  parabolic  to  locally  hyperbolic.  We 
found  that  in  such  a  case  the  use  of  adaptive  refinement  helped  the  solvers  to  better  handle  the  steep  wave 
front. 

In  future  work  we  want  to  include  more  physics  in  the  equations.  The  PFs  for  the  ASCI  II  project  are 
David  Keyes  and  Alex  Pothen  (Old  Dominion  University  and  ICASE)  and  Dimitri  Mavriplis  (ICASE). 
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FLUID  MECHANICS 


Receptivity  and  transition  in  laminar  boundary  layers 

W.O.  Criminale  and  D.G.  Lasseigne 

This  work  has  been  the  principal  area  of  research  and  a  comprehensive  and  thorough  treatise  based  on 
first  principles  has  been  made.  This  report  is  entitled  “A  new  paradigm  for  investigating  receptivity  and 
bypass  dynamics  in  laminar  boundary  layers,”  by  W.O.  Criminale,  R.D.  Joslin,  and  T.L.  Jackson  and  will 
be  submitted  to  the  Journal  of  Fluid  Mechanics  for  publication. 

Flow  control  and  drag  reduction 

W.  0.  Criminale  and  D.  G.  Lasseigne 

This  problem  is  a  natural  extension  of  the  work  done  on  boundary  layer  receptivity.  More  specifically,  by 
using  the  framework  of  initial- value  problems  for  vorticity  in  the  boundary  layer,  it  is  possible  to  determine 
the  full  range  of  the  temporal  dynamics.  As  a  result  it  is  then  possible  to  design  feedback  in  such  a  way 
as  to  suppress  any  damaging  response  in  this  system.  The  question  remaining  now  is  how  to  transpose 
such  results  into  physically  realizeable  mechanisms  for  the  boundary  layer.  Not  only  is  the  analysis  able 
to  investigate  such  behavior,  but  it  has  led  to  a  more  rational  method  for  analytically  solving  the  stability 
equations  in  this  field.  A  second  report  is  in  progress  that  deals  with  the  general  technique  (W.O.  Criminale 
and  D.G.  Lasseigne). 

Laser-induced  thermal  acoustics  for  flow  diagnostics 

Roger  C .  Hart 

Off-body  nonintrusive  measurement  of  flow  velocities  is  currently  routinely  done  with  laser  doppler 
velocimetry  (LDV)  and  particle  interval  velocimetry  (PIV).  Both  techniques  require  seeding  the  flow  with 
silica  particles  or  droplets  of  liquid  which  serve  as  scatterers  for  the  incident  laser  light.  Seeding  is  undesirable 
in  some  circumstances  (e.g.,  highly  polished  high-Reynolds  number  facilities  or  supersonic  facilities,  where 
the  seed  may  literally  sandblast  models  and  tunnel  surfaces),  and  impossible  in  others  (there  may  be  locations 
on  models,  behind  rearward  facing  steps  for  instance,  where  you  simply  cannot  get  the  seed  to  go).  Laser- 
induced  thermal  acoustics  (LITA)  is  a  relatively  new  technique  being  developed  at  Langley  and  a  few  other 
places  which  allows  truly  seedless  point  velocimetry,  and  additionally  measures  temperature  and  pressure. 
Our  goal  is  to  explore  and  demonstrate  the  potential  of  this  method  by  a  combination  of  laboratory  research 
and  tunnel  measurements. 

In  LITA  a  powerful  pulsed  laser  creates  a  pair  of  counterpropagating  acoustic  wave  packets  at  the 
measurement  point.  A  second  CW  or  long-pulse  probe  laser  illuminates  the  wave  packets  at  the  Bragg 
diffraction  angle.  Doppler-shifted  light  diffracted  from  each  packet  is  mixed  on  a  detector  with  a  small 
fraction  of  the  unshifted  probe  beam.  Analysis  of  the  three  resulting  beat  frequencies  allows  determination 
of  both  a  component  of  flow  velocity  and  the  speed  of  sound  (and  thus  temperature).  Additionally,  pressure 
can  be  found  by  analysis  of  the  temporal  decay  of  the  signal.  A  measurement  is  made  in  a  few  microseconds 
over  a  volume  of  100’s  of  microns  (transverse)  by  ~  a  centimeter  (longitudinally) .  A  first-generation  prototype 
LITA  apparatus  has  been  constructed  using  compact,  fieldable,  “user-friendly”  lasers  and  employing  a  novel 
technique  for  performing  the  heterodyne  signal  detection  that  is  vastly  more  stable  than  our  earlier  efforts. 
The  apparatus  consists  of  transmitter  and  receiver  units,  each  with  a  volume  of  about  three  cubic  feet  and 
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weighing  less  than  100  pounds  each,  plus  one  or  two  racks  of  electronics,  power  supplies,  computers,  and 
other  ancillary  equipment.  Development  of  this  instrument  is  continuing,  but  currently  we  observe  one-sigma 
single-shot  uncertainties  of  0.7  to  1.0  m/s  over  the  range  of  8  m/s  to  35  m/s.  Performance  is  limited  largely 
by  signal/noise  level. 

We  currently  expect  to  begin  a  40  day  test /demonstration  of  LITA  at  BART  during  the  second  week  of 
October  2000.  We  will  compare  LITA  velocity  measurements  with  pitot  tube  measurements  of  free  stream 
velocity,  and  then  attempt  to  map  the  separation,  reattachment,  and  recirculation  behind  a  rearward-facing 
step. 

This  research  was  conducted  in  collaboration  with  R.J.  Balia  and  G.C.  Herring  (NASA  Langley). 

The  probability  density  function  (PDF)  approach  to  turbulence  modeling 

Guowei  He 

The  primary  motivation  of  this  research  is  to  develop  a  PDF  method  for  turbulence  modeling.  In 
comparison  with  the  moment  method,  the  PDF  method  takes  advantage  to  exactly  model  nonlinear  terms, 
but  needs  to  model  molecular  diffusion  and  pressure.  We  will  develop  a  mapping  closure  approximation 
(MCA)  for  scalar  diffusion  in  turbulent  shear  flows.  The  second  motivation  of  this  research  is  to  investigate 
the  effect  of  shear  on  small-scale  statistics  in  turbulence.  The  current  turbulence  models  are  mainly  based 
on  the  assumption  of  isotropy.  We  will  develop  a  model  to  incorporate  small  scale  anisotropy. 

We  have  developed  a  two-dimensional  MCA  approach  to  calculate  the  PDFs  of  scalar  gradients  in 
turbulent  shear  flows,  where  velocity  is  designed  to  be  a  superposition  of  isotropic  Gaussian  random  filed 
and  a  directed  shear  flow.  We  derived  and  analyzed  the  governing  equation  of  mapping  function  by  exactly 
modeling  “rapid-change”  shear  flows.  The  result  obtained  shows  that  shear  does  enhance  more  intense 
intermittency  of  scalar  gradient  in  the  direction  of  shear  than  the  one  transversal  to  the  direction  of  shear 
and  induces  difference  of  their  PDFs.  Therefore,  scalar  field  in  shear  flow  is  anisotropic  at  the  level  of 
dissipation. 

We  plan  to  numerically  verify  the  above  theoretical  results. 

This  work  is  accomplished  in  collaboration  with  Dr.  R.  Rubinstein. 

Direct  numerical  simulation  (DNS)  of  isotropic  homogeneous  turbulence 

Guowei  He 

In  a  collaborative  effort  with  L.P.  Wang  (University  of  Delaware),  we  have  developed  a  spectrum  code 
specially  for  the  cluster  computer  at  ICASE  for  DNS  of  isotropic  homogeneous  turbulence.  Several  test  cases 
have  been  successfully  run  at  the  grid  size  643,  1283  and  2563.  The  largest  case  at  the  grid  size  5123  allowed 
by  current  capacity  of  the  cluster  computer  is  running  for  the  largest  database  which  can  resolve  small  scale 
fluctuations  in  turbulence. 

The  database  from  DNS  can  be  used  to  investigate  the  mechanism  and  modeling  of  turbulence.  In 
collaboration  with  R.  Rubinstein,  we  have  calculated  the  time  correlations  of  velocity  fields  from  the  DNS 
database  and  will  compare  them  with  the  results  obtained  from  LES  subgrid  models.  The  comparison  allows 
us  to  estimate  the  predictability  of  the  current  LES  subgrid  models,  which  has  not  been  investigated  so  far. 

We  plan  to  apply  the  comparison  methodology  to  investigate  a  variety  of  turbulence  models. 
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Eigen-solution  analysis  of  the  discontinuous  Galerkin  method 

Fang  Hu 

It  is  well  known  that  for  a  discontinuous  Galerkin  scheme  employing  basis  polynomials  of  order  p,  the 
rate  of  convergence  has  been  shown  to  be  hp+1/2  in  general  and  hp+1  in  some  special  cases.  In  contrast  to 
the  many  studies  on  the  convergence  rates  there  have  been  relatively  few  works  on  the  wave  propagation 
properties  of  DGM,  especially  with  non-uniform  grids.  Our  work  is  motivated  by  the  need  to  understand 
DGM  with  non-uniform  grids. 

A  detailed  study  on  spatially  propagating  waves  in  a  discontinuous  Galerkin  scheme  has  been  carried  out. 
We  are  especially  interested  in  wave  propagation  properties  when  there  is  an  abrupt  change  of  element  sizes, 
as  well  as  the  influences  the  flux  formula  has  on  these  properties.  By  analyzing  the  semi-discrete  system,  it 
is  found  that,  for  any  given  order  of  the  basis  functions,  there  are  only,  at  most,  two  spatially  propagating 
numerical  wave  modes  for  each  physical  wave  of  the  PDE.  One  of  the  modes  can  accurately  represent  the 
physical  wave  and  the  other  is  spurious.  The  directions  of  propagation  of  these  two  numerical  modes  are 
opposite.  In  most  cases,  the  spurious  mode  has  a  large  damping  rate.  When  full  (exact)  upwind  flux  formula 
is  used,  however,  the  spurious  mode  is  non-existent.  It  is  also  shown  analytically  that,  for  the  physically 
accurate  mode,  the  numerical  dispersion  relation  is  accurate  to  order  2p  +  2  where  p  is  the  highest  order  of 
the  basis  polynomials.  The  results  of  the  eigen-solutions  are  then  used  to  study  the  propagation  of  waves 
through  an  interface  of  abrupt  change  of  element  sizes.  Closed  form  numerical  reflection  and  transmission 
coefficients  are  derived  and  analyzed. 

Works  on  2-D  and  3-D  equations  will  be  carried  out  in  the  future. 

This  work  is  in  collaboration  with  H.L.  Atkins  (NASA  Langley). 

Dispersion,  dissipation,  isotropy,  Galilean  invariance,  and  stability  of  the  lattice  Boltzmann 
equation 

Li-Shi  Luo 

Like  any  numerical  schemes  to  solve  partial  differential  equations,  the  lattice  Boltzmann  equation  (LBE) 
has  a  number  of  numerical  artifacts  due  to  the  discretization.  These  artifacts  are  yet  to  be  fully  understood. 
Our  objective  is  to  analyze  the  artifacts  of  the  lattice  Boltzmann  equation,  such  as  its  dissipation  for  non¬ 
zero  wavenumber  (hyper- viscosities),  its  isotropy  (wave-number  dependence  of  the  transport  coefficients), 
and  non-Galilean  invariant  effects  due  to  discretization. 

We  construct  a  2D  generalized  lattice  Bhatnagar-Gross-Krook  (LBGK)  Boltzmann  model  with  multiple 
relaxation  parameters  in  moment  space.  We  analyze  the  linearized  dispersion  equation  of  the  generalized 
LBGK  model  to  obtain  its  generalized  hydrodynamics. 

We  have  obtained  the  wave-number  dependence  of  transport  coefficients  theoretically,  and  verified  the 
results  by  various  numerical  simulations.  By  introducing  multiple  relaxation  parameters  in  the  collision 
process  of  the  LBE  model,  we  can  drastically  improve  the  numerical  stability  of  the  LBE  method.  We  have 
also  quantitatively  analyzed  the  impact  of  interpolations  to  the  LBE  method.  This  work  has  been  published 
in  Physical  Review  E. 

We  will  extend  the  model  to  3D  and  for  complex  fluids. 

This  work  is  accomplished  in  collaboration  with  Prof.  P.  Lallemand  (Laboratoire  C.N.R.S.-A.S.C.L, 
Prance),  and  it  is  partially  supported  by  the  program  of  Innovative  Algorithms  for  Aerospace  Engineering 
Analysis  and  Optimization  at  NASA  LaRC. 
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Stochastic  description  of  transitional  flows 

Robert  Rubinstein 

There  is  still  no  well-established  procedure  for  incorporating  transition  in  turbulence  calculations.  While 
aerodynamic  flows  can  be  computed  successfully  using  any  of  several  different  turbulence  models  if  the 
transition  location  is  prescribed  in  advance,  no  single  turbulence  model  can  reliably  predict  transition.  If 
transition  is  computed  incorrectly,  the  entire  flow  calculation  is  generally  unsatisfactory. 

The  analysis  of  the  growth  of  subharmonics  due  to  three- wave  interactions  was  applied  to  an  experiment 
performed  by  Kachanov  and  collaborators.  Analysis  of  this  experiment  had  suggested  an  important  effect  of 
initial  phase  correlations;  the  present  work  quantifies  this  effect  by  evaluating  the  joint  probability  density 
of  modal  amplitudes  and  phases.  The  probability  that  the  subharmonic  amplitude  exceeds  the  primary  has 
also  been  evaluated  as  a  function  of  local  Reynolds  number;  it  is  possible  that  this  condition  can  serve  as  a 
simple  transition  criterion. 

Work  on  the  stochastic  formulation  of  the  parabolized  stability  equations  is  continuing. 

This  work  was  performed  in  collaboration  with  Meelan  Choudhari  (NASA  Langley). 
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STRUCTURES  AND  MATERIALS 


Decohesion  elements  for  the  simulation  of  progressive  delamination  in  composites 

Pedro  P.  Camanho 

The  fracture  process  of  high  performance  composite  laminates  is  quite  complex,  involving  both  intralam¬ 
inar  damage  mechanisms  (e.g.,  matrix  cracking,  fibre  fracture)  and  interlaminar  damage  (delamination).  Al¬ 
though  progress  has  been  made  lately  in  the  development  of  accurate  tools  for  the  prediction  of  intralaminar 
damage  growth,  similar  tools  for  delamination  are  still  difficult  to  implement,  and  thus  delamination  is  gen¬ 
erally  not  considered.  Without  delamination,  the  predictive  capabilities  of  progressive  failure  analyses  will 
remain  limited.  The  analysis  of  delamination  is  commonly  divided  into  the  study  of  the  initiation  and  the 
analysis  of  the  propagation  of  an  already  initiated  area.  Delamination  initiation  analysis  are  usually  based  on 
stresses  and  use  criteria  such  as  the  quadratic  interaction  of  the  interlaminar  stresses  in  conjunction  with  a 
characteristic  distance.  This  distance  is  a  function  of  geometry  and  material  properties,  so  its  determination 
always  requires  extensive  testing. 

One  approach  to  simulate  both  delamination  onset  and  growth  in  composite  materials  is  the  use  of 
interfacial  decohesion  elements  placed  between  the  composite  layers.  This  approach  is  based  on  a  Dudgale- 
Barenblatt  type  cohesive  zone.  Cohesive  zone  models  are  particularly  attractive  when  interfacial  strengths 
are  relatively  weak  when  compared  with  the  adjoining  material,  as  in  the  case  of  composite  laminates.  The 
proposed  constitutive  equations  for  the  interface  are  phenomenological  mechanical  relations  between  the 
tractions  and  interfacial  separations.  The  work  of  normal  and  tangential  separation  can  be  related  to  the 
critical  values  of  energy  release  rates.  In  fact,  it  can  be  shown  that  the  cohesive  zone  formulation  is  identical 
to  Griffith’s  theory  of  fracture. 

In  order  to  predict  delamination  onset  and  growth,  an  eight-node  decohesion  element  with  interpolating 
shape  functions  for  the  top  and  bottom  faces  compatible  with  those  of  the  elements  being  connected  will 
be  developed.  The  material  response  built  into  the  element  represents  damage  using  Barenblatt’s  model  in 
which  a  cohesive  zone  ahead  of  crack  tip  is  used  to  predict  delamination  growth.  A  bilinear  strain  softening 
constitutive  equation  will  be  implemented  such  that  the  element  does  not  carry  any  tensile  or  shear  loads 
after  the  interfacial  fracture  energy  has  been  consumed.  Such  an  element  will  be  used  to  model  the  interface 
between  sublaminates  or  between  bonded  components. 

Update  of  Dr.  Jer-Nan  Juang’s  system  observer /controller  identification  toolbox 

Joshua  Grohs 

Due  to  the  increasing  use  of  Dr.  Juang’s  System  Observer/ Controller  Identification  Toolbox  (SOCIT) 
the  necessity  of  expanding  its  compatibility  and  support  became  imperative.  With  the  current  reliance  on 
the  internet,  Dr.  Juang  felt  it  necessary  to  provide  online  support  for  his  toolbox.  In  addition  to  text-based 
support,  he  also  believed  that  it  needed  to  support  a  couple  of  different  platforms,  therefore  translation  of 
the  toolbox  began.  These  projects  were  the  focus  of  my  summer  research. 

In  order  to  satisfy  the  need  for  online  support  I  designed  a  webpage  dedicated  to  the  toolbox.  This 
website  offers  support  for  his  Graphical  User  Interface  (GUI)  version  of  SOCIT  as  well  as  a  comprehensive 
user  guide.  This  user  guide  describes  each  function  in  depth  and  was  modeled  after  MATLAB’s  help  pages. 
In  addition  to  these  online  aids,  Dr.  Juang’s  toolbox  was  translated  from  MATLAB  to  SCILAB  in  order  to 
increase  its  platform  support.  SCILAB  is  a  matrix-based  scientific  software  package  resembling  MATLAB 
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with  one  very  big  difference.  SCILAB  is  freeware  and  can  be  downloaded  from  the  French  National  Institute 
for  Research  in  Computer  Science  and  Control  SCILAB  Homepage.  By  increasing  the  support  of  Dr.  Juang’s 
toolbox  it  will  become  more  user-friendly  and  this  will  result  in  more  people  using  it. 

Future  applications  of  my  work  will  include  connecting  this  support  webpage  to  the  NASA  Research  site, 
instead  of  its  current  link  from  Dr.  Juang’s  personal  webpage.  Perhaps,  in  the  future  even  more  translations 
of  SOCIT  will  be  created,  thus  furthering  its  software  support. 

Debonding  in  stringer  reinforced  composite  components 

Ronald  Krueger 

Many  composite  components  in  aerospace  structures  are  made  of  flat  or  curved  panels  with  co-cured 
or  adhesively  bonded  frames  and  stiffeners.  Testing  of  stiffened  panels  designed  for  pressurized  aircraft 
fuselage  has  shown  that  bond  failure  at  the  tip  of  the  frame  flange  is  an  important  and  very  likely  failure 
mode.  Comparatively  simple  specimens  consisting  of  a  stringer  flange  bonded  onto  a  skin  were  developed  in 
previous  investigations.  The  failure  that  initiates  at  the  tip  of  the  flange  in  these  specimens  is  identical  to 
the  failure  observed  in  the  full-scale  panels  and  the  frame  pull-off  specimens.  The  objective  of  this  work  is 
to  investigate  the  damage  mechanisms  in  composite  bonded  skin/stringer  constructions  under  uniaxial  and 
biaxial  (in-plane/out-of-plane)  loading  conditions  using  experimental  and  numerical  approaches. 

Based  upon  experimental  findings,  a  two-dimensional  nonlinear  plane-strain  finite  element  (FE)  analysis 
was  performed  using  the  ABAQUS(r)  FE  code.  A  stress  analysis  was  used  to  investigate  the  onset  of 
failure.  This  approach  showed  that  the  location  and  orientation  of  the  initial  transverse  ply  crack  in  the 
flange  are  dependent  on  the  stress  distribution  in  the  critical  area  near  the  flange  tip.  For  all  three  loading 
conditions,  computed  maximum  principal  tensile  stresses  were  almost  identical  and  exceeded  the  transverse 
tension  strength  of  the  material.  A  fracture  mechanics  approach  was  used  to  determine  the  potential  for 
delamination  growth  from  the  initial  transverse  crack.  In  this  approach,  delaminations  of  various  lengths 
originating  from  the  transverse  crack  as  observed  in  the  experiments  were  simulated.  Mode  I  and  mode 
II  strain  energy  release  rate  contributions  were  calculated  for  all  load  cases  using  the  virtual  crack  closure 
technique.  Computed  total  strain  energy  release  rates  were  compared  to  critical  values  obtained  from  an 
existing  mixed-mode  failure  criterion.  The  results  suggest  that  once  a  matrix  ply  crack  has  initiated  in  the 
flange,  a  delamination  will  form  and  grow  in  an  unstable  manner  between  the  adhesive  film  and  the  top  0° 
skin  ply  as  observed  in  the  micrographs. 

Future  work  will  focus  on  the  development  of  a  methodology  for  determining  the  fatigue  life  of  bonded 
composite  skin/stringer  structures  using  delamination  fatigue  characterization  data  and  a  geometric  nonlin¬ 
ear  finite  element  analysis.  Results  will  be  compared  to  fatigue  tests  on  stringer  flange/skin  specimens  to 
verify  the  methodology. 

This  work  is  done  in  collaboration  with  Isabelle  L.  Paris  (National  Research  Council),  T.  Kevin  O’Brien 
(Army  Research  Laboratory,  Vehicle  Technology  Directorate  located  at  NASA  Langley)  and  Pierre  J. 
Minguet  (The  Boeing  Company,  Philadelphia). 

Molecular  weight  effects  on  the  mechanical  response  of  a  glassy  thermoplastic  polyimide 

Lee  M.  Nicholson 

The  objective  of  this  work  is  to  detail  and  summarize  the  physical  and  mechanical  testing  of  an  advanced 
polymer  (LaRC™-SI).  Using  five  known  variations  in  molecular  weight,  the  results  of  mechanical  tests  were 
compared  over  a  range  of  temperatures  below  the  glass  transition.  Results  from  these  tests  will  be  presented 
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along  with  descriptions  of  the  material,  test  methods  and  an  interpretation  of  the  polymer’s  mechanical 
response. 

The  molecular  weights,  molecular  weight  distributions,  glass  transition  temperatures  and  the  mechan¬ 
ical  response  of  an  advanced  thermoplastic,  glassy  polyimide  have  been  determined.  The  notched  tensile 
strength  is  a  strong  function  of  both  temperature  and  molecular  weight.  A  critical  molecular  weight  (Mc) 
is  observed  to  occur  at  a  ~MW  ~  22000  g/mol  below  which,  the  notched  tensile  strength  decreases  rapidly. 
This  critical  molecular  weight  transition  is  not  temperature  dependent.  The  observed  microstructures  and 
nonlinear  (inelastic)  behavior  have  helped  to  characterize  further,  the  brittle  to  ductile  transition  as  a  func¬ 
tion  of  molecular  weight  and  temperature.  Independent  of  temperature,  high  molecular  weight  materials 
(Mw  >  Mc)  will  tend  to  fail  in  a  ductile  manner  and  low  molecular  weight  materials  (Mw  <  Mc)  will  tend 
to  exhibit  brittle  failure.  The  Young’s  Modulus  and  Shear  Modulus  are  strong  functions  of  temperature  only. 
Using  time-temperature  superposition  analysis  of  the  creep  compliance,  low  molecular  weight  materials  have 
increased  creep  compliance  and  increased  creep  rate,  over  that  seen  in  the  high  molecular  weight  materials. 
This  fundamental  understanding  of  how  the  molecular  structure  affects  the  macroscopic  mechanical  perfor¬ 
mance  of  a  polymer  matrix  composite  is  extremely  important  for  the  efficiency  of  design  of  new  materials 
and  resulting  engineered  components.  It  not  only  reduces  material  developmental  costs,  but  also  reduces 
time  to  market  for  components. 

Future  research  areas  are,  to  establish  the  connectivity  between  nano-meso-micro  scale  of  materials 
through  experimentation  and  analysis,  validation  of  computational  tools,  and  establishment  of  key  parame¬ 
ters  that  link  synthesis  and  assembly  to  macro-level  engineering  performance. 

This  research  was  conducted  in  collaboration  with  Karen  S.  Whitley  and  Thomas  S.  Gates  (NASA 
Langley). 

High  temperature  piezoelectric  films 

Zoubeida  Ounaies 

Developing  high  performance,  high  temperature  flexible  piezoelectric  materials  for  active  flow  control 
sensors  and  airframe  health  monitoring  sensors  is  a  critical  component  of  NASA’s  Morphing  Program,  which 
is  focused  on  investigating  smart  and  biomimetic  material  technologies  that  will  enable  self-adaptive  flight 
with  improved  performance  and  safety  air  and  space  travel. 

Early  molecular  modeling  by  NASA  scientists  indicated  the  potential  of  forcing  piezoelectricity  in  amor¬ 
phous  polyimides  by  adding  dipoles  to  the  polymer  structure.  This  was  the  impetus  for  the  current  study  in 
which  we  have  synthesized  a  series  of  polyimides  to  investigate  the  potential  piezoelectric  response  resulting 
from  incorporating  various  dipoles  into  the  structure. 

Structure-property  Studies :  Structure-property  relationship  in  piezoelectric  polyimides  shed  some  light 
into  the  effects  of  dipole  nature,  concentration,  location,  chain  flexibility  and  glass  transition  temperature 
on  piezoelectricity  in  polyimides. 

Thermal  Stability :  We  have  shown  that  the  piezoelectricity  of  this  polymer  was  stable  under  both 
dynamic  and  static  thermal  stimuli  and  statistically  no  loss  was  observed  up  to  100°C.  Since  this  amorphous 
piezoelectric  polyimide  can  generate  a  piezoelectric  response  at  elevated  temperatures,  its  thermal  stability 
should  be  beneficial  for  high  temperature  aerospace  applications. 

Piezocomposite  Processing :  We  propose  to  combine  the  polyimides  with  piezoelectric  powder,  namely 
lead  zirconate  titanate  (PZT)  to  make  a  0-3  piezocomposite.  Such  a  composite  would  exhibit  the  piezoelectric 
coefficients  of  the  ceramic  and  the  flexibility,  strength  and  lightweight  of  the  polymer.  The  relative  ease  of 
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processing  and  potential  high  temperature  use  make  this  material  an  attractive  candidate  for  various  sensing 
applications. 

The  feasibility  of  developing  new  monomers  and  polymers  with  larger  dipole  concentrations  while  main¬ 
taining  chain  flexibility  and  favorable  geometries  is  still  pursued  in  an  attempt  to  increase  the  piezoelectric 
behavior  and  the  use  temperature.  Future  work  will  further  explore  the  processing  and  properties  of  the 
polyimide/PZT  composite. 

This  work  is  done  in  collaboration  with  Joycelyn  S.  Harrison  (NASA  Langley)  and  Cheol  Park  (ICASE). 

Piezoelectric  ceramics  for  use  as  actuators 

Zoubeida  Ounaies 

The  goal  is  to  implement  piezoelectric  actuators  in  vibration  and  noise  suppression  applications  steered 
towards  improving  aircraft  operation.  This  will  lead  to  reduction  of  control  electronics  size  and  cost,  and  a 
lowering  of  the  voltage  drive  by  using  thinner,  more  efficient,  better  characterized  actuators. 

We  have  completed  the  determination  of  capacitive  behavior  and  power  consumption  of  piezoelectric 
actuators.  This  is  an  important  step  in  order  to  optimize  the  system  configuration  and  to  provide  efficient 
driving  electronics  when  using  the  piezoelectric  wafers  as  actuators.  Empirical  relations  were  developed  from 
experimental  data  to  predict  the  capacitance  and  loss  tangent  of  a  piezoelectric  ceramic  as  nonlinear  functions 
of  both  applied  peak  voltage  and  driving  frequency.  It  was  demonstrated  that  by  incorporating  the  variation 
of  capacitance  and  power  loss  with  voltage  and  frequency,  satisfactory  estimates  of  power  requirements  can 
be  obtained.  These  relations  allow  general  guidelines  in  selection  and  application  of  piezoelectric  actuators 
and  driving  electronics  for  active  control  applications. 

To  gain  a  fundamental  understanding  of  how  processing  variables  affect  actuator  performance,  we  have 
recently  completed  displacement  measurements  as  a  function  of  temperature  (25  to  200°  C),  voltage  and 
frequency  on  pre-stressed  piezoelectric  ceramic  transducers.  Results  show  that  the  pre-stressed  ceramics  can 
have  large  displacements  stable  up  to  these  high  temperatures. 

Future  work  includes  understanding  the  effect  of  pre-stressing  inherent  to  the  material  on  nonlinear 
behavior  as  a  function  of  applied  field  and  frequency.  Towards  that  end,  we  will  determine  the  power 
efficiency,  the  electro-mechanical  conversion  coefficient  as  well  as  the  resonance  behavior  as  a  function  of 
geometry. 

We  also  plan  on  continuing  the  modeling  of  the  hysteresis  behavior  of  piezoelectric  ceramics  as  a  func¬ 
tion  of  frequency  and  temperature.  Modeling  the  doming  behavior  and  the  displacement  characteristics  of 
THUNDER  actuators  is  done  in  collaboration  with  Dr.  R.C.  Smith  and  Mr.  R.E.  Wieman. 

This  research  was  conducted  in  collaboration  with  Ralph  Smith  (North  Carolina  State  University)  and 
Thomas  Jordan  (NASA  Langley). 

Molecular  dynamics  simulations  for  mechanical  and  electrical  properties  of  metallic  and 
polymer  systems 

Sun  Mok  Paik 

For  a  quantitative  understanding  of  the  underlying  mechanisms  for  mechanical  and  electrical  failures, 
it  is  necessary  to  investigate  and  quantify  the  effects  of  atomistic  structures  and  defects  on  the  materials 
properties.  The  Molecular  Dynamics  (MD)  simulation  technique  is  an  exact  (within  numerical  accuracy  of 
the  computer)  deterministic  numerical  simulation  method,  and  one  of  the  best-suited  methods  to  study  local 
structures  of  a  system.  In  this  investigation,  we  will  study  the  mechanical  properties  of  metallic  systems,  and 
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the  dielectric  properties  of  polymers  using  the  Molecular  Dynamics  simulation  technique  providing  failure 
sensing  and  preventing  mechanisms. 

The  materials  simulations  fall  into  three  different  categories:  (1)  deterministic  simulations  such  as  the 
MD  simulation,  (2)  stochastic  such  as  a  Monte  Carlo  (MC)  simulation,  and  (3)  ab  initio  calculation.  These 
methods  have  both  pros  and  cons  but  the  MD  method  is  best  suited  for  our  purpose.  Using  the  MD 
technique,  we  will  simulate  a  stable  structure  of  metallic  alloy  systems,  and  the  effects  of  the  structure 
factors,  like  alloy  composition  ratio,  segregation,  grain  boundaries,  impurities,  defects,  etc.,  on  the  crack 
formation  and  propagation.  A  many-body  Embedded  Atomic  Method  (EAM)  potential  is  used  for  the 
interatomic  interactions  between  the  metal  atoms.  We  have  completed  coding  the  MD  program  in  C  for  a 
single  metallic  (Au)  system,  and  are  extending  it  to  the  alloy  systems.  For  polymer  systems,  we  are  in  the 
stage  of  developing  a  suitable  empirical  molecular  interaction  potential. 

In  the  future,  we  will  simulate  and  calculate  various  mechanical  and  dielectric  properties  of  more  com¬ 
plicated  metallic  alloys  and  polymer  systems.  We  also  plan  to  write  an  MD  code  for  a  parallel  computer. 

This  work  is  done  in  collaboration  with  Dr.  Min  Namkung  (NASA  Langley). 


Theoretical  study  of  pulsed  eddy  current  decay 

Sun  Mok  Paik 


The  pulsed  eddy  current  testing  method  has  been  widely  used  in  Non-destructive  Evaluation  (NDE)  and 
depth  measurement  of  the  conducting  plates.  The  method  measures  the  change  of  the  field  near  the  metal 
surface,  and  can  predict  the  depth  and  flaw  of  the  plate  quite  accurately.  However,  there  are  two  serious 
problems  in  real  applications:  (1)  the  fact  that  the  result  is  very  sensitive  to  experimental  conditions  such 
as  the  sensor  position  and  direction  has  made  the  experiments  difficult,  and  (2)  the  thickness  of  the  plates 
that  can  be  measured  are  limited  to  less  than  about  6  mm.  In  this  investigation,  we  study  the  eddy  current 
decay  mechanism  theoretically  and  find  a  way  to  circumvent  the  above  mentioned  problems. 

The  second-order  differential  equation  for  vector  potential  A ,  is  derived  from  Maxwell’s  equations, 


a 2  a  .  dA  a2 a  A/r 
a-a  =  -m>o  +  ^  +  +<*(- 


and  solved  for  a  system  of  two-conducting  plates  as  a  function  of  sensor  position,  r,  sensor  height,  z,  and 
the  lift-off,  l  of  the  drive-coil  and  sensor.  We  find  that  the  field,  A ,  decays  as  time  increases  after  applying 
a  voltage  pulse  with  a  decay  constant  r  that  depends  on  the  thickness  of  the  first  conducting  plate.  For 
the  different  lift-off  /,  the  signal,  A,  changes  but  the  decay  behavior  does  not  change.  Thus,  with  a  suitable 
scaling,  the  tilting  and  lift-off  problems  can  be  avoided  in  an  eddy  current  measurement. 

In  the  future  we  will  investigate  a  method  to  increase  the  maximum  thickness  of  the  plate  that  can  be 
measured  using  the  eddy  current  decay  method. 

This  work  is  done  in  collaboration  with  Dr.  Min  Namkung  (NASA  Langley). 


Sensor-actuator  dual  functional  electroactive  polymeric  hybrid  systems 

Ji  Su 

Sensor  and  actuation  functions  are  two  important  concerns  in  smart  structures  and  biomimetic  tech¬ 
nologies.  Piezoelectric  and  electrostrictive  materials  are  two  promising  candidates  for  being  employed  in 
these  areas.  Piezoelectric  and  electrostrictive  polymeric  materials  have  been  drawing  more  and  more  atten¬ 
tion  due  to  their  lightweight,  excellent  processability,  and  electrical/mechanical  toughness  as  well  as  very 
promising  properties  exhibited  in  some  polymer  materials.  Among  them,  poly  (vinylidene  fluoride)  (PVDF) 
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and  its  copolymers  are  widely  used  as  sensor  materials.  Recently,  several  electrostrictive  polymers  have 
been  developed  for  actuation  applications  due  to  their  large  electric-field-induced  strain  and  high  mechanical 
energy  density,  superior  to  PVDF  and  its  copolymers.  For  actuation  applications,  both  the  piezoelectric 
mechanism  and  electrostrictive  mechanism  are  utilizable.  However,  the  piezoelectric  mechanism  needs  a 
remanent  polarization  and  the  electrostrictive  mechanism  does  not  need  a  remanent  polarization,  these  two 
mechanisms  usually  do  not  appear  in  the  same  material.  The  recently  developed  NASA-LaRC  electrostric¬ 
tive  graft  elastomer  offers  not  only  promising  electromechanical  properties,  but  also  excellent  processability. 
This  allows  the  development  of  a  piezoelectric-electrostrictive  polymeric  molecular  hybrid  composite  system 
using  conventional  processing  methods. 

A  piezoelectric  PVDF-based  copolymer-electrostrictive  graft  elastomer  molecular  composite  system  was 
developed  using  the  solution  casting  method.  The  measurements  of  mechanical  modulus,  dielectric  constant, 
and  piezoelectric  strain  coefficient  were  conducted  as  a  function  of  temperature  and  frequency.  The  electric 
field  induced  strain  response  of  the  molecular  composite  films  was  also  tested.  Both  piezoelectric  and 
electrostrictive  contributions  were  observed  in  the  strain  measurement.  The  research  results  offer  a  material 
which  can  be  utilized  for  unified  sensor/actuator  devices  and  structure. 

More  characterization  work  on  the  newly  developed  materials  is  needed  for  understanding  the  mech¬ 
anisms  further  and  developing  prototype  devices  using  their  materials.  The  characterization  work  will  be 
conducted  in  the  near  future. 

This  research  was  conducted  in  collaboration  with  Joycelyn  S.  Harrison  (NASA  Langley)  and  Zoubeida 
Ounaies  (ICASE). 
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Gottlieb,  Sigal,  Chi- Wang  Shu,  and  Eitan  Tadmor:  Strong  stability  preserving  high-order  time  discretization 
methods .  ICASE  Report  No.  2000-15,  (N AS A/CR-2000-2 10093),  April  7,  2000,  28  pages.  Submitted  to 
SIAM  Review. 

In  this  paper  we  review  and  further  develop  a  class  of  strong-stability  preserving  (SSP)  high-order  time 
discretizations  for  semi-discrete  method-of-lines  approximations  of  partial  differential  equations.  Termed 
TVD  (total  variation  diminishing)  time  discretizations  before,  this  class  of  high-order  time  discretization 
methods  preserves  the  strong-stability  properties  of  first-order  Euler  time  stepping  and  has  proved  very 
useful  especially  in  solving  hyperbolic  partial  differential  equations.  The  new  contributions  in  this  paper 
include  the  development  of  optimal  explicit  SSP  linear  Runge-Kutta  methods,  their  application  to  the  strong 
stability  of  coercive  approximations,  a  systematic  study  of  explicit  SSP  multi-step  methods,  and  a  study  of 
the  strong-stability  preserving  property  of  implicit  Runge-Kutta  and  multi-step  methods. 

Munoz,  Cesar,  and  Victor  Carreno:  Aircraft  trajectory  modeling  and  alerting  algorithm  verification.  ICASE 
Report  No.  2000-16,  (NASA/CR-2000-210097),  April  12,  2000,  25  pages.  Submitted  to  the  International 
Conference  on  Theorem  Proving  in  Higher  Order  Logics  (TPHOL2000). 

The  Airborne  Information  for  Lateral  Spacing  (AILS)  program  at  NASA  Langley  Research  Center  aims 
at  giving  pilots  the  information  necessary  to  make  independent  approaches  to  parallel  runways  with  spacing 
down  to  2500  feet  in  Instrument  Meteorological  Conditions.  The  AILS  concept  consists  of  accurate  traffic 
information  visible  at  the  navigation  display  and  an  alerting  algorithm  which  warns  the  crew  when  one  of 
the  aircraft  involved  in  a  parallel  landing  is  diverting  from  the  intended  flight  path.  In  this  paper  we  present 
a  model  of  aircraft  approaches  to  parallel  runways.  Based  on  this  model,  we  analyze  the  alerting  algorithm 
with  the  objective  of  verifying  its  correctness.  The  formalization  is  conducted  in  the  general  verification 
system  PVS. 

Lallemand,  Pierre,  and  Li-Shi  Luo:  Theory  of  the  lattice  Boltzmann  method :  Dispersion ,  dissipation ,  isotropy , 
Galilean  invariance,  and  stability.  ICASE  Report  No.  2000-17,  (NASA/CR-2000-2 10103),  April  19,  2000, 
33  pages.  To  appear  in  Physical  Review  E. 

The  generalized  hydrodynamics  (the  wave  vector  dependence  of  the  transport  coefficients)  of  a  gener¬ 
alized  lattice  Boltzmann  equation  (LBE)  is  studied  in  detail.  The  generalized  lattice  Boltzmann  equation 
is  constructed  in  moment  space  rather  than  in  discrete  velocity  space.  The  generalized  hydrodynamics  of 
the  model  is  obtained  by  solving  the  dispersion  equation  of  the  linearized  LBE  either  analytically  by  using 
perturbation  technique  or  numerically.  The  proposed  LBE  model  has  a  maximum  number  of  adjustable 
parameters  for  the  given  set  of  discrete  velocities.  Generalized  hydrodynamics  characterizes  dispersion,  dis¬ 
sipation  (hyper- viscosities) ,  anisotropy,  and  lack  of  Galilean  invariance  of  the  model,  and  can  be  applied  to 
select  the  values  of  the  adjustable  parameters  which  optimize  the  properties  of  the  model.  The  proposed 
generalized  hydrodynamic  analysis  also  provides  some  insights  into  stability  and  proper  initial  conditions  for 
LBE  simulations.  The  stability  properties  of  some  2D  LBE  models  are  analyzed  and  compared  with  each 
other  in  the  parameter  space  of  the  mean  streaming  velocity  and  the  viscous  relaxation  time.  The  procedure 
described  in  this  work  can  be  applied  to  analyze  other  LBE  models.  As  examples,  LBE  models  with  various 
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interpolation  schemes  are  analyzed.  Numerical  results  on  shear  flow  with  an  initially  discontinuous  veloc¬ 
ity  profile  (shock)  with  or  without  a  constant  streaming  velocity  are  shown  to  demonstrate  the  dispersion 
effects  in  the  LBE  model;  the  results  compare  favorably  with  our  theoretical  analysis.  We  also  show  that 
whereas  linear  analysis  of  the  LBE  evolution  operator  is  equivalent  to  Chapman-Enskog  analysis  in  the  long 
wave-length  limit  (wave  vector  k  =  0),  it  can  also  provide  results  for  large  values  of  k.  Such  results  are 
important  for  the  stability  and  other  hydrodynamic  properties  of  the  LBE  method  and  cannot  be  obtained 
through  Chapman-Enskog  analysis. 

Povitsky,  Alex:  Improving  jet  reactor  configuration  for  production  of  carbon  nanotubes.  ICASE  Report  No. 
2000-18.  (NASA/CR-2000-210108),  April  20,  2000,  25  pages.  To  be  submitted  to  Computers  and  Fluids. 

The  jet  mixing  reactor  has  been  proposed  for  the  industrial  production  of  fullerene  carbon  nanotubes. 
Here  we  study  the  flowfield  of  this  reactor  using  the  SIMPLER  algorithm.  Hot  peripheral  jets  are  used 
to  enhance  heating  of  the  central  jet  by  mixing  with  the  ambiance  of  reactor.  Numerous  configurations  of 
peripheral  jets  with  various  number  of  jets,  distance  between  nozzles,  angles  between  the  central  jet  and 
a  peripheral  jet,  and  twisted  configuration  of  nozzles  are  considered.  Unlike  the  previous  studies  of  jet 
mixing,  the  optimal  configuration  of  peripheral  jets  produces  strong  non-uniformity  of  the  central  jet  in  a 
cross-section.  The  geometrical  shape  of  reactor  is  designed  to  obtain  a  uniform  temperature  of  a  catalyst. 

Allan,  Brian  G.:  A  reduced  order  model  of  the  linearized  incompressible  Navier-Stokes  equations  for  the 
sens  or /actuator  placement  problem.  ICASE  Report  No.  2000-19,  (NASA/CR-2000-210110),  April  20,  2000, 
30  pages.  Presented  at  the  3rd  ASME/JSME  Joint  Fluids  Engineering  Conference,  FEDSM99-6927,  San 
Francisco,  CA,  July  1999. 

A  reduced  order  modeling  approach  of  the  Navier-Stokes  equations  is  presented  for  the  design  of  a 
distributed  optimal  feedback  kernel.  This  approach  is  based  on  a  Krylov  subspace  method  where  significant 
modes  of  the  flow  are  captured  in  the  model.  This  model  is  then  used  in  an  optimal  feedback  control  design 
where  sensing  and  actuation  is  performed  on  the  entire  flow  field.  This  control  design  approach  yields  an 
optimal  feedback  kernel  which  provides  insight  into  the  placement  of  sensors  and  actuators  in  the  flow  field. 
As  an  evaluation  of  this  approach,  a  two-dimensional  shear  layer  and  driven  cavity  flow  are  investigated. 

Hart,  Roger  C.,  R.  Jeffrey  Balia,  and  G.C.  Herring:  Optical  measurement  of  the  speed  of  sound  in  air  over 
the  temperature  range  300-650  K.  ICASE  Report  No.  2000-20,  (NASA/CR-2000-210114),  April  25,  2000,  9 
pages.  Submitted  to  the  Journal  of  the  Acoustical  Society  of  America. 

Using  laser-induced  thermal  acoustics  (LITA),  the  speed  of  sound  in  room  air  (1  atm)  is  measured 
over  the  temperature  range  300-650  K.  Since  the  LITA  apparatus  maintains  a  fixed  sound  wavelength  as 
temperature  is  varied,  this  temperature  range  simultaneously  corresponds  to  a  sound  frequency  range  10- 
15  MHz.  The  data  are  compared  to  a  published  model  and  typically  agree  within  0. 1-0.4%  at  each  of  21 
temperatures. 
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Roberts,  T.W.,  D.  Sidilkover,  and  S.V.  Tsynkov:  On  the  combined  performance  of  non-local  artificial  bound¬ 
ary  conditions  with  the  new  generation  of  advanced  multigrid  flow  solvers .  ICASE  Report  No.  2000-21, 
(NASA/CR-2000-210115),  April  28,  2000,  43  pages.  To  be  submitted  to  Computers  and  Fluids. 

We  develop  theoretically  and  implement  numerically  a  unified  flow  solution  methodology  that  combines 
the  advantages  relevant  to  two  independent  groups  of  methods  in  CFD  that  have  recently  proven  successful: 
The  new  factorizable  schemes  for  the  equations  of  hydrodynamics  that  facilitate  the  construction  of  optimally 
convergent  multigrid  algorithms,  and  highly  accurate  global  far-field  artificial  boundary  conditions  (ABCs). 
The  primary  result  that  we  have  obtained  is  the  following.  Global  ABCs  do  not  hamper  the  optimal 
(i.e.,  unimprovable)  multigrid  convergence  rate  pertinent  to  the  solver.  At  the  same  time,  contrary  to  the 
standard  local  ABCs,  the  solution  accuracy  provided  by  the  global  ABCs  deteriorates  very  slightly  or  does 
not  deteriorate  at  all  when  the  computational  domain  shrinks,  which  clearly  translates  into  substantial 
savings  of  computer  resources. 

Hart,  Roger  C.,  R.  Jeffrey  Balia,  and  G.C.  Herring:  Simultaneous  velocimetry  and  thermometry  of  air  using 
nonresonant  heterodyned  laser-induced  thermal  acoustics.  ICASE  Report  No.  2000-22,  (NASA/CR-2000- 
210118),  May  5,  2000,  12  pages.  Submitted  to  Applied  Optics. 

Non-resonant  laser-induced  thermal  acoustics  (LITA)  is  employed  with  heterodyne  detection  to  measure 
temperature  (285-295K)  and  a  single  component  of  velocity  (20-150  m/s)  in  an  atmospheric  pressure,  sub¬ 
sonic,  unseeded  air  jet.  Good  agreement  is  found  with  pitot-tube  measurements  of  velocity  (0.2%  at  150  m/s 
and  2%  at  20  m/s)  and  the  isentropic  expansion  model  for  temperature  (0.3%). 

Hysom,  David,  and  Alex  Pothen:  Efficient  parallel  computation  of  ILU(k)  preconditioners.  ICASE  Report 
No.  2000-23.  (NASA/CR-2000-210120),  May  5,  2000,  21  pages.  Proceedings  of  SC’99,  ACM,  November 
1999. 

We  report  the  development  of  a  parallel  algorithm  for  computing  ILU  preconditioners.  The  algorithm 
attains  a  high  degree  of  parallelism  through  employment  of  a  two-level  ordering  strategy,  coupled  with  a 
subdomain  graph  constraint  that  regulates  the  location  of  nonzeros  in  the  Schur  complement.  Experimental 
results  include  timings  on  four  parallel  platforms,  for  problems  with  up  to  20  million  unknowns  running  on 
up  to  216  processors.  The  results  support  our  theoretic  analysis  that  the  algorithm  is  highly  scalable,  for 
both  preconditioner  computation  (factorization)  and  application  (triangular  solve)  stages. 

Hysom,  David,  and  Alex  Pothen:  Parallel  ILU  ordering  and  convergence  relationships:  Numerical  experi¬ 
ments.  ICASE  Report  No.  2000-24,  (NASA/CR-2000-210119),  May  5,  2000,  23  pages.  To  be  submitted  to 
SIAM. 

We  recently  developed  a  parallel  algorithm  for  computing  ILU  preconditioners,  which  was  presented  at 
Super  Computing  1999.  The  algorithm  has  been  shown  to  be  highly  scalable,  in  terms  of  execution  time 
required  for  preconditioner  factorization  and  application,  for  problems  with  up  to  20  million  unknowns  run¬ 
ning  on  up  to  216  processors.  However,  since  the  algorithm  reorders  the  matrix,  and  it  is  widely  known  that 
ordering  can  significantly  affect  convergence,  questions  were  raised  concerning  the  quality  of  the  computed 
preconditioners.  In  this  report  we  present  experimental  results  demonstrating  that  the  orderings  imposed  by 
the  algorithm  do  not  significantly  degrade  convergence,  as  long  as  the  number  of  unknowns  per  subdomain 
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is  not  too  small.  We  report  on  two  model  problems,  Poisson’s  equation,  and  a  special  case  of  the  convection- 
diffusion  equation,  which  other  researchers  have  used  for  ordering  and  convergence  studies.  We  show  that 
convergence  behavior  is  fairly  flat  as  long  as  subdomains  contain  at  least  512  nodes. 

Arian,  E.,  M.  Fahl,  and  E.W.  Sachs:  Trust-region  proper  orthogonal  decomposition  for  flow  control.  ICASE 
Report  No.  2000-25,  (NASA/CR-2000-210124),  May  11,  2000,  23  pages.  To  be  submitted  to  the  Journal  of 
Computational  Physics. 

The  proper  orthogonal  decomposition  (POD)  is  a  model  reduction  technique  for  the  simulation  of  physical 
processes  governed  by  partial  differential  equations,  e.g.  fluid  flows.  It  can  also  be  used  to  develop  reduced 
order  control  models.  Fundamental  is  the  computation  of  POD  basis  functions  that  represent  the  influence 
of  the  control  action  on  the  system  in  order  to  get  a  suitable  control  model.  We  present  an  approach  where 
suitable  reduced  order  models  are  derived  successively  and  give  global  convergence  results. 

Lewis,  Robert  Michael,  Virginia  Torczon,  and  Michael  W.  Trosset:  Direct  search  methods:  Then  and  now. 
ICASE  Report  No.  2000-26,  (NASA/CR-2000-210125),  May  11,  2000,  21  pages.  To  appear  in  a  special  issue 
of  the  Journal  of  Computational  and  Applied  Mathematics. 

We  discuss  direct  search  methods  for  unconstrained  optimization.  We  give  a  modern  perspective  on 
this  classical  family  of  derivative-free  algorithms,  focusing  on  the  development  of  direct  search  methods 
during  their  golden  age  from  1960  to  1971.  We  discuss  how  direct  search  methods  are  characterized  by  the 
absence  of  the  construction  of  a  model  of  the  objective.  We  then  consider  a  number  of  the  classical  direct 
search  methods  and  discuss  what  research  in  the  intervening  years  has  uncovered  about  these  algorithms. 
In  particular,  while  the  original  direct  search  methods  were  consciously  based  on  straightforward  heuristics, 
more  recent  analysis  has  shown  that  in  most  -  but  not  all  -  cases  these  heuristics  actually  suffice  to  ensure 
global  convergence  of  at  least  one  subsequence  of  the  sequence  of  iterates  to  a  first-order  stationary  point  of 
the  objective  function. 

Montero,  Ruben  S.,  and  Ignacio  M.  Llorente:  Robust  multigrid  algorithms  for  the  incompressible  Navier- 
Stokes  equations.  ICASE  Report  No.  2000-27,  (NASA/CR-2000-210126),  May  12,  2000,  24  pages.  To  be 
submitted  to  the  Journal  of  Computational  Physics. 

Anisotropies  occur  naturally  in  CFD  where  the  simulation  of  small  scale  physical  phenomena,  such  as 
boundary  layers  at  high  Reynolds  numbers,  causes  the  grid  to  be  highly  stretched  leading  to  a  slow  down 
in  convergence  of  multigrid  methods.  Several  approaches  aimed  at  making  multigrid  a  robust  solver  have 
been  proposed  and  analyzed  in  literature  using  the  scalar  diffusion  equation.  However,  they  have  been  rarely 
applied  to  solving  more  complicated  models,  like  the  incompressible  Navier-Stokes  equations.  This  paper 
contains  the  first  published  numerical  results  of  the  behavior  of  two  popular  robust  multigrid  approaches 
(alternating-plane  smoothers  combined  with  standard  coarsening  and  plane  implicit  smoothers  combined 
with  semi-coarsening)  for  solving  the  3-D  incompressible  Navier-Stokes  equations  in  the  simulation  of  the 
driven  cavity  and  a  boundary  layer  over  a  flat  plate  on  a  stretched  grid.  The  discrete  operator  is  obtained 
using  a  staggered-grid  arrangement  of  variables  with  a  finite  volume  technique  and  second-order  accuracy  is 
achieved  using  defect  correction  within  the  multigrid  cycle.  Grid  size,  grid  stretching  and  Reynolds  number 
are  the  factors  considered  in  evaluating  the  robustness  of  the  multigrid  methods.  Both  approaches  yield 
large  increases  in  convergence  rates  over  cell-implicit  smoothers  on  stretched  grids.  The  combination  of 
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plane  implicit  smoothers  and  semi-coarsening  was  found  to  be  fully  robust  in  the  flat  plate  simulation  up 
to  Reynolds  numbers  106  and  the  best  alternative  in  the  driven  cavity  simulation  for  Reynolds  numbers 
above  103.  The  alternating-plane  approach  exhibits  a  better  behavior  for  lower  Reynolds  numbers  (below  to 
103)  in  the  driven  cavity  simulation.  A  parallel  variant  of  the  smoother,  tri-plane  ordering,  presents  a  good 
trade-off  between  convergence  and  parallel  properties. 

Liittgen,  Gerald,  and  Michael  Mendler:  The  intuitionism  behind  Statecharts  steps.  ICASE  Report  No. 
2000-28.  (NASA/CR-2000-210302),  July  31,  2000,  45  pages.  To  appear  in  the  27th  International  Collo¬ 
quium  on  Automata,  Languages,  and  Programming  (ICALP  2000). 

The  semantics  of  Statecharts  macro  steps,  as  introduced  by  Pnueli  and  Shalev,  lacks  compositionality. 
This  report  first  analyzes  the  compositionality  problem  and  traces  it  back  to  the  invalidity  of  the  Law  of  the 
Excluded  Middle.  It  then  characterizes  the  semantics  via  a  particular  class  of  linear,  intuitionistic  Kripke 
models,  namely  stabilization  sequences.  This  yields,  for  the  first  time  in  the  literature,  a  simple  fully-abstract 
semantics  which  interprets  Pnueli  and  Shalev’s  concept  of  failure  naturally.  The  results  not  only  give  insight 
into  the  semantic  subtleties  of  Statecharts,  but  also  provide  a  basis  for  an  implementation,  for  developing 
algebraic  theories  for  macro  steps,  and  for  comparing  different  Statecharts  variants. 

Llorente,  Ignacio  M.,  Manuel  Prieto-Matias,  and  Boris  Diskin:  An  efficient  parallel  multigrid  solver  for  3-d 
convection- dominated  problems.  ICASE  Report  No.  2000-29,  (NASA/CR-2000-210319),  August  24,  2000, 
33  pages.  To  be  submitted  to  the  Parallel  Computing  Journal. 

Multigrid  algorithms  are  known  to  be  highly  efficient  in  solving  systems  of  elliptic  equations.  However, 
standard  multigrid  algorithms  fail  to  achieve  optimal  grid-independent  convergence  rates  in  solving  non- 
elliptic  problems.  In  many  practical  cases,  the  non-elliptic  part  of  a  problem  is  represented  by  the  convection 
operator.  Downstream  marching,  when  it  is  viable,  is  the  simplest  and  most  efficient  way  to  solve  this 
operator.  However,  in  a  parallel  setting,  the  sequential  nature  of  marching  degrades  the  efficiency  of  the 
algorithm.  The  aim  of  this  report  is  to  present,  evaluate  and  analyze  an  alternative  highly  parallel  multigrid 
method  for  3-D  convection-dominated  problems.  This  method  employs  semicoarsening,  a  four-color  plane- 
implicit  smoother,  and  discretization  rules  allowing  the  same  cross-characteristic  interactions  on  all  the  grids 
involved  to  be  maintained.  The  resulting  multigrid  solver  exhibits  a  fast  grid-independent  convergence  rate 
for  solving  the  convection-diffusion  operator  on  cell-centered  grids  with  stretching.  The  load  imbalance  below 
the  critical  level  is  the  main  source  of  inefficiency  in  its  parallel  implementation.  A  hybrid  smoother  that 
degrades  the  convergence  properties  of  the  method  but  improves  its  granularity  has  been  found  to  be  the 
best  choice  in  a  parallel  setting.  The  numerical  and  parallel  properties  of  the  multigrid  algorithm  with  the 
four-color  and  hybrid  smoothers  are  studied  on  SGI  Origin  2000  and  Cray  T3E  systems. 

Thomas,  James  L.,  Boris  Diskin,  Achi  Brandt,  and  Jerry  C.  South,  Jr.:  General  framework  for  achieving 
textbook  multigrid  efficiency:  One- dimensional  Euler  example.  ICASE  Report  No.  2000-30,  (NASA/CR- 
2000-210320),  August  24,  2000,  21  pages.  Submitted  to  Frontiers  of  Computational  Fluid  Dynamics  -  2000. 

A  general  multigrid  framework  is  discussed  for  obtaining  textbook  efficiency  to  solutions  of  the  com¬ 
pressible  Euler  and  Navier-Stokes  equations  in  conservation  law  form.  The  general  methodology  relies  on  a 
distributed  relaxation  procedure  to  reduce  errors  in  regular  (smoothly  varying)  flow  regions;  separate  and 
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distinct  treatments  for  each  of  the  factors  (elliptic  and/or  hyperbolic)  are  used  to  attain  optimal  reduc¬ 
tions  of  errors.  Near  boundaries  and  discontinuities  (shocks),  additional  local  relaxations  of  the  conservative 
equations  are  necessary.  Example  calculations  are  made  for  the  quasi-one-dimensional  Euler  equations;  the 
calculations  illustrate  the  general  procedure. 

Allan,  Brian  G.,  Jer-Nan  Juang,  David  L.  Raney,  Avi  Seifert,  Latunia  G.  Pack,  and  Donald  E.  Brown: 
Closed-loop  separation  control  using  oscillatory  flow  excitation.  ICASE  Report  No.  2000-32,  (NASA/CR- 
2000-210324),  August  30,  2000,  27  pages. 

Design  and  implementation  of  a  digital  feedback  controller  for  a  flow  control  experiment  was  performed. 
The  experiment  was  conducted  in  a  cryogenic  pressurized  wind  tunnel  on  a  generic  separated  configuration 
at  a  chord  Reynolds  number  of  16  million  and  a  Mach  number  of  0.25.  The  model  simulates  the  upper 
surface  of  a  20airfoil  at  zero  angle-of-attack.  A  moderate  favorable  pressure  gradient,  up  to  55pressure 
gradient  which  is  relaxed  towards  the  trailing  edge.  The  turbulent  separation  bubble,  behind  the  adverse 
pressure  gradient,  is  then  reduced  by  introducing  oscillatory  flow  excitation  just  upstream  of  the  point  of 
flow  separation.  The  degree  of  reduction  in  the  separation  region  can  be  controlled  by  the  amplitude  of  the 
oscillatory  excitation.  A  feedback  controller  was  designed  to  track  a  given  trajectory  for  the  desired  degree 
of  flow  reattachment  and  to  improve  the  transient  behavior  of  the  flow  system.  Closed-loop  experiments 
demonstrated  that  the  feedback  controller  was  able  to  track  step  input  commands  and  improve  the  transient 
behavior  of  the  open-loop  response. 

Fibich,  Gadi,  and  Semyon  Tsynkov:  High-order  two-way  artificial  boundary  conditions  for  nonlinear  wave 
propagation  with  backscattering.  ICASE  Report  No.  2000-33,  (NASA/CR-2000-210326),  September  6,  2000, 
44  pages.  To  be  submitted  to  the  Journal  of  Computational  Physics. 

When  solving  linear  scattering  problems,  one  typically  first  solves  for  the  impinging  wave  in  the  absence 
of  obstacles.  Then,  using  the  linear  superposition  principle,  the  original  problem  is  reduced  to  one  which 
involves  only  the  scattered  wave  (which  is  driven  by  the  values  of  the  impinging  field  at  the  surface  of 
the  obstacles).  When  the  original  domain  is  unbounded,  special  artificial  boundary  conditions  (ABCs) 
have  to  be  set  at  the  outer  (artificial)  boundary  of  the  finite  computational  domain,  in  order  to  guarantee 
the  reflectionless  propagation  of  waves  through  this  external  artificial  boundary.  The  situation  becomes 
conceptually  different  when  the  propagation  equation  is  nonlinear.  In  this  case  the  impinging  and  scattered 
waves  can  no  longer  be  separated,  and  the  problem  has  to  be  solved  in  its  entirety.  In  particular,  the  boundary 
on  which  the  incoming  field  values  are  prescribed,  should  transmit  the  given  incoming  waves  in  one  direction 
and  simultaneously  be  transparent  to  all  the  outgoing  waves  that  travel  in  the  opposite  direction.  We  call 
this  type  of  boundary  condition  two-way  ABCs.  In  the  paper,  we  construct  the  two-way  ABCs  for  the 
nonlinear  Helmholtz  equation,  which  models  a  continuous-wave  (CW)  laser  beam  propagation  in  a  medium 
with  nonlinear  index  of  refraction.  In  this  case,  the  forward  propagation  of  the  beam  is  accompanied  by 
backscattering,  i.e.,  generation  of  waves  in  the  opposite  direction  to  that  of  the  incoming  signal.  Our  two-way 
ABCs  generate  no  reflection  of  the  backscattered  waves  and  at  the  same  time  impose  the  correct  values  of 
the  incoming  wave.  The  ABCs  are  obtained  in  the  framework  of  a  fourth-order  accurate  discretization  to 
the  Helmholtz  operator  inside  the  computational  domain.  The  fourth-order  convergence  of  our  methodology 
is  corroborated  experimentally  by  solving  linear  model  problems.  We  also  present  solutions  in  the  nonlinear 
case  using  the  two-way  ABC  which,  unlike  the  traditional  Dirichlet  boundary  condition  approach,  allows  for 
direct  calculation  of  the  magnitude  of  backscattering. 


«  -  Al-Theneyan,  Ahmed,  Piyush  Mehrotra,  and  Mohammad  Zubair:  Enhancing  Jini  for  use  across  non¬ 
multicastable  networks .  ICASE  Report  No.  2000-34,  (NASA/CR-2000-210329),  September  20,  2000,  13 
pages.  To  appear  in  Operating  Systems  Review. 

Distributed  heterogeneous  systems  are  being  increasingly  used  to  execute  a  variety  of  large  size  simu¬ 
lation  and  computational  problems.  Resource  management  is  one  of  the  most  important  issues  in  building 
such  systems.  Recently,  Sun  introduced  the  Jini  connection  technology  for  building  plug-and-play  networks 
of  resources.  Jini  relies  on  multicasting  across  the  network  for  its  internal  protocols.  However,  in  a  dis¬ 
tributed  environment,  such  as  the  one  under  consideration  here,  multicasting  may  not  be  supported  across 
the  subnets  for  various  reasons.  In  this  paper  we  describe  enhancements  to  Jini  required  to  use  it  for  building 
a  middleware  resource  management  system  in  a  distributed  environment  that  does  not  support  multicasting. 
In  particular,  we  introduce  a  lightweight  service,  called  the  Tunneling  Service,  which  tunnels  multicast  mes¬ 
sages  across  the  subnets.  We  have  implemented  our  mechanism  and  used  it  to  successfully  tunnel  messages 
between  subnets  in  a  single  domain  and  also  between  different  domains.  In  this  paper  we  describe  our  design 
choices  and  our  implementation  of  the  system. 

Povitsky,  Alex:  Numerical  study  of  wave  propagation  in  a  non-uniform  flow .  ICASE  Report  No.  2000-35, 
(NASA /CR-2000-210533) ,  September  26,  2000,  31  pages.  To  be  submitted  to  the  Journal  of  Computational 
Physics. 

The  propagation  of  acoustic  waves  originating  from  cylindrical  and  spherical  pulses,  in  a  non-uniform 
mean  flow,  and  in  the  presence  of  a  reflecting  wall  is  investigated  by  Hardin  and  Pope  approach  using  compact 
approximation  of  spatial  derivatives.  The  2-D  and  3-D  stagnation  flows  and  a  flow  around  a  cylinder  are 
taken  as  prototypes  of  real  world  flows  with  strong  gradients  of  mean  pressure  and  velocity.  The  intensity  and 
directivity  of  acoustic  wave  patterns  appear  to  be  quite  different  from  the  benchmark  solutions  obtained  in 
a  static  environment  for  the  same  geometry.  The  physical  reasons  for  amplification  and  weakening  of  sound 
are  discussed  in  terms  of  dynamics  of  wave  profile  and  redistribution  of  acoustic  energy  and  its  potential  and 
kinetic  components.  For  an  acoustic  wave  in  the  flow  around  a  cylinder,  the  observed  mean  acoustic  pressure 
is  approximately  doubled  (upstream  pulse  position)  and  halved  (downstream  pulse  position)  in  comparison 
with  the  sound  propagation  in  static  ambient  conditions. 

Dolan,  Elizabeth  D.,  Robert  Michael  Lewis,  and  Virginia  Torczon:  On  the  local  convergence  of  pattern 
search .  ICASE  Report  No.  2000-36,  (NASA/CR-2000-210535),  September  26,  2000,  20  pages.  Submitted 
to  the  SIAM  Journal  on  Optimization. 

We  examine  the  local  convergence  properties  of  pattern  search  methods,  complementing  the  previously 
established  global  convergence  properties  for  this  class  of  algorithms.  We  show  that  the  step-length  control 
parameter  which  appears  in  the  definition  of  pattern  search  algorithms  provides  a  reliable  asymptotic  measure 
of  first-order  stationary.  This  gives  an  analytical  justification  for  a  traditional  stopping  criterion  for  pattern 
search  methods.  Using  this  measure  of  first-order  stationarity,  we  analyze  the  behavior  of  pattern  search 
in  the  neighborhood  of  an  isolated  local  minimizer.  We  show  that  a  recognizable  subsequence  converges 
r-linearly  to  the  minimizer. 


39 


ICASE  COLLOQUIA 


N  ame  /  Affiliat  ion  /  Title _ _ _ 

Walter  Vogler,  Universitaet  Augsburg,  Germany 
“Efficiency  of  Asynchronous  Systems” 

Bruce  Hendrickson,  Sandia  National  Labs 
“Devising  Effective  Parallel  Algorithms” 

Bruce  M.  Boghosian,  Center  for  Computational  Science,  Boston  University 
ICASE  Series  on  Computational  Nanotechnology:  “A  Macroscopic,  Mesoscopic  and 
Microscopic  Look  at  Amphiphilic  Fluids” 

Deepak  Srivastava,  CSC/NAS,  NASA  Ames  Research  Center 
ICASE  Series  on  Computational  Nanotechnology:  “Nanotechnology  of  Molecular 
Electronics  and  Machines:  Carbon  Nanotubes” 

Reto  Luginbuehl,  University  of  Washington,  Seattle 
“Bioengineering  the  Way  to  Smart  Nanomaterials:  Phase  Transitions  at  the  Nanoscale” 

Mario  Rodriguez,  Massachusetts  Institute  of  Technology 
“Identifying  Mode  Confusion  Potential  in  Software  Design” 

Michael  Mendler,  Sheffield  University,  United  Kingdom 
“What  is  in  a  Statecharts  Step,  Exactly?” 

Suvranu  De,  Massachusetts  Institute  of  Technology 
“The  Method  of  Finite  Spheres” 

Michael  Barton,  Analytical  Methods,  Inc.,  Redmond,  WA 
“The  Role  of  Instability  in  Modeling  Wall-bounded  Turbulent  Flows” 

Peter  Habermehl,  University  Paris  7,  Paris,  France 
“On  the  Verification  of  Infinite-state  Systems” 

Sivaram  Arepalli,  G.B.  Tech/Lockheed  Martin 
“Progress  in  Production,  Purification  and  Applications  of  Single  Wall  Carbon 
Nanotubes  at  NAS  A/ JSC” 

Alfons  Geser,  University  of  Tubingen,  Germany 
“More  Correct  Hardware” 


Date 


April  13 


April  18 


April  19 


May  3 


May  5 


May  8 


May  18 


May  26 


June  8 


June  9 


June  13 


June  19 
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N  ame  /  Affiliat  ion/Title 


Date 


Bassam  Bamieh,  University  of  California,  Santa  Barbara 
“Modeling  and  Control  of  Transition  to  Turbulence  in  Wall  Bounded  Shear  Flows” 

Andrew  Adams,  University  of  St.  Andrews,  Scotland 
“Real  Number  Theorem  Proving  in  PVS,  with  Transcendental  Function  Support” 

Ana  Cavalli,  Institut  National  des  Telecomunications,  France 
“Application  of  Formal  Methods  to  Conformance  Test  Generation” 

Raymond  Ristorcelli,  Los  Alamos  National  Laboratory 
“Turbulent  Dispersion  and  Turbulent  Diffusivities  in  the  Presence  of  Stationary  Mean 
Deformations” 

Gilles  Dowek,  INRIA  -  Rocquencourt,  France 
“What  are  Explicit  Proof-terms  Needed  For?” 

Steven  Cox,  University  of  Florida,  Gainesville 
“Global  Optimization  of  the  High  Speed  Civil  Transport” 

Myron  Ginsberg,  HPC  Research  &  Education  and  University  of  Michigan 
“Opportunity  to  Utilize  NASA-supported  Computer  Technology  to  Improve  the 
Performance  for  Large  Industrial  Applications” 

Alex  Povitsky,  ICASE 

“Aeroacoustics  of  a  Non-uniform  Mean  Flow” 

Peter  Bernard,  University  of  Maryland,  College  Park 
“The  Structure  of  Turbulent  Boundary  Layers” 

David  Montgomery,  Dartmouth  College 
“Turbulence  Computations  with  Viscous  Boundary  Conditions” 

Murali  Rangarajan,  University  of  Cincinnati 
“Analysis  of  Requirements  Through  Automated  Application  of  Proof  Obligations” 

Sharath  Girimaji,  Texas  A&M  University 
“Turbulence  Modeling  of  Complex  Flows:  An  Integrated  Approach” 

Yukio  Kaneda,  Nagoya  University,  Japan 

“Some  Results  of  High-resolution  DNS  of  Homogeneous  and  ‘Isotropic’  Turbulence” 


June  20 


June  26 


June  30 


July  5 


July  6 


July  14 


July  26 


July  27 


August  2 


August  2 


August  7 


August  8 


August  9 
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N  ame  /  Affiliat  ion/Tit  le 


Date 


Semyon  Tsynkov,  Tel  Aviv  University,  Israel 
“Non-deteriorating  Numerical  Methods  for  the  Propagation  of  Waves  Over 
Unbounded  Domains” 

Tom  Girardeau,  Georgia  Institute  of  Technology 
“Synthesis  and  Dynamics  of  Poly  (ethylene  glycol) /Cyclodextrin  Rotaxanes” 

Ching  Loh,  Taitech  Inc.  and  NASA  Glenn  Research  Center 
“Computational  Aeroacoustics  by  the  Space-time  CE/SE  Method” 

Eli  Turkel,  Tel  Aviv  University,  Israel 
“Robust  Preconditioning  Methods” 

Saul  Abarbanel,  Tel  Aviv  University,  Israel 
“On  the  Analysis  and  Construction  of  Absorbing  Layers  in  CEM  and  CFD” 

Jason  Rouse,  Lehigh  University 

“Understanding  and  Controlling  the  Swellability  of  Polyelectrolyte/ Clay  Multilayer 
Films” 

Carlo  Bottasso,  Politecnico  di  Milano,  Italy 
“Aerodynamic  and  Structural  Modeling  of  Rotorcraft  Systems” 

Pedro  Camanho,  University  of  Porto,  Portugal 
“A  Progressive  Damage  Model  for  Mechanically  Fastened  Joints  in  Composite 
Laminates” 

Jean-Pierre  Bertoglio,  Ecole  Centrale  de  Lyon,  France 
“Towards  Applications  of  Spectral  Turbulence  Models  to  CFD” 

Peter  Bakker,  Delft  University  of  Technology,  The  Netherlands 
“Topology  of  Three-dimensional  Separations” 

Robert  Yates,  Lawrence  Livermore  National  Laboratory 
“Denotational  Semantics  of  Dataflow  Networks  with  Nonmonotonic  Processes” 

Cord  Rossow,  German  Aerospace  Center 
“On  the  Extension  of  Compressible  Flow  Solvers  to  Incompressible  Flows” 

Tim  Clark,  Los  Alamos  National  Laboratory 
“A  Lattice  Boltzmann  Method  Applied  to  the  Turbulent  Rayleigh-Taylor  Mixing 
Layer” 


August  10 


August  14 


August  15 


August  23 


August  24 


August  25 


August  28 


August  29 


September  5 


September  7 


September  8 


September  13 


September  18 
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Name/Affiliation/Title 


Date 


Gerhard  Goldbeck-Wood,  Polymer  Consortium  Molecular  Simulations,  Inc.,  UK 
“Mesoscale  Simulations  of  Polymers  and  Soft  Materials:  Software  Tools  and 
Applications” 

Kiyosi  Horiuti,  Tokyo  Institute  of  Technology 
“Analysis  of  the  Structures  Responsible  for  Subgrid-scale  Energy  Transfer  in  Large 
Eddy  Simulation” 

Paul  Hovland,  Argonne  National  Laboratory 
“Matrix-free  Newton-Krylov  Methods  Using  a  Hybrid  Automatic  Differentiation- 
finite  Difference  Strategy” 


September  21 


September  26 


September  28 
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ICASE  SUMMER  ACTIVITIES 


The  summer  program  for  2000  included  the  following  visitors: 


VISITOR  and 

AREA  OF  RESEARCH 

AFFILIATION 

DATE  OF 

VISIT 

Abarbanel,  Saul 

Applied  &  Numerical  Math 

Tel  Aviv  University,  Israel 

7/03  -  9/01 

Attar,  Peter  J. 

Structures  and  Materials 

Duke  University 

6/01  -  7/28 

Avihar,  Edan  A. 

Fluid  Mechanics 

Tel  Aviv  University,  Israel 

7/31  -  9/01 

Baggag,  Abdelkader 

Computer  Science 

Purdue  University 

6/05  -  9/08 

Banks,  H.  Thomas 

Applied  &  Numerical  Math 

North  Carolina  State  University 

6/05  -  6/09 

Barton,  J.  Michael 

Fluid  Mechanics 

Analytical  Methods,  Inc. 

6/05  -  6/09 

Bernard,  Peter  S. 

Fluid  Mechanics 

University  of  Maryland 

8/01  -  8/04 

Bertoglio,  Jean-Pierre  E. 

Fluid  Mechanics 

Universite  Claude  Bernard  -  Lyon  1, 

France 

8/28  -  9/08 

Bottasso,  Carlo  L. 

Applied  &  Numerical  Math 

Politecnico  di  Milano,  Italy 

8/21  -  9/29 

Brandt,  Achi 

Applied  &  Numerical  Math 

The  Weizmann  Institute  of  Science, 

Israel 

7/10  -  7/28 

Camanho,  Pedro  M. 

Structures 

University  of  Portugal 

7/31  -  9/29 
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VISITOR  and 

AREA  OF  RESEARCH 

AFFILIATION 

DATE  OF 

VISIT 

Chapman,  Barbara  M. 

Computer  Science 

University  of  Houston 

5/15-6/09 
7/24  -  8/04 

Chepfer,  Helene 

Atmospheric  Sciences 

University  of  Paris  6,  France 

7/03  -  9/30 

Ciardo,  Gianfranco 

Computer  Science 

The  College  of  William  &  Mary 

5/31  -  6/09 
8/14-8/18 

Clark,  Timothy 

Fluid  Mechanics 

Los  Alamos  National  Laboratory 

9/06  -  9/20 

Cox,  Steven  E. 

Applied  &  Numerical  Math 

University  of  Florida,  Gainesville 

5/15  -  7/28 

Criminale,  William 

Fluid  Mechanics 

University  of  Washington,  Seattle 

6/19-6/30 

d’Humieres,  Dominique  G. 

Computer  Science 

Laboratoire  de  Physique  Statistique, 

France 

7/31  -  9/01 

Ditkowski,  Adi 

Applied  &  Numerical  Math 

Brown  University 

8/07-8/18 

Dowek,  Gilles  A. 

Computer  Science 

INRIA,  France 

6/28  -  7/07 

Dowell,  Earl 

Structures  and  Materials 

Duke  University 

7/10-7/12 

Girimaji,  Sharath  S. 

Fluid  Mechanics 

Texas  A&M  University 

7/03 -.8/18 
9/19-9/21 

Goozee,  Richard 

Fluid  Mechanics 

The  University  of  Queensland,  Australia 

9/25  -  10/06 

Gotoh,  Toshiyuki 

Fluid  Mechanics 

Nagoya  Institute  of  Technology,  Japan 

7/31  -  8/04 
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VISITOR  and 

AREA  OF  RESEARCH 

AFFILIATION 

DATE  OF 

VISIT 

Gottlieb,  David 

Applied  &  Numerical  Math 

Brown  University 

7/10  -  7/21 

Griffith,  D.  Todd 

Structures  and  Materials 

University  of  Kentucky 

5/30  -  8/04 

Grohs,  Joshua  W. 

Structures  and  Materials 

Virginia  Polytechnic  Institute  and 

State  University 

5/15  -  8/16 

Habermehl,  Peter  Franz- Joseph 

Computer  Science 

University  of  Paris  7,  France 

5/30  -  6/  09 

Hafez,  Mohammed 

Applied  &  Numerical  Math 

University  of  California,  Davis 

7/31  -  8/04 

Hattori,  Yuji 

Applied  &  Numerical  Math 

Kyushu  Institute  of  Technology,  Japan 

8/21  -  8/25 

Hesthaven,  Jan 

Applied  &  Numerical  Math 

Brown  University 

7/24  -  7/28 
8/07  -  8/18 

Horiuchi,  Kiyoshi 

Fluid  Mechanics 

Tokyo  Institute  of  Technology,  Japan 

9/18  -  9/29 

Hu,  Fang 

Applied  &  Numerical  Math 

Old  Dominion  University 

5/08  -  6/16 

Huedo-Cuesta,  Eduardo 

Applied  &  Numerical  Math 

Universidad  Complutense,  Spain 

7/03  -  7/28 

Jakatdar,  Amol  S. 

Computer  Science 

Old  Dominion  University 

5/22  -  8/25 

Kandil,  Osama 

Applied  &  Numerical  Math 

Old  Dominion  University 

5/23  -  6/29 
Intermittently 

Kaneda,  Yukio 

Fluid  Mechanics 

Nagoya  University,  Japan 

7/31  -  8/11 
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VISITOR  and 

AREA  OF  RESEARCH 

AFFILIATION 

DATE  OF 

VISIT 

Keyes,  David 

Computer  Science 

Old  Dominion  University 

6/01  -  8/31 
Intermittently 

Kholodar,  Denis  B. 

Structures  and  Materials 

Duke  University 

6/26-8/18 

Kulkarni,  Prashanth  P. 

Computer  Science 

Old  Dominion  University 

5/22-8/25 

Lallemand,  Pierre 

Computer  Science 

Universita  Paris  Sud,  France 

9/25-11/03 

Lee,  Min  Hee 

Materials 

Inha  University,  Korea 

7/03  -  8/30 

Li,  Yawei 

Computer  Science 

University  of  Florida,  Gainesville 

5/22-8/10 

Llorente,  Ignacio  M. 

Applied  &  Numerical  Math 

Universidad  Complutense,  Spain 

7/03  -  7/28 

Loh,  Ken 

Fluid  Mechanics 

Taitech,  Inc. 

8/14-8/18 

Longman,  Richard 

Structures  and  Materials 

Columbia  University 

5/01  -  5/05 
5/15-6/16 
7/03  -  9/01 

Mattsson,  Ken  R. 

Applied  &  Numerical  Math 

The  Aeronautical  Research  Institute 

of  Sweden 

8/21-9/15 

Mei,  Renwei 

Computer  Science 

University  of  Florida,  Gainesville 

7/21  -  8/18 

Mendler,  Michael  V. 

Computer  Science 

Sheffield  University,  United  Kingdom 

5/15  -  5/26 

Milder,  Seth 

Computer  Science 

George  Mason  University 

5/23  -  8/04 
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VISITOR  and 

AREA  OF  RESEARCH 

AFFILIATION 

DATE  OF 

VISIT 

Montero,  Ruben  S. 

Applied  &  Numerical  Math 

Universidad  Complutense,  Spain 

7/03  -  7/28 

Naguib,  Ahmed 

Structures  and  Materials 

Michigan  State  University 

6/05  -  6/16 

Nordstrom,  Jan  M. 

Applied  &  Numerical  Math 

The  Aeronautical  Research  Institute 

of  Sweden 

6/10  -  7/15 

Patil,  Amit  J. 

Computer  Science 

University  of  Houston 

5/15  -  9/01 

Pothen,  Alex 

Computer  Science 

Old  Dominion  University 

6/01  -  8/31 
Intermittently 

Prabhakar,  Achal 

Computer  Science 

University  of  Houston 

5/15  -  9/01 

Prieto-Matias,  Manuel 

Applied  &  Numerical  Math 

Universidad  Complutense,  Spain 

7/03  -  7/28 

Qi,  Dewei 

Fluid  Mechanics 

Western  Michigan  University 

7/10  -  8/25 

Ristorcelli,  J.  Raymond 

Fluid  Mechanics 

Los  Alamos  National  Laboratory 

7/03  -  7/07 

Rossow,  Cord  C. 

Applied  &  Numerical  Math 

German  Aerospace  Center 

8/18-9/15 

Ryaben’kii,  Viktor  S. 

Applied  &  Numerical  Math 

Russian  Academy  of  Sciences 

6/26  -  7/10 

Seifert,  Avraham 

Fluid  Mechanics 

Tel  Aviv  University,  Israel 

6/15  -  9/01 

Shu,  Chi- Wang 

Applied  &  Numerical  Math 

Brown  University 

7/17  -  8/04 
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VISITOR  and 

AREA  OF  RESEARCH 

AFFILIATION 

DATE  OF 

VISIT 

Svard,  L.  Magnus 

Applied  &  Numerical  Math 

The  Aeronautical  Research  Institute 

of  Sweden 

8/21-9/15 

Torczon,  Virginia 

Computer  Science 

The  College  of  William  &  Mary 

5/15  -  5/26 

Touil,  Hatem 

Fluid  Mechanics 

Universite  Claude  Bernard  -  Lyon  1, 

France 

8/28  -  9/08 

Tsynkov,  Semyon 

Applied  &  Numerical  Math 

Tel  Aviv  University,  Israel 

6/26-8/15 

Turkel,  Eli 

Applied  &  Numerical  Math 

Tel  Aviv  University,  Israel 

6/19  -  8/25 

Wang,  Lian-Ping 

Fluid  Mechanics 

University  of  Delaware 

6/05  -  6/23 
7/31  -  8/11 

Warburton,  Timothy  C. 

Applied  &  Numerical  Math 

Brown  University 

5/08  -  5/12 
6/12  -  8/11 

Woodruff,  Stephen  L. 

Fluid  Mechanics 

Florida  State  University 

6/19-6/30 

Xu,  Kun 

Fluid  Mechanics 

The  Hong  Kong  University  of  Science 
and  Technology 

6/05  -  6/30 

Yom-Tov,  Jonathan 

Fluid  Mechanics 

Tel  Aviv  University,  Israel 

7/31  -  9/01 

Zhou,  Ye 

Fluid  Mechanics 

Lawrence  Livermore  National  Laboratory 

9/11  -  9/15 
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OTHER  ACTIVITIES 


On  August  3-4,  2000,  ICASE  sponsored  a  Workshop  on  Kinetic  Methods  for  CFD  and  Parallel  Com¬ 
puting  at  ICASE.  The  purposes  of  the  Workshop  were  (1)  to  review  the  LBE  method,  (2)  to  evaluate  the 
range  of  its  potential  applications  to  CFD,  (3)  to  discuss  the  important  issues  which  remain  to  be  solved  in 
the  future  and  (4)  to  assess  the  future  of  the  method.  A  number  of  experts  on  the  LBE  method  and  CFD 
gave  presentations. 


50 


ICASE  STAFF 


I.  ADMINISTRATIVE 

Manuel  D.  Salas,  Director,  M.S.,  Aeronautics  and  Astronautics,  Polytechnic  Institute  of  Brooklyn,  1970. 
Fluid  Mechanics  and  Numerical  Analysis. 

Linda  T.  Johnson,  Office  and  Financial  Administrator 

Barbara  A.  Cardasis,  Administrative  Secretary 

Etta  M.  Morgan,  Accounting  Supervisor 

Emily  N.  Todd,  Conference  Manager/Executive  Assistant 

Shannon  K.  Verstynen,  Information  Technologist 

Gwendolyn  W.  Wesson,  Contract  Accounting  Clerk 

Shouben  Zhou,  Systems  Manager 

Peter  J.  Kearney,  Coral  Manager  /Assistant  Systems  Manager 

II.  SCIENCE  COUNCIL 

David  Gottlieb,  (Chair)  Professor,  Division  of  Applied  Mathematics,  Brown  University. 

Lee  Beach,  Professor,  Department  of  Physics,  Computer  Science  &  Engineering,  Christopher  Newport  Uni¬ 
versity. 

Jack  Dongarra,  Distinguished  Professor,  Department  of  Computer  Science,  University  of  Tennessee. 

Joseph  E.  Flaherty,  Amos  Eaton  Professor,  Departments  of  Computer  Science  and  Mathematical  Sciences, 
Rensselaer  Polytechnic  Institute. 

Forrester  Johnson,  Aerodynamics  Research,  Boeing  Commercial  Airplane  Group. 

John  C.  Knight,  Professor,  Department  of  Computer  Science,  School  of  Engineering  and  Applied  Science, 
University  of  Virginia. 

Robert  W.  MacCormack,  Professor,  Department  of  Aeronautics  and  Astronautics,  Stanford  University. 

Stanley  G.  Rubin,  Professor,  Department  of  Aerospace  Engineering  and  Engineering  Mechanics,  University 
of  Cincinnati. 

Manuel  D.  Salas,  Director,  ICASE,  NASA  Langley  Research  Center. 
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III.  RESEARCH  FELLOWS 


Dimitri  Mavriplis  -  Ph.D.,  Mechanical  and  Aerospace  Engineering,  Princeton  University,  1988.  Applied  k 
Numerical  Mathematics  [Grid  Techniques  for  Computational  Fluid  Dynamics].  (February  1997  to  August 
2003) 

Piyush  Mehrotra  -  Ph.D.,  Computer  Science,  University  of  Virginia,  1982.  Computer  Science  [Programming 
Languages  for  Multiprocessor  Systems].  (January  1991  to  August  2000) 


IV.  SENIOR  STAFF  SCIENTISTS 

Guowei  He  -  Ph.D.,  Theoretical  and  Applied  Mechanics,  Northwestern  Polytechnic  University,  Xian,  China, 
1991.  Fluid  Mechanics  [Turbulence  Modeling  and  Direct  Numerical  Simulation],  (July  2000  to  June  2003) 

Josip  Loncaric  -  Ph.D.,  Applied  Mathematics,  Harvard  University,  1985.  Applied  k  Numerical  Mathematics 
[Multidisciplinary  Design  Optimization],  (March  1996  to  August  2002) 

Li-Shi  Luo  -  Ph.D.,  Physics,  Georgia  Institute  of  Technology,  1993.  Computer  Science  [Parallel  Algorithms]. 
(November  1996  to  August  2002) 

Zoubeida  Ounaies  -  Ph.D.,  Engineering  Science  and  Mechanics,  The  Pennsylvania  State  University,  1995. 
Structures  k  Materials  [Characterization  of  Advanced  Piezoelectric  Materials].  (March  1999  to  November 
2002) 

Alexander  Povitsky  -  Ph.D.,  Mechanical  Engineering,  Moscow  Institute  of  Steel  and  Alloys  Technology 
(MISA),  Russia,  1988.  Computer  Science  [Parallelization  and  Formulation  of  Higher  Order  Schemes  for 
Aeroacoustics  Noise  Propagation].  (October  1997  to  August  2001) 

Robert  Rubinstein  -  Ph.D.,  Mathematics,  Massachusetts  Institute  of  Technology,  1972.  Fluid  Mechanics 
[Turbulence  Modeling,  Turbulence  Management,  and  Acoustics].  (May  1998  to  October  2000) 

David  Sidilkover  -  Ph.D.,  Applied  Mathematics,  The  Weizmann  Institute” of  Science,  1989.  Applied  k 
Numerical  Mathematics  [Numerical  Analysis  and  Algorithms].  (November  1994  to  May  2000) 

Ji  Su  -  Ph.D.,  Polymer  Materials  Science,  Rutgers  University,  1995.  Structures  k  Materials  [Electroactive 
Materials].  (January  2000  to  November  2000) 


V.  SCIENTIFIC  STAFF 

Brian  G.  Allan  -  Ph.D.,  Mechanical  Engineering,  University  of  California  at  Berkeley,  1996.  Applied  k 
Numerical  Mathematics  [Multidisciplinary  Design  Optimization],  (February  1996  to  November  2001) 

Theodorus  Dingemans  -  Ph.D.,  Organic  Chemistry,  University  of  North  Carolina  at  Chapel  Hill,  1998. 
Structures  k  Materials.  (September  2000  to  August  2002) 
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Boris  Diskin  -  Ph.D.,  Applied  Mathematics,  The  Weizmann  Institute  of  Science,  Israel,  1998.  Teaching 
Assistant,  The  Weizmann  Institute  of  Science,  Israel.  Applied  &  Numerical  Mathematics  [Convergence 
Acceleration].  (July  1998  to  September  2002) 

Roger  C.  Hart  -  Ph.D.,  Physics,  University  of  Tennessee,  1991.  Fluid  Mechanics  [Measurement  Science  and 
Technology].  (December  1998  to  October  2001) 

Luc  Huyse  -  Ph.D.,  Civil  Engineering,  Structures,  University  of  Calgary,  Canada,  1999.  Applied  &  Numerical 
Mathematics  [Managing  Uncertainties  in  Multidisciplinary  Research].  (October  1999  to  August  2001) 

Ronald  Krueger  -  Ph.D.,  Aerospace  Engineering,  University  of  Stuttgart,  Germany,  1996.  Structures  Sz 
Materials  [Analysis  of  Composite  Delamination  of  Structures].  (August  2000  to  August  2002) 

Gerald  Luttgen  -  Ph.D.,  Computer  Science,  University  of  Passau,  Germany,  1998.  Computer  Science  [Formal 
Methods  Research  for  Safety  Critical  Systems].  (October  1998  to  August  2000) 

Cesar  A.  Munoz  -  Ph.D.,  Computer  Science,  University  of  Paris  7,  1997.  Computer  Science  [Formal  Methods 
Research  for  Safety  Critical  Systems].  (May  1999  to  April  2002) 

Lee  M.  Nicholson  -  Ph.D.,  Materials  Science,  University  of  Cambridge,  United  Kingdom,  1997.  Structures 
&  Materials  [Computational  Nanotechnolgy].  (May  2000  to  May  2003) 


VI.  VISITING  SCIENTISTS 

Yeon-Gon  Mo  -  Ph.D.,  Electrical  Engineering,  University  of  Nebraska  at  Lincoln,  1999.  Graduate  Research 
Assistant,  Department  of  Electrical  Engineering,  University  of  Nebraska  at  Lincoln.  Structures  &  Materials 
[Smart  Materials  and  Flow  Control].  (October  1999  to  September  2000) 

Sun  Mok  Paik  -  Ph.D.,  Physics,  University  of  Maryland,  1988.  Assistant  Associate  Professor,  Department  of 
Physics,  Kangwon  National  University,  Korea.  Structures  &  Materials  [Computational  Materials].  (February 
2000  to  January  2001) 

J.  Raymond  Skarda-  Ph.D.,  Mechanical  and  Aerospace  Engineering,  Case  Western  Reserve  University,  1996. 
Fluid  Mechanics  [Gravity  Flows].  (August  2000  -  February  2001) 

Linda  Stals  -  Ph.D.,  Mathematics,  Australian  National  University,  1996.  Post-Doc,  Department  of  Computer 
Science,  Old  Dominion  University.  Computer  Science  [ParallelTmplicit  Multilevel  Algorithms].  (November 
1998  to  October  2001) 

Michael  Wagner  -  Ph.D.,  Mathematical  Programming,  Cornell  University,  2000.  Assistant  Professor,  Depart¬ 
ment  of  Mathematics  &  Statistics,  Old  Dominion  University.  Applied  &  Numerical  Mathematics  [Applied 
Optimization].  (September  2000  to  September  2001) 
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VII.  SHORT-TERM  VISITING  SCIENTISTS 


J.  Michael  Barton  -  Ph.D.,  Engineering  Science,  University  of  Tennessee,  2000.  Director  of  Technology, 
Analytical  Methods  Inc.,  Redmond,  WA.  Physical  Sciences  [Fluid  Mechanics].  (June  2000) 

Peter  S.  Bernard  -  Ph.D.,  Mechanical  Engineering,  University  of  California-Berkeley,  1977.  Professor,  De¬ 
partment  of  Mechanical  Engineering,  University  of  Maryland.  Fluid  Mechanics.  (August  2000) 

Jean-Pierre  Bertoglio  -  Ph.D.,  Mechanics,  Ecole  Centrale  de  Lyon,  France,  1980.  Directeur  de  Recherche 
CNRS,  Fluid  Mechanics  Laboratory,  Ecole  Centrale  de  Lyon,  France.  Fluid  Mechanics.  (August  2000  to 
September  2000) 

Hester  Bijl  -  Ph.D.,  Computational  Fluid  Dynamics,  Delft  University  of  Technology,  The  Netherlands.  As¬ 
sistant  Professor,  Department  of  Aerospace  Engineering,  Delft  University  of  Technology,  The  Netherlands. 
Fluid  Mechanics  [Unsteady  Aerodynamics].  (March  2000  to  August  2000) 

Carlo  L.  Bottasso  -  Ph.D.,  Aerospace  Engineering,  Politecnico  di  Milano,  Italy,  1993.  Assistant  Profes¬ 
sor,  Dipartimento  di  Ingegneria  Aerospaziale,  Politecnico  di  Milano,  Milano,  Italy.  Applied  &  Numerical 
Mathematics.  (August  2000  to  September  2000) 

Pedro  Manuel  Camanho  -  Ph.D.,  Composite  Materials,  Imperial  College  of  Science,  Technology  and  Medicine, 
1999.  Assistant  Professor,  Faculty  of  Engineering,  Department  of  Mechanical  Engineering,  University  of 
Porto,  Portugal.  Structures  &  Materials.  (July  2000  to  September  2000) 

Barbara  Chapman  -  Ph.D.,  Computer  Science,  Queen’s  University  of  Belfast,  Northern  Ireland,  U.K.,  1998. 
Associate  Professor,  Department  of  Computer  Science,  University  of  Houston.  Computer  Science.  (May 
2000  to  August  2000) 

Helene  Chepfer  -  Ph.D.,  Applied  Physics,  University  of  Lille,  France,  1997.  Assistant  Professor,  Laboratoire 
de  Meteorologic  Dynamique,  University  of  Paris  6,  France.  Atmospheric  Sciences.  (July  2000  to  September 
2000) 

Gianfranco  Ciardo  -  Ph.D.,  Computer  Science,  Duke  University,  1989.  Assistant  Professor,  Department  of 
Computer  Science,  The  College  of  William  &  Mary.  Computer  Science  [Formal  Methods].  (May  2000  to 
August  2000) 

William  0.  Criminale  -  Ph.D.,  Aeronautics,  The  Johns  Hopkins  University,  1960.  Professor,  Department  of 
Applied  Mathematics,  University  of  Washington.  Fluid  Mechanics.  (June  2000) 

Gilles  A.  Dowek  -  Ph.D.,  Computer  Science,  University  of  Paris  7, 1991.  Research  Scientist,  Institut  National 
de  Recherche  en  Informatique  et  en  Automatique,  INRIA.  Computer  Science.  (June  2000  to  July  2000) 

Toshiyuki  Gotoh  -  Ph.D.,  Applied  Physics,  Nagoya  University,  Japan,  1985.  Professor,  Department  of 
Systems  Engineering,  Nagoya  Institute  of  Technology,  Japan.  Fluid  Mechanics.  (August  2000) 

Peter  Habermehl  -  Ph.D.,  Computer  Science,  VERIMAG,  University  Joseph  Fourier,  Grenoble,  France,  1998. 
Assistant  Professor,  LIAFA,  University  of  Paris  7,  France.  Computer  Science.  (June  2000) 
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Mohammed  Hafez  -  Ph.D.,  Aerospace  Engineering,  University  of  Southern  California,  1972.  Professor,  De¬ 
partment  of  Mechanical  and  Aeronautical  Engineering,  University  of  California- Davis.  Applied  k  Numerical 
Mathematics.  (August  2000) 

Yuji  Hattori  -  Ph.D.,  Physics,  University  of  Tokyo,  1994.  Associate  Professor,  Department  of  Computer 
Aided  Science,  Kyushu  Institute  of  Technology,  Japan.  Applied  k  Numerical  Mathematics  [Higher  Order 
Numerical  Methods].  (August  2000) 

Kiyoshi  Horiuchi  -  Ph.D.,  Applied  Physics,  University  of  Tokyo,  Japan,  1982.  Associate  Professor,  De¬ 
partment  of  Mechano- Aerospace  Engineering,  Tokyo  Institute  of  Technology.  Fluid  Mechanics.  (September 
2000) 

Yukio  Kaneda  -  Ph.D.,  Physics,  Tokyo  University,  Japan,  1976.  Professor,  Department  of  Computational 
Science  and  Engineering,  Graduate  School  of  Engineering,  Nagoya  University.  Fluid  Mechanics.  (July  2000 
to  August  2000) 

Pierre  Lallemand  -  Ph.D.,  Physics,  Universite  de  Paris,  1966.  Director  of  Research,  Centre  National  de  la 
Recherche  Scientifique,  A.S.C.I.,  Universite  Paris-Sud.  Computer  Science.  (September  2000  to  November 
2000) 

Min  Hee  Lee  -  Ph.D.,  Physics,  Korea  Advanced  Institute  of  Science,  1979.  Professor,  Department  of  Physics, 
Inha  University,  Inchon,  Korea.  Structures  k  Materials.  (July  2000  to  August  2000) 

Ignacio  M.  Llorente  -  Ph.D.,  Computer  Science,  Complutense  University  of  Madrid,  Spain,  1995.  Professor, 
Department  of  Computer  Architecture,  Complutense  University  of  Madrid,  Spain.  Applied  &  Numerical 
Mathematics.  (July  2000) 

Ken  Loh  -  Ph.D.,  Applied  Mathematics,  University  of  Western  Ontario,  Canada,  1986.  Research  Scientist, 
Taitech  Inc.,  NASA  John  H.  Glenn  Research  Center.  Fluid  Mechanics.  (August  2000) 

Richard  W.  Longman  -  Ph.D.,  Aerospace  Engineering,  University  of  California- San  Diego,  1969.  Professor, 
Department  of  Mechanical  Engineering,  Columbia  University.  Structures  k  Materials  [Controls].  (May  2000 
to  September  2000) 

Renwei  Mei  -  Ph.D.,  Theoretical  and  Applied  Mechanics,  University  of  Illinois-Urbana,  1990.  Associate 
Professor,  Department  of  Aerospace  Engineering,  Mechanics  k  Engineering  Science  (AeMES),  University  of 
Florida.  Computer  Science.  (July  2000  to  August  2000) 

Michael  Mendler  -  Ph.D.,  Computer  Science,  University  of  Edinburgh,  United  Kingdom,  1993.  Associate 
Professor,  Department  of  Computer  Science,  University  of  Sheffield,  United  Kingdom.  Computer  Science 
[Formal  Methods].  (May  2000) 

Ahmed  M.  Naguib  -  Ph.D.,  Mechanical  and  Aerospace  Engineering,  Illinois  Institute  of  Technology,  1992. 
Assistant  Professor,  Department  of  Mechanical  Engineering,  Michigan  State  University.  Structures  k  Ma¬ 
terials  [Controls].  (June  2000) 
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Dewei  Qi  -  Ph.D.,  Physics,  University  of  Waterloo,  Ontario,  Canada,  1992.  Assistant  Professor,  Department 
of  Paper  &  Printing  Science  &  Engineering,  Western  Michigan  University.  Fluid  Mechanics.  (July  2000  to 
August  2000) 

Cord- Christian  Rossow  -  Ph.D.,  Aerospace  Engineering,  Technical  University  of  Braunschweig,  Germany, 
1988.  Branch  Head,  Dr.-Ing,  DLR,  Institute  of  Design  Aerodynamics,  Germany.  Applied  &  Numerical 
Mathematics.  (August  2000  to  September  2000) 

Avi  Seifert  -  Ph.D.,  Fluid  Mechanics,  Tel-Aviv  University,  1990.  Senior  Lecturer,  Department  of  Fluid 
Mechanics  &  Heat  Transfer,  Faculty  of  Engineering,  Tel  Aviv  University.  Fluid  Mechanics  [Flow  Control 
Development].  (June  2000  to  August  2000) 

Walter  Vogler  -  Ph.D.,  Mathematics,  University  of  Hamburg,  Germany,  1984.  Professor,  Institut  fuer  Infor- 
matik,  Universitaet  Augsburg,  Germany.  Computer  Science  [Formal  Methods].  (April  2000) 

Lian-Ping  Wang  -  Ph.D.,  Mechanical  Engineering,  Washington  State  University,  1990.  Assistant  Professor, 
Department  of  Mechanical  Engineering,  University  of  Delaware.  Fluid  Mechanics.  (June  2000  to  August 
2000) 

Tim  Warburton  -  Ph.D.,  Applied  Mathematics,  Brown  University,  1998.  Visiting  Research  Associate,  Divi¬ 
sion  of  Applied  Math,  Brown  University.  Applied  &  Numerical  Mathematics.  (May  2000  and  June  2000  to 
August  2000) 

Stephen  L.  Woodruff  -  Ph.D.,  Aerospace  Engineering,  The  University  of  Michigan,  1986.  Associate  Professor, 
School  of  Computational  Science  &  Information  Technology,  Florida  State  University.  Fluid  Mechanics. 
(June  2000) 

Kun  Xu  -  Ph.D.,  Astrophysics,  Columbia  University,  1993.  Assistant  Professor,  Department  of  Mathematics, 
The  Hong  Kong  University  of  Science  and  Technology,  Hong  Kong.  Fluid  Mechanics.  (June  2000) 


VIII.  ASSOCIATE  RESEARCH  FELLOW 

David  E.  Keyes  -  Ph.D.,  Applied  Mathematics,  Harvard  University,  1984.  Computer  Science  [Parallel  Nu¬ 
merical  Algorithms] 


IX.  CONSULTANTS 

Saul  Abarbanel  -  Ph.D.,  Theoretical  Aerodynamics,  Massachusetts  Institute  of  Technology,  1959.  Professor, 
Department  of  Applied  Mathematics,  Tel  Aviv  University,  Israel.  Applied  &  Numerical  Mathematics  [Global 
Boundary  Conditions  for  Aerodynamics  and  Aeroacoustic  Computations] 

H.  Thomas  Banks  -  Ph.D.,  Applied  Mathematics,  Purdue  University,  1967.  Professor,  Department  of  Math¬ 
ematics,  Center  for  Research  in  Scientific  Computations,  North  Carolina  State  University.  Applied  fe  Nu¬ 
merical  Mathematics  [Control  Theory] 
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Richard  W.  Barnwell  -  Ph.D.,  Engineering  Mechanics,  Virginia  Polytechnic  Institute  and  State  University, 
1968.  Professor,  Department  of  Aerospace  and  Ocean  Engineering,  Engineering  Science  and  Mechanics. 
Virginia  Polytechnic  Institute  and  State  University.  Fluid  Mechanics  [Turbulence  Modeling] 

Oktay  Baysal  -  Ph.D.,  Mechanical  Engineering,  Louisiana  State  University,  1982.  Eminent  Scholar  and  Pro¬ 
fessor,  Department  of  Aerospace  Engineering,  Old  Dominion  University.  Applied  &  Numerical  Mathematics 

Achi  Brandt  -  Ph.D.,  Mathematics,  The  Weizmann  Institute  of  Science,  1965.  Professor,  Department  of 
Applied  Mathematics,  The  Weizmann  Institute  of  Science,  Israel.  Applied  &  Numerical  Mathematics  [Con¬ 
vergence  Acceleration] 

Thomas  W.  Crockett  -  B.S.,  Mathematics,  The  College  of  William  &  Mary,  1977.  Senior  Research  Associate, 
Computational  Science  Cluster,  The  College  of  William  &  Mary.  Computer  Science  [Scientific  Visualization] 

Ayodeji  0.  Demuren  -  Ph.D.,  Mechanical  Engineering,  Imperial  College  London,  United  Kingdom,  1979. 
Associate  Professor,  Department  of  Mechanical  Engineering  and  Mechanics,  Old  Dominion  University.  Fluid 
Mechanics  [Numerical  Modeling  of  Turbulent  Flows] 

Earl  H.  Dowell  -  Ph.D.,  Aeronautics  &  Astronautics,  Massachusetts  Institute  of  Technology,  1964.  J.A.  Jones 
Professor,  Department  of  Mechanical  Engineering,  Duke  University.  Structures  &  Materials  [Autonomous 
System  Identification  and  Control]  . 

Sharath  Girimaji  -  Ph.D.,  Mechanical  Engineering,  Cornell  University,  1990.  Associate  Professor,  Depart¬ 
ment  of  Aerospace  Engineering,  Texas  A&M  University.  Fluid  Mechanics 

David  Gottlieb  -  Ph.D.,  Numerical  Analysis,  Tel-Aviv  University,  Israel,  1972.  Ford  Foundation  Professor  & 
Chair,  Division  of  Applied  Mathematics,  Brown  University.  Applied  &  Numerical  Mathematics  [Boundary 
Conditions  for  Hyperbolic  Systems] 

Jan  S.  Hesthaven  -  Ph.D.,  Applied  Mathematics/Numerical  Analysis,  Technical  University  of  Denmark,  1995. 
Visiting  Assistant  Professor,  Division  of  Applied  Mathematics,  Brown  University.  Applied  &  Numerical 
Mathematics  [Computational  Electromagnetics] 

Fang  Q.  Hu  -  Ph.D.,  Applied  Mathematics,  Florida  State  University,  1990.  Assistant  Professor,  Department 
of  Mathematics  and  Statistics,  Old  Dominion  University.  Fluid  Mechanics  [Aeroacoustics] 

Osama  A.  Kandil  -  Ph.D.,  Engineering  Mechanics,  Virginia  Polytechnic  Institute  and  State  University,  1974. 
Professor  and  Eminent  Scholar  of  Aerospace  Engineering  and  Chair,  Department  of  Aerospace  Engineering, 
Old  Dominion  University.  Applied  &  Numerical  Analysis  [Computational  Fluid  Dynamics] 

Frank  Kozusko  -  Ph.D.,  Computational  and  Applied  Mathematics,  Old  Dominion  University,  1995.  Assistant 
Professor,  Department  of  Mathematics,  Hampton  University.  Fluid  Mechanics  [Airfoil  Design] 

David  G.  Lasseigne  -  Ph.D.,  Applied  Mathematics,  Northwestern  University,  1985.  Assistant  Professor, 
Department  of  Mathematics  and  Statistics,  Old  Dominion  University.  Fluid  Mechanics  [Asymptotic  and 
Numerical  Methods  for  Computational  Fluid  Dynamics] 
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R.  Michael  Lewis  -  Ph.D.,  Mathematical  Sciences,  Rice  University,  1989.  Assistant  Professor,  Department  of 
Applied  Mathematics,  The  College  of  William  &  Mary.  Applied  &  Numerical  Mathematics  [Multidisciplinary 
Optimization  and  Managing  Uncertainties] 

James  E.  Martin  -  Ph.D.,  Applied  Mathematics,  Brown  University,  1991.  Assistant  Professor,  Department 
of  Mathematics,  Christopher  Newport  University.  Fluid  Mechanics  [Turbulence  and  Computation] 

Jan  Nordstrom  -  Ph.D.,  Numerical  Analysis,  Uppsala  University,  Sweden,  1993.  Senior  Scientist,  The  Aero¬ 
nautical  Research  Institute  of  Sweden.  Applied  &  Numerical  Mathematics  [Global  Boundary  Conditions  for 
Aerodynamic  and  Aeroacoustic  Computations] 

Merrell  L.  Patrick  -  Ph.D.,  Mathematics,  Carnegie  Institute  of  Technology,  1964.  Retired.  Computer  Science 

Alex  Pothen  -  Ph.D.,  Applied  Mathematics,  Cornell  University,  1984.  Professor,  Department  of  Computer 
Science,  Old  Dominion  University.  Computer  Science  [Parallel  Numerical  Algorithms] 

Viktor  Ryaben’kii  -  Ph.D.,  Stability  of  Difference  Equations,  Moscow  State  University,  1953.  Leading  Re¬ 
search  Scientist,  Keldysh  Institute  for  Applied  Mathematics,  Russian  Academy  of  Sciences  and  Full  Professor, 
Department  of  Control  and  Applied  Mathematics,  Moscow  Institute  of  Physics  and  Technology.  Applied  & 
Numerical  Mathematics  [Global  Boundary  Conditions  for  Aerodynamic  and  Aeroacoustic  Computations] 

Chi- Wang  Shu  -  Ph.D.,  Mathematics,  University  of  California-Los  Angeles,  1986.  Associate  Professor,  Divi¬ 
sion  of  Applied  Mathematics,  Brown  University.  Fluid  Mechanics  [Computational  Aeroacoustics] 

Ralph  C.  Smith  -  Ph.D.,  Mathematics,  Montana  State  University,  1990.  Professor,  Department  of  Mathe¬ 
matics,  North  Carolina  State  University.  Applied  &  Numerical  Mathematics  [Optimal  Control  Techniques 
for  Structural  Acoustics  Problems] 

Virginia  Torczon  -  Ph.D.,  Mathematical  Sciences,  Rice  University,  1989.  Assistant  Professor,  Department  of 
Computer  Science,  The  College  of  William  &  Mary.  Computer  Science  [Parallel  Algorithms  for  Optimization 
Including  Multidisciplinary  Optimization] 

Michael  W.  Trosset  -  Ph.D.,  Statistics,  University  of  California-Berkeley,  1983.  Assistant  Professor,  Depart¬ 
ment  of  Mathematics,  The  College  of  William  &  Mary.  Applied  &  Numerical  Mathematics  [Multidisciplinary 
Optimization] 

Semyon  V.  Tsynkov  -  Ph.D.,  Computational  Mathematics,  Keldysh  Institute  for  Applied  Mathematics,  Rus¬ 
sian  Academy  of  Sciences,  1991.  Senior  Lecturer,  Department  of  Applied  Mathematics,  Tel  Aviv  University, 
Israel.  Applied  &  Numerical  Mathematics  [Global  Boundary  Conditions  for  Aerodynamic  and  Aeroacoustic 
Computations] 

Ziaodong  Zhang  -  Ph.D.,  Computer  sciences,  University  of  Colorado-Boulder,  1989.  Professor,  Department 
of  Computer  Science,  The  College  of  William  &  Mary.  Computer  Science 

Mohammad  Zubair  -  Ph.D.,  Computer  Science,  Indian  Institute  of  Technology,  Delhi,  India,  1987.  Professor, 
Department  of  Computer  Science,  Old  Dominion  University.  Computer  Science  [Performance  of  Unstruc¬ 
tured  Flow-solvers  on  Multi-processor  Machines] 
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X.  GRADUATE  STUDENTS 


Ahmed  H.  Al-Theneyan  -  Department  of  Computer  Science,  Old  Dominion  University.  (May  1999  to  Present) 

Peter  Attar  -  Department  of  Engineering  and  Materials  Science,  Duke  University.  (June  2000  to  July  2000) 

Edan  Avihar  -  Department  of  Fluid  Mechanics,  Tel  Aviv  University,  Israel.  (August  2000) 

Orlando  M.  Ayala  -  Department  of  Mechanical  Engineering,  University  of  Delaware.  (August  2000  -  October 

2000) 

Abdelkader  Baggag  -  Department  of  Computer  Science,  Purdue  University.  (September  1995  to  Present) 

Steven  Cox  -  Department  of  Aerospace  Engineering,  Mechanics,  and  Engineering  Science,  University  of 
Florida.  (May  2000  to  July  2000) 

D.  Todd  Griffith  -  Department  of  Mechanical  Engineering,  University  of  Kentucky.  (May  2000  to  August 

2000) 

Gregory  Hicks  -  Department  of  Applied  Mathematics,  Center  for  Research  in  Scientific  Computations,  North 
Carolina  State  University.  (September  2000  to  Present) 

Jianing  Huang  -  Department  of  Computer  Science,  Old.  Dominion  University.  (September  2000  to  Present) 

Eduardo  Huedo-Cuesta  -  Departamento  de  Arquitectura  de  Computadores  y  Automatica,  Universidad  Com- 
..  plutense,  Madrid,  Spain.  (July  2000) 

Amol  Jakatdar  -  Department  of  Computer  Science,  Old  Dominion  University.  (May  2000  to  Present) 

Michelle  Joyner  -  Department  of  Applied  Mathematics,  Center  for  Research  in  Scientific  Computations, 
North  Carolina  State  University.  (September  2000  to  Present) 

Mohamed  Kholief  -  Department  of  Computer  Science,  Old  Dominion  University.  (August  2000  to  Present) 

Denis  B.  Kholodar  -  Department  of  Mechanical  Engineering  and  Material  Science,  Duke  University.  (June 
2000  to  August  2000) 

Hye- Young  Kim  -  Department  of  Aerospace  Engineering,  Texas  A&M  University.  (January  1999  to  August 

2000) 

Brahmadatt  Koodallur  -  Department  of  Computer  Science,  Old  Dominion  University.  (August  2000  to 
Present) 

Prashanth  Kulkarni  -  Department  of  Computer  Science,  Old  Dominion  University.  (May  2000  to  August 

2000) 

Yawei  Li  -  Department  of  Aerospace  Engineering,  Mechanics  and  Engineering  Sciences,  University  of  Florida. 
(January  2000  to  Present) 
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Ken  Mattsson  -  Department  of  Scientific  Computing,  Uppsala  University,  Sweden.  (August  2000  to  Septem¬ 
ber  2000) 

Seth  D.  Milder  -  Department  of  Physics  and  Astronomy,  George  Mason  University.  (September  1997  to 
Present) 

Ruben  Montero  -  Departamento  de  Arquitectura  de  Computadores  y  Automatica,  Universidad  Complutense, 
Madrid,  Spain.  (July  2000) 

Kara  Schumacher  Olson  -  Department  of  Computer  Science,  Old  Dominion  University.  (January  1999  to 
Present) 

Amit  Patil  -  Department  of  Computer  Science,  University  of  Houston.  (May  2000  to  August  2000) 

Juan  A.  Pelaez  -  Department  of  Aerospace  Engineering,  Old  Dominion  University.  (March  1999  to  Present) 

Achal  Prabhakar  -  Department  of  Computer  Science,  University  of  Houston.  (May  2000  to  September  2000) 

Manuel  Prieto-Matias  -  Departamento  de  Arquitectura  de  Computadores  y  Automatica,  Universidad  Com¬ 
plutense,  Madrid,  Spain.  (July  2000) 

Magnus  Svard  -  Department  of  Information  Technology,  Uppsala  University,  Sweden.  (August  2000  to 
September  2000) 

Hatem  Touil  -  Department  of  Fluid  Mechanics,  Ecole  Centrale  de  Lyon,  France.  (August  2000  to  September 

2000) 

Jonathan  Yom-Tov  -  Department  of  Fluid  Mechanics,  Tel  Aviv  University,  Israel.  (August  2000) 

Dazhi  Yu  -  Department  of  Aerospace  Engineering,  Mechanics  and  Engineering  Sciences,  University  of  Florida. 
(January  2000  to  Present) 

XI.  STUDENT  ASSISTANTS 

Joshua  W.  Grohs  -  Department  of  Mechanical  Engingeering,  Virginia  Polytechnic  Institute  and  State  Uni¬ 
versity.  (May  2000  to  August  2000) 
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